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Turbine  Engine  Control  Using  MEMS  for 
Reduction  of  High-Cycle-Fatigue 

(ABSTRACT) 

The  research  effort  presented  in  this  dissertation  consists  of  employing  active  trailing  edge 
blowing  control  to  reduce  the  unsteady  stator-rotor  interaction  in  a  turbofan  simulator.  Two 
active  flow  control  systems  With  different  wake  sensing  approaches  are  successfully 
implemented  on  the  engine  simulator. 

The  first  flow  control  system  utilizes  Pitot  probes  as  flow  sensors.  Use  of  Pitot  probes  as  sensors 
is  appropriate  as  a  first  step  toward  a  more  in  depth  investigation  of  active  trailing  edge  blowing 
control.  An  upper  performance  limit  in  terms  of  wake-filling  can  be  obtained  and  serves  as  the 
baseline  in  evaluating  other  control  systems  with  indirect  wake  sensors.  The  ability  of  the 
system  to  achieve  effective  wake  filling  when  subjected  to  a  change  in  inlet  flow  conditions 
demonstrates  the  feasibility  and  advantage  of  active  flow  control.  Significant  tonal  noise 
reductions  in  the  far  field  are  also  obtained. 

The  second  control  system  involves  using  microphones  as  indirect  wake  sensors.  The 
significance  of  these  acoustic  sensing  approaches  is  to  provide  a  practical  TEB  approach  for 
realistic  engines  implementations.  Microphones  are  flush  mounted  on  the  inlet  case  to  sense  the 
tonal  noise  at  the  blade  passing  frequency.  The  first  sensing  approach  only  uses  the  tone 
magnitude  while  the  second  novel  sensing  approach  utilizes  both  the  tone  magnitude  and  phase 
as  error  information.  The  convergence  rate  of  the  second  sensing  approach  is  comparable  with 
that  of  the  Pitot-prove  based  experiments.  The  acoustic  results  obtained  from  both  sensing 
approaches  agree  well  with  those  obtained  using  Pitot  probes  as  sensors. 

In  addition  to  the  experimental  part  of  this  research,  analytical  studies  are  also  conducted  on  the 
trailing  edge  blowing  modeling  using  an  aeroacoustic  code.  An  analytical  model  for  trailing 
edge  blowing  is  first  proposed.  This  model  is  then  introduced  into  the  two-dimensional 
aeroacoustic  code  to  investigate  effect  of  various  trailing  edge  blowing  managements  in  the  tonal 
sound  generation. 
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Chapter  1 


Introduction 


This  chapter  presents  the  background  information  of  the  research  performed  in  this  dis¬ 
sertation.  First,  high  cycle  fatigue  (HCF)  problem  in  gas  turbine  engines  is  addressed 
in  section  1.1.  It  consists  of  the  description  of  the  HCF  impact  on  aircraft  engines,  the 
HCF  causes,  and  the  current  HCF  reduction  methods.  Trailing  edge  blowing  (TEB), 
known  as  a  wake  management  technique  to  reduce  HCF,  is  discussed  at  great  length.  It 
is  followed  by  the  motivations  and  objectives  of  this  research.  Finally,  the  organization 
of  this  dissertation  is  outlined. 
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1.1  High  Cycle  Fatigue 


Fatigue  is  a  type  of  failure  that  results  in  cracking  or  fracture  when  metal  components 
are  subjected  to  repeated  or  fluctuating  loads.  A  fatigue  failure  starts  when  a  crack 
forms  at  an  initiation  site.  The  crack  leads  to  a  localized  stress  concentration,  which,  if 
high  enough,  will  cause  the  crack  to  propagate  with  each  load  cycle.  As  the  crack  prop¬ 
agates.  the  stress  in  the  component  will  continue  to  rise  until  the  remaining  material 
fails  abruptly.  High  cycle  fatigue  refers  to  fatigue  failures  which  occur  after  the  first 
1000  load  cycles,  while  low  cycle  fatigue  failures  (LCF)  are  regarded  as  those  which 
occur  within  that  number  of  cycles. 


High  cycle  fatigue  is  the  main  cause  of  failure  in  aircraft  engines,  especially  military 
engines.  After  thousands  of  hours  of  flying,  the  engine  blades  may  disintegrate  due  to 
HCF  which  may  cause  an  aircraft  loss  with  absolutely  no  warning.  In  the  past  decades, 
high  cycle  fatigue  has  accounted  for  55%  of  the  US  Air  Force's  engine-related  mishaps 
and  30%;  of  the  total  maintenance  bill.  About  66%  of  US  Air  Force's  maintenance 
actions  are  spent  to  ensure  HCF  does  not  impact  mission  readiness.  The  U.S.  Air 
Force  estimates  an  expenditure  of  about  $100  million  each  year  to  inspect  and  fix  HCF 
related  problems  [1].  In  addition,  HCF  is  also  a  primary  concern  for  future  engines. 
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Figure  1.1  shows  a  typical  turbofan  engine.  The  compressor  of  turbofan  engines  is 
usually  composed  of  multiple  stator-rotor  stages.  The  rotor  blades,  especially  for  the 
first  stage  fan  blades,  are  the  components  which  are  most  vulnerable  to  the  HCF  failure. 
As  will  be  discussed  later  in  detail  in  this  chapter,  one  of  the  causes  of  HCF  failure 
in  the  first  stage  rotor  blades  is  from  the  upstream  stator  wakes,  i.e.,  inlet  guide  vane 
(IGY)  wakes.  Inlet  guide  vane  wake  is  a  flow  velocity  loss  whose  typical  velocity  profile 
is  shown  in  Figure  1.2.  To  reduce  the  HCF  failure  in  the  rotor  blades  by  filling  the 
wake  through  injection  of  air  from  the  trailing  edge  of  the  IGY  is  the  primary  purpose 
of  this  research. 


1.1.1  Causes  of  high  cycle  fatigue  failure 

In  gas  turbine  engines,  the  LCF  failures  occurring  on  the  rotor  blades  are  usually  caused 
by  large  cyclic  excursions  between  zero  stress  and  operating  stress.  For  example,  a 
rotor  disk  would  experience  one  LCF  cycle  as  it  accelerates  from  rest  to  operating 
speed,  then  decelerates  to  shutdown.  On  the  other  hand.  HCF  failures  generally  occur 
due  to  numerous,  smaller  cyclic  excursions  between  two  or  more  stress  states.  These 
excursions  could  result  from  vibrations,  thermal  cycling,  or  other  fluctuating  loads  [2]. 


Figure  1.1:  High  cycle  fatigue  failure  on  the  rotor  blades  of  aircraft 


Figure  1.2:  IGV  wake  velocity  profile. 
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Rotor  blades  are  subject  to  all  of  these  types  of  loads  and  particularly  vulnerable  to 
HCF  failures. 

There  are  two  types  of  rotor  blade  vibrations  which  cause  HCF  failure:  flutter  and 
forced  vibrations.  The  former  is  a  self-excited  oscillation  or  instability  at  or  close  to 
the  natural  frequency  of  the  mechanical  system  which  does  not  require  any  external 
disturbance  to  excite  it  [3,  4],  Since  flutter  is  an  aeroelastic  instability,  it  cannot  be 
stopped  in  most  cases  once  it  is  initiated.  Flutter  measurements  made  about  a  quarter 
century  ago  showed  vibration  amplitudes  of  about  1/8  inch  (0.32  cm)  for  a  blade 
thickness  of  0.15  inch  (0.38  cm).  This  is  precisely  the  reason  why  it  receives  significant 
attention  from  industry.  Forced  vibrations  primarily  arise  from  the  movement  of  the 
rotor  through  flow  defects  or  flow  nonuniformities  [3,  5,  6].  The  three  primary  types  of 
flow  nonuniformities  or  flow  defects  are  inlet  distortion,  wake  disturbances  introduced 
by  the  upstream  IGVs  and  potential  disturbances  from  adjacent  blade  rows.  It  can 
also  arise  from  the  excitation  produced  by  rotor  blades  passing  through  the  rotating 
stall  and  surge,  and  from  the  separated  flow  on  the  rotor  blades.  Large  unsteady 
aerodynamic  loads  can  be  experienced  by  the  blades  when  they  pass  through  these 
flow  defects,  especially  when  the  frequency  of  the  unsteady  loading  is  near  one  of  the 
blade  natural  frequencies.  In  addition  to  the  aerodynamic  sources  described  above, 
forced  vibration  may  also  result  from  mechanical  aspects  such  as  blade  unbalance, 
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mistiming  and  so  forth  [3].  The  lower  modes  of  the  blade  forced  response,  such  as  first 
bending  or  torsion,  may  result  in  fracture  of  an  entire  blade  and  its  attachment.  Such 
a  failure  may  cause  large  unbalance  and  severe  secondary  damage.  The  damage  could 
even  be  greater  if  the  blade  is  not  contained.  At  higher  modes,  such  as  second  torsion 
or  second  or  third  bending,  only  outer  portion  of  a  blade  may  fracture,  causing  some 
unbalance  and  significant  damage.  At  even  higher  modes,  only  a  portion  of  the  tip 
may  break  off.  causing  minimal  damage. 


Inlet  distortion  is  the  physical  distortion  of  inlet  which  results  in  the  nonuniformity  of 
the  inlet  core  flow.  Usually  IGV  wake  disturbances  and  potential  disturbances  are  called 
aerodynamic  gust  forcing  functions  [7.  8,  9],  The  IGY  wake  disturbance  (or  vortices 
disturbance),  a  result  of  viscous  action  at  the  upstream  IGY  surface,  is  generated  by 
velocity  losses  in  the  flow  over  the  downstream  rotor  blades  (see  Figure  1.1).  A  typical 
wake  velocity  profile  of  an  IGY  is  shown  in  Figure  1.2.  The  figure  indicates  that  a 
velocity  deficit  as  compared  to  the  free  stream  occurs  behind  the  upstream  IGY.  The 
IGY  wakes  decay  slowly  with  distance  travelled  [10,  11].  On  the  other  hand,  potential 
disturbances  (sometimes  also  named  as  static  pressure  or  acoustic  disturbance)  are 
an  inviscid  effect  which  is  resulted  from  finite  blade  thickness  and  lifting  properties 
[10,  11].  They  extend  both  upstream  and  downstream,  decaying  exponentially  in  a 
subsonic  flow  (but  do  not  decay  in  a  supersonic  flow).  Therefore,  potential  disturbance 
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Figure  1.3:  Causes  of  HCF  in  the  rotor  blades  of  gas  turbine  engines. 


is  significant  for  closely  spaced  blade  rows  while  the  wake  disturbance  is  dominant  for 


not-so-dosely  spaced  rows. 


A  summary  diagram  for  the  various  causes  of  HCF  is  shown  in  Figure  1.3.  In  this 
dissertation,  the  focus  is  on  the  investigations  of  reducing  the  IG Y  wakes,  to  mitigate 
blade  forced  vibrations,  and  thus  HCF. 
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1.1.2  Current  Research  to  reduce  HCF 


In  the  past,  a  primary  blade  design  consideration  for  bladed-disk  assemblies  was  LCF. 
The  data  base  developed  over  the  past  decades  has  served  well  in  bringing  LCF  related 
failures  under  control.  However,  the  changes  that  have  occurred  during  this  period 
in  the  design  and  development  of  blades  and  disks  has  led  to  failures  due  to  HCF. 
The  problem  looms  large  at  the  present  time  and  the  attempted  fixes  have  been  very 
expensive.  However,  the  current  theoretical  understanding  of  HCF  is  limited  and  design 
for  HCF  tends  to  be  based  on  experience.  Therefore,  there  is  an  overwhelming  demand 
within  the  engine  community  that  a  rational  basis  needs  to  be  developed  to  address 
this  serious  problem. 

Preventing  HCF  failure  of  engine  blades  requires  a  multi-displinary  research  work.  Fig¬ 
ure  1.4  shows  the  basic  categories  of  HCF  research  [1,  2,  3].  It  includes  researches  on 
forced  blade  response,  blade  vibration  damping,  active  control,  test  techniques,  health 
monitoring ,  materials  characterization,  robust  design,  etc.  Forced  blade  response  in¬ 
volves  investigations  of  the  external  forcing  functions  and  blade  structural  dynamics. 
The  external  forcing  functions  include  blade-row  wakes,  potential  field,  and  inlet  dis¬ 
tortion.  It  may  also  include  engine  rotating  stall  and  surge  when  engine  is  working  in 
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an  off-dcsign  condition.  Blade  passive  damping  involves  structural  damping  and  aero¬ 
dynamic  damping  [12].  Structural  damping  refers  to  the  damping  in  the  mechanical 
system  while  aerod;,  namic  damping  results  from  damping  from  the  aero-mechanical 
coupling.  Active  control  means  some  amount  of  external  energy  is  employed  to  reduce 
the  blade  vibration  as  compared  to  the  passive  damping.  Test  techniques  means  the 
measurement  approaches  to  the  aerodynamic  forces,  pressure,  velocity,  blade  vibrations 
when  engine  is  in  an  operating  condition.  Health  monitoring  is  HCF  failure  monitoring 
when  aircraft  is  in  service.  The  purpose  of  materials  characterization  research  is  to 
select  and  improve  blade  material  life.  Robust  Design  involves  blade  design  research 
on  blade  flutter  reduction,  blade  forced  response,  part-to-part  and  engine-to-engine 
variability,  structural  analysis,  and  so  forth. 


It  is  known  that  the  forced  vibratory  response  of  a  system  can  be  reduced  by  three 
possible  approaches.  The  most  effective  way  is  the  alteration  of  the  system  resonant 
frequencies  so  that  they  do  not  match  the  external  forcing  frequencies.  Clearly,  such 
management  strategy  may  include  the  changes  in  blade  geometry  (taper,  shrouds, 
aspect  ratio,  hub  tip  ratio,  etc)  [13.  14].  The  second  approach  is  to  increase  system 
damping  to  reduce  blade  vibrations  by  means  of  designing  dampers  such  as  plat  form  or 
blade-to-blade  dampers.  The  reduction  of  forcing  functions  involves  the  investigation 
of  the  axial  spacing  of  stator-rotor,  blade  and  vane  counts,  inlet  distortion,  potential 
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Figure  1.4:  HCF  research  categories. 


flows,  and  re-energizing  the  stator  wakes  [15,  16]. 


This  dissertation  attempts  to  attenuate  the  forcing  function  by  re-energizing  the  up¬ 
stream  IGY  wakes,  shown  in  Figure  1.2,  using  TEB  technique.  More  specifically,  the 
research  effort  will  be  focused  on  the  active  optimization  of  TEB  mass  flow  rate  for 
best,  wake-filling  using  automatic  control  logic.  The  information  obtained  from  wake 
sensors  is  input  to  a  digital  processor.  According  to  some  control  algorithms,  the  digi¬ 
tal  processor  generates  control  signals  which  are  fed  to  the  flow  actuators  to  adjust  the 
TEB  mass  flow  rate.  The  topic  involving  the  TEB  technique  will  be  elaborated  in  the 


next  section. 
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1.1.3  TEB  technique  to  reduce  HCF 


Supersonic  aircraft  inlets  typically  have  IGVs  which  are  placed  directly  in  front  of  the 
fan  face  and  normally  turn  the  core  flow  into  the  fan  rotor.  As  was  stated  above, 
the  vortex  wakes  are  created  when  the  boundary  layers,  that  are  developed  along  the 
surface  of  the  IGVs  due  to  the  viscous  effect,  meet  at  the  trailing  edge  resulting  in  a 
low  velocity  region  compared  to  the  core  flow.  The  TEB  technique  involves  adding  the 
external  mass  flow  at  the  trailing  edges  to  offset  the  velocity  deficit.  By  optimizing 
the  TEB  rate,  the  flow  into  the  rotor  blades  is  made  more  uniform.  Consequently,  the 
interaction  between  rotor  and  IGV  wakes  is  attenuated,  which  leads  to  the  reduction 
of  forced  blade  vibrations,  and  thus  the  reduction  of  HCF.  The  concept  of  TEB  is 
shown  in  Figure  1.5  which  is  obtained  on  an  F109  turbofan  engine  [17].  It  is  noted 
that  the  y-axis  is  the  flow  velocity  normalized  by  the  free  flow  v^'ocity.  As  we  can 
see,  the  wake  velocity  profile  with  TEB  is  almost  the  same  as  the  core  flow  velocity, 
i.e..  the  wake  is  re-energized.  It,  is  also  worth  to  mention  that  if  we  view  this  problem 
from  the  momentum  perspective,  the  wake  region  constitutes  a  momentum  deficit. 
Trailing  edge  blowing  tries  to  remove  this  momentum  deficit  by  injecting  an  additional 
momentum  so  that  the  wake  region  is  rendered  momentum-deficit-less,  or  briefly  called 
limomentumless  wake’’ .  Therefore,  a  fully  re-energized  wake  is  sometimes  also  called  a 


momentumless  wake. 
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Figure  1.5:  The  wake  profile  before  and  after  TEB. 

The  IGY  wakes  interacting  with  downstream  rotors  are  known  to  contribute  not  only 
to  forced  vibrations  in  rotor  blades,  but  they  generate  tonal  noise  as  well.  This  noise 
source  is  generated  as  the  rotor  blades  chop  through  the  wakes  of  the  upstream  IGYs. 
Hence.  TEB  lias  an  additional  benefit  of  reducing  tonal  fan  noise  through  re-energizing 
the  IGY  wakes. 

Previous  work  on  the  wake  management  using  TEB  is  somewhat  limited.  Probably 
the  earliest  work  on  momentum  injection  from  a  two-dimensional  body  was  done  by 
Wood  [18].  His  work  was  conducted  using  a  standard  airfoil  which  was  equipped  to 
bleed  air  from  its  trailing  edge.  He  discovered  that  the  drag  of  the  airfoil  decreased 
with  an  increased  amount  of  air  flow  from  the  trailing  edge.  He  realized  that  the  drag 
reduction  was  caused  by  the  blowing  flow  to  offset  the  formation  of  the  vorticitv  caused 
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by  the  shear  layers.  However,  he  did  not  address  the  concept  of  moment  uniless  wake 
and  the  concept  of  adding  momentum  to  the  flow. 

Several  experiments  were  performed  recently  to  examine  the  effectiveness  of  different 
wake  management  strategies.  Potential  for  wake  management  using  TEB  to  reduce 
unsteady  blade  loading  has  been  reported  by  Waitz  et.  al  [19].  They  investigated  the 
effect  of  rotor  wake  management  by  removing  and/or  adding  fluid  from  the  flow  around 
the  rotor  blade.  This  study  concluded  that  wake  management  is  feasible  and  that  TEB 
is  more  effective  and  easier  to  implement  than  boundary  layer  suction.  Xaumann  [20], 
Corcoran  [21],  and  Park  and  Cimbala  [22.  23]  had  experimentally  demonstrated  that 
TEB  reduces  significantly  both  the  time-mean  wake  deficit  and  unsteadiness  in  wakes 
shed  by  flat  plates  in  a  water  channel.  Xaumann  showed  a  reduction  of  about  90% 
in  the  time-mean  wake  deficit  and  a  reduction  of  greater  than  50%  in  the  turbulent 
velocity  fluctuations.  Xaumann  and  Corcoran  experimented  with  various  blowing  con¬ 
figurations.  In  particular,  Xaumann  compared  a  continuous  slit  and  a  set  of  discrete 
jets,  with  and  without  vortex  generators.  This  work  showed  that  the  most  effective 
method  of  TEB  was  to  use  a  set  of  discrete  jets,  and  that  the  presence  of  vortex  gen¬ 
erators  enhanced  the  mixing  of  the  wake  with  the  jets.  More  recently.  Brookfield  and 
Waitz  [24]  demonstrated  on  a  1/6-scale  high-bypass-ratio  fan  stage  that  the  TEB  was 
effective  for  reducing  the  rotor  wakes  and  their  mean  harmonic  amplitudes. 
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Previous  experimental  research  at  Virginia  Tech  by  Lcitcli  [25]  and  Saunders  [2G]  has 
demonstrated  the  effectiveness  of  TEB  in  reducing  the  blade  passing  frequency  (BPF) 
tone.  In  the  work  by  Leitch,  different  diameters  for  the  TEB  holes  were  used  to  achieve 
uniform  wake  filling  along  the  span  of  the  strut.  Hence,  optimum  wake  filling  could  only 
be  achieved  at  one  particular  engine  operating  condition.  In  addition,  the  TEB  pressure 
was  adjusted  manually.  Far-field  sound  pressure  reading  showed  noise  reduction  at  the 
BPF  tone  of  up  to  7  dB.  His  research  differed  from  the  current  research  in  that  blowing 
mass  flow  rate  is  optimized  by  wake  control  systems  here. 


In  addition,  few  analytical  studies  on  the  reduction  of  tone  noise  using  TEB  can  are 
found  in  the  open  literature.  Waitz  et  al  [19]  performed  some  calculations  using  the 
computational  codes  to  determine  the  sensitivity  of  the  tone  noise  to  specified  modi¬ 
fication  of  the  wake.  A  two-dimensional  Navier-Stokes  cascade  solver,  UNSFLO  was 
used  to  parameterize  the  effects  of  varying  wake  width  and  deficit  on  unsteady  forces 
on  the  stator  vane.  The  effects  on  the  acoustic  modes  were  determined  using  LINSUB 
(  a  two-dimensional  linearized  panel  method,  using  flatplatc  airfoils).  The  modification 
of  the  wake  in  their  calculation  was  made  by  varying  the  wake  width  and  wake  depth. 
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1.2  Motivations  and  Objectives  of  the  Research 


1.2.1  Motivations 


As  described  above,  TEB  techniques  have  been  successfully  demonstrated  to  be  effec¬ 
tive  in  re-energizing  the  wakes  in  realistic  rotating  environments.  In  addition,  it  is  also 
important  to  mention  that  the  application  of  TEB  for  wake  management  is  feasible 
in  real  engines,  since  the  technology  necessary  to  employ  mass  addition  strategies  is 
mature.  State  of  the  art  high  pressure  turbines  require  blades  with  internal  passages 
and  small  injection  ports  on  the  surface  for  film  cooling.  Turbine  cooling  also  requires 
bleed  air  from  the  compressor,  along  with  tubing  to  duct  the  compressor  air  to  turbine 
stages.  Therefore,  the  application  of  TEB  the  realistic  engine  environments  in  the 
future  is  promising. 

However,  the  realistic  engine  situation  may  require  the  use  of  active  control.  This 
research  is  motivated  by  the  potential  advantages  of  using  an  active  control  system. 


Thev  are: 
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•  Optimum  wake-filling  can  be  achieved  along  the  span  of  IGYs  when  the  flow  rate 
of  each  TEB  hole  is  adjusted  individually  by  using  flow  actuators. 

•  Optimum  wake-filling  can  be  maintained  as  the  engine  working  condition  changes 
due  to  the  tracking  capabilities  of  active  controllers. 

•  The  amount  of  air  for  TEB  can  be  minimized  due  to  the  blowing  rate  optimiza¬ 
tion,  which  is  an  important  factor  for  TEB  technology  to  be  applied  to  actual 
engines. 


To  the  author’s  knowledge,  there  is  no  research  reported  in  the  open  literature  on 
the  active  control  of  TEB  either  experimentally  or  analytically.  The  traditional  active 
engine  control  techniques  are  employed  to  increase  engine  operating  range,  such  as  to 
increase  the  stall  margin  [27.  28].  Here,  active  TEB  control  is  another  novel  application 
in  the  arena  of  active  engine  control.  This  research  is  devoted  to  the  investigation 
of  active  TEB  control  for  the  optimization  of  the  wake-filling  in  a  realistic  rotating 
turbomachinery  environment  to  reduce  HCF  failure  and  fan  noise.  Basically,  an  active 
TEB  control  system  will  involves  developing  wake  sensors,  control  algorithms  and  flow 
actuators.  These  three  control  elements  will  be  the  topics  of  the  following  chapters  in 


this  dissertation. 
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1.2.2  Objectives 


Since  this  is  the  first  research  effort  on  active  TEB  control,  one  of  the  main  objectives 
here  is  to  investigate  the  feasibility  of  implementing  active  flow  control  systems  in 
a  realistic  rotating  environment.  This  includes  the  development  of  each  element  of 
a  flow  control  system:  error  sensors,  flow  actuators  and  control  algorithms.  Another 
objective  also  includes  the  integration  of  these  components  with  a  digital  processor  and 
experimental  demonstrations  of  the  system  validity  on  a  small  scale  engine  simulator. 

The  development  of  practical  wake  sensing  techniques  plays  a  key  role  for  the  imple¬ 
mentation  of  TEB  control  system  to  realistic  situations.  Therefore,  another  objective 
in  this  research  is  to  investigate  and  implement  non-intrusive  wake  sensing  approaches. 
Specifically,  acoustic  approaches  are  going  to  be  pursued  to  sense  indirectly  the  IGY 
wakes. 


Mechanical-Electro-Machine  System  (MEMS)  technology  may  have  considerable  po¬ 
tential  to  affect  and  manipulate  flows  and  thus  to  be  used  in  flow  control  systems.  So 
it  is  also  an  objective  of  this  dissertation  to  integrate  MEMS  based  microvalves  into 
TEB  control  systems  and  evaluate  its  performance  as  flow  actuators. 


Jimvei  Feng 


Chapter  1.  Introduction 


18 


In  addition  to  the  experimental  studies,  active  TEB  is  also  investigated  analytically. 
The  effects  of  mass  flow  rate,  wake-filling  shape  and  axial  spacings  of  IGY-rotor  on 
TEB  and  inlet  noise  are  investigated.  The  results  could  help  to  develop  more  advanced 
acoustic  wake  sensing  strategies. 


1.3  Dissertation  Outline 


The  rest  of  the  dissertation  is  organized  in  the  following  manner. 


Chapter  2  describes  the  experimental  setup  and  research  facilities  used  in  the  TEB 
control  experiments.  A  general  experimental  setup  is  presented  to  introduce  the  basic 
idea  of  this  TEB  control  research.  Research  facilities  include  a  small  scale  turbofan 
propulsion  simulator.  IGVs.  MEMS  based  microvalves,  digital  controllers  and  so  forth. 

Chapter  3  reviews  the  fan  noise  theory.  Basic  knowledge  of  fan  noise  theory  is  required 
for  the  acoustic  measurements  in  the  experiments,  and  for  the  design  of  microphone- 
based  TEB  control  systems  developed  in  the  latter  part  of  this  dissertation. 
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Chapter  4  presents  an  active  flow  control  system  with  Pitot  probes  as  wake  sensors. 
The  probe-based  control  system  with  a  PID  control  algorithm  can  serve  as  the  baseline 
case  in  evaluation  of  other  possible  control  systems  with  indirect  wake  sensors. 


Chapter  5  presents  active  control  systems  using  microphones  flush-mounted  on  the 
inlet  wall  as  the  indirect  wake  sensors.  This  non-intrusive  sensing  technique  provides 
a  more  feasible  approach  to  control  TEB  in  the  realistic  engine  environments. 


In  Chapter  6.  analytical  studies  are  conducted  on  the  active  TEB  control.  The  novelty 
of  the  investigation  is  that  TEB  is  taken  into  account  in  the  noise  prediction  in  the 
inlet.  The  results  from  this  chapter  are  expected  to  provide  theoretical  guidance  to 
design  advanced  acoustic  sensing  approaches. 

Chapter  7  draws  up  the  overall  conclusions  and  outlines  possible  future  work. 


Chapter  2 


Research  Facilities 


This  chapter  describes  the  experimental  setup  and  research  facilities  used  in  the  TEB 
control  experiments.  Firstly,  the  general  experimental  setup  is  described  to  introduce 
the  basic  idea  of  this  TEB  control  research.  Then,  each  component  of  the  setup  is 
elaborated  one  by  one  in  the  following  sections.  A  small  scale  turbofan  propulsion  sim¬ 
ulator  is  introduced  in  section  2.2.  This  simulator  is  used  to  simulate  the  aerodynamic 
and  acoustic  signature  of  the  bypass  fan  in  a  jet  engine.  Section  2.3  describes  an  inlet 
used  with  the  simulator.  Section  2.4  illustrates  IGYs  which  serve  as  the  wake  gener- 
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ators.  It  is  followed  by  the  description  of  MEMS  based  microvalves  and  the  digital 
control  hardware.  Finally,  a  complete  setup  in  the  Virginia  Tech  anechoic  chamber  is 
depicted. 


2.1  The  experimental  setup 

As  stated  in  the  previous  chapter  and  shown  in  Figure  1.1.  the  first  stage  fan  blades 
in  turbofan  engines  are  most  vulnerable  to  the  HCF  failure  due  to  the  unsteady  IG\  - 
rotor  interaction.  In  this  study,  this  first  stage  rotating  environment  is  simulated  by 
employing  a  small  scale  turbofan  engine  simulator.  Since  this  research  is  a  preliminary 
step  toward  a  more  in-depth  investigation  of  using  TEB  control  to  reduce  unsteady 
IGY-rotor  interaction,  a  reduced  number  of  non-turning  IGYs  (four)  is  used  in  the 
inlet  to  serve  as  the  wake  generator.  Xon-turning  IGYs  are  used  because  the  simulator 
is  designed  to  ingest  an  axial  flow  only.  In  addition,  the  non-turning  wakes  generated 
from  the  straight  IGYs  provide  a  simple  bench  for  the  preliminary  TEB  control  study 
as  compared  to  the  use  of  turning  IGYs.  The  first  stage  is  formed  by  the  four  IGYs 
upstream  and  simulator  fan  blades  downstream. 
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Figure  2.1  shows  a  general  schematic  of  the  experimental  setup  used  in  this  research. 
The  main  components  of  the  experiments  are  the  inlet.,  fan  rotor.  IGYs  with  TEB. 
MEMS-based  microvalve  actuators,  error  sensors,  and  the  digital  controller.  As  the 
fan  rotates,  it  draws  the  air  into  the  inlet.  The  flow  velocity  deficit  occurs  behind  the 
four  IGYs  when  the  flow  travels  through  them.  When  the  rotor  chops  the  four  wakes, 
it  experiences  the  unsteady  loading  which  causes  blade  vibration  and  fan  noise.  The 
MEMS-based  microvalves  are  selected  to  serve  as  the  flow  actuators.  A  tube  providing 
each  TEB  hole  with,  the  blowing  air  is  connected  to  a  microvalve.  The  goal  of  these 
microvalves  is  to  adjust  the  TEB  rate  to  achieve  the  best  wake-filling  for  each  IGY. 

This  figure  also  depicts  the  general  idea  of  TEB  control  approach.  The  feedback  control 
loop  begins  with  the  sensing  of  the  wakes.  The  kernel  part  is  the  digital  signal  processor 
(DSP)  which  implements  the  control  algorithms.  The  controller  takes  the  input  from 
the  wake  sensors  and  generates  the  output  according  to  the  control  algorithms.  The 
output  provides  the  control  signal  for  the  microvalves  to  adjust  blowing  rate  to  re¬ 
energizing  the  wakes.  Hence,  a  feedback  loop  is  formed.  This  research  is  focused  on 
investigating  different  wake  sensing  approaches,  developing  of  control  algorithms,  and 
the  implementation  and  validation  of  these  control  concepts.  These  control  approaches 
will  be  described  in  detail  in  the  following  chapters.  The  various  components  of  the 
experimental  setup  are  now  described  in  the  following  sections. 
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Figure  2.1:  The  experimental  setup  for  TEB  control. 

2.2  Turbofan  propulsion  simulator 


The  turbofan  propulsion  simulator  (TPS)  (Model  460)  from  Tech  Development  Inc., 
provides  the  rotor  for  investigating  IGY-rotor  unsteady  interaction  effects.  It  is  repre¬ 
sented  by  the  Fan  Rotor  in  Figure  2.1.  The  TPS  simulates  the  acoustics  and  aerody¬ 
namics  of  a  high  bypass  fan  in  an  aircraft  gas  turbine  engine.  Figure  2.2  is  a  photograph 
of  the  simulator  and  its  rotating  components.  The  TPS  works  in  a  fashion  similar  to 
an  ordinary  gas  turbine,  except  without  a  combustor.  The  power  turbine  is  a  single 
stage  with  29  blades  and  driven  by  high-pressure  air.  The  single  stage  turbine  can 
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Figure  2.2:  Turbofan  propulsion  simulator  (TPS). 

handle  a  maximum  mass  flow  rate  of  1.50  lbm/sec  (0.68kg/sec.)  and  is  connected  on  a 
shaft  to  a  single  stage  rotor.  The  18-blade  rotor  has  a  diameter  of  4.1  inches  (10.4  cm) 
with  a  hub  to  tip  ratio  of  0.44  at  the  fan  inlet.  It  is  capable  of  developing  a  maximum 
pressure  ratio  of  1..6  at  a  mass  flow  rate  of  2.72  lbm/sec  (1.23kg/sec.).  and  maximum 
speed  of  80.000  rpm.  A  row  of  26  Exit  Guide  Vanes  (EGA')  is  present  one  rotor  chord 
length  downstream  of  the  fan.  These  guide  vanes  also  cause  far-field  noise  radiation 
due  to  the  unsteady  lift  generated  by  the  passage  of  rotor  wakes.  Instrumentation  on 
the  TPS  turbine  includes  thermocouples  to  measure  the  turbine  bearing  temperatures 
and  a  magnetic  speed  pickup  to  measure  the  rotational  speed.  The  fan  blade  tip  speed 
becomes  supersonic  at  the  speed  of  65,000  RPM. 
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2.3  Inlet 


Many  different  inlet  configurations  can  be  mounted  on  the  TPS  body.  This  research 
requires  an  inlet  to  be  designed  such  that  it  focuses  on  the  rotor-IGY  interaction  effects. 
In  this  study,  the  symmetric  inlet  design  used  in  previous  research  by  T.  Lcitch  [29]  is 
used  as  shown  in  Figure  2.3. 

Figure  2.3  compares  the  inlet  used  in  this  study  and  representative  supersonic  and 
subsonic  inlets.  The  inlet  is  longer  than  a  conventional  subsonic  inlet  and  equipped 
with  a  centerbodv  supported  by  the  IGYs.  On  the  other  hand,  the  inlet  is  shorter  than 
a  typical  supersonic  inlet.  The  centerbodv  is  greatly  simplified,  and  the  outer  profile 
does  not  accelerate  the  flow  through  the  inlet. 

Figure  2.4  shows  a  three  dimensional  view  of  the  inlet  alone.  The  inlet  is  a  constant 
area  with  a  diameter  of  4.1  inches  (10.4  cm).  A  larger  diameter  step  at  the  discharge 
end  allows  the  inlet  to  be  mounted  on  the  simulator.  The  mouth  of  the  inlet  is  pro¬ 
vided  with  a  bellmouth  to  avoid  inlet  flow  separation.  Five  tapped  holes  are  provided 
circumferentially.  The  centers  of  four  of  these  holes  line  up  with  the  trailing  edge  of 
each  IGV.  and  the  remaining  one  is  positioned  in  the  free  stream  (not  shown  in  Figure 
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Bellmouth 


Figure  2.4:  The  cutaway  view  of  the  inlet. 

2.4).  These  holes  are  used  for  mounting  either  Pitot  probe  flow  sensors  or  microphone 
sensors.  The  centerbody  is  a  cylinder  with  a  conical  section  at  the  inlet.  This  allows 
for  smooth  turning  of  the  flow  entering  the  inlet.  The  armulus  area  between  the  cen¬ 
terbody  and  inside  diameter  of  the  inlet  is  of  constant  cross  section.  Four  slots  are  cut 
on  the  surface  of  the  centerbody,  90  degree  apart.  The  IGYs  are  placed  in  these  slots 
and  hence  a  smooth  geometric  transition  is  obtained.  Upon  assembly,  care  must  be 
taken  to  fill  in  all  cavities,  to  avoid  the  interaction  of  the  rotor  with  flow  distortions 
generated  by  these  cavities. 
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2.4  Blowing  IGVs 


The  four  IGVs  (also  referred  as  stators  in  the  text  where  appropriate),  which  support 
the  centerbody,  also  serve  as  the  wake  generators.  They  introduce  spatially  periodic 
distortions,  i.e.  four  wakes  in  the  flow  field  ingested  by  the  rotor.  The  IGVs  are  faceted, 
with  sharp  leading  and  trailing  edges  and  provide  no  turning  to  the  flow.  The  geometry 
of  these  struts  is  very  similar  to  those  used  by  T.  Leitch  [30],  except  for  the  trailing 
edge  geometry.  The  trailing  edge  used  here  is  sharp  as  compared  to  the  blunt  edge 
used  by  Leitch  [30].  The  sharp  trailing  edge  has  the  effect  of  reducing  the  width.  This 
also  results  in  a  reduction  of  the  far-field  sound  levels  generated  for  the  baseline  case 
with  no  trailing  edge  blowing. 


The  trailing  edge  of  the  IGVs  is  provided  with  six  equally-spaced  blowing  holes.  The 
size  and  spacing  of  the  trailing  edge  blowing  holes  are  critical  to  provide  uniform  span- 
wise  re-energizing  of  the  IGVs  wakes.  The  goal  in  designing  the  trailing  edge  hole 
configuration  is  to  produce  a  blowing  profile  from  the  tip  to  the  hub  that  matches  the 
wake  profile.  A  secondary  consideration  is  that  the  hole  configuration  be  designed  to 
use  a  minimal  amount  of  blowing  air.  After  a  series  of  bench  tests.  Leitch  [25]  con¬ 
cluded  that  six  equally-spaced  holes  along  the  span  provided  the  best  configuration. 
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In  addition,  Lcitch  used  different  hole  diameters  to  achieve  uniform  wake  filling  along 
the  span  of  the  IGY.  In  this  research,  the  trailing  edge  holes  have  the  same  diame¬ 
ters  because  the  supply  air  is  adjusted  for  each  hole  individually  using  MEMS  based 
microvalves.  The  IGY  used  in  the  research  is  shown  in  Figure  2.5  while  Figure  2.6 
gives  the  detailed  dimensions.  The  hole  centers  are  spaced  two  hole  diameters  apart. 
The  blowing  holes  are  1/16  inches  (0.16  cm)  in  diameter.  Each  hole  is  'L'  shaped  and 
opens  out  at  the  top  of  the  IGY  in  a  1/8  inches  (0.32  cm)  diameter  hole.  Because  of 
the  sharp  trailing  edge,  the  blowing  holes  Create  a  crenellation  as  they  open  out  on 
the  trailing  edge.  This  is  expected  to  aid  in  mixing  of  the  discrete  jets  and  the  wake. 
Leitch  [25]  pointed  out  that  reduction  in  tone  was  obtained  in  spite  of  lack  of  adequate 
wake  filling  at  the  hub.  Hence  the  blowing  holes  are  moved  up  so  that  in  the  bottom 
25%  of  IGY  height  the  wake  is  not  completely  re-energized. 


The  IGYs  are  positioned  so  that  the  fan  face  is  an  IGY  chord  length  downstream  of 
the  trailing  edges,  i.e.,  2.088  inches  (5.303  cm).  Bench  tests  showed  that  for  uniform 
pitchwise  mixing  in  the  wake  at  least  0.5  IGY  chord  length  is  required  as  illustrated  in 
Figure  2.7.  This  means  that  if  a  flow  sensor  is  used,  it  has  to  be  placed  no  less  than 
0.5  IGY  chord  length  downstream  of  the  IGY  trailing  edge  to  ensure  that  the  sensor 
would  measure  the  mixed  out  flow  from  the  jet-wake  system,  rather  than  the  jet  issuing 
out  of  the  blowing  holes.  This  is  an  important  consideration  for  the  implementation  of 
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Figure  2.7:  The  jet  mixing  process. 


active  flow  control  using  flow  sensor  as  described  in  Chapter  4. 


2.5  MEMS  technology  and  MEMS  based  microvalve 


Since  one  of  the  objectives  in  this  research  is  to  investigate  the  feasibility  of  integrat¬ 
ing  the  MEMS  technology  into  the  TEB  control  systems,  a  basic  background  review 
of  MEMS  is  introduced  first  in  the  following  subsection.  MEMS  is  a  field  which  is 
rich  in  scientific  challenges  as  well  as  an  enabling  technology  for  broad  engineering 
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applications. 


2.5.1  MEMS  Technology 


MEMS  differ  from  IC’s  (integrated  circuits)  in  that  they  consist  of  both  mechanical  and 
electrical  devices.  While  the  electronic  devices  provide  the  MEMS  with  the  function¬ 
alities  of  IC’s,  the  integration  of  the  mechanical  devices  such  as  sensors  and  actuators 
enable  MEMS  to  interact  with  the  external  environment.  Although  electrical  devices 
were  well  established  technically,  the  study  of  mechanical  devices  such  as  microactua¬ 
tors  began  only  about  two  decades  ago.  In  the  past,  almost  all  the  mechanical  parts 
had  sizes  much  larger  than  the  electronic  elements  because  the  traditional  manufac¬ 
turing  technology  lacked  the  capability  to  fabricate  small  mechanical  parts.  In  1983  u 
new  manufacturing  process,  the  surface  micromachining  technology  [31]  became  avail¬ 
able  which  could  be  used  to  manufacture  the  micron  sized  mechanical  devices  with 
sophisticated  geometry.  Since  then,  a  remarkable  progress  has  been  made  and  some 
MEMS  products  have  been  marketed  [32,  33].  MEMS  technology  is  so  attractive  due 
to  its  three  characteristic  features  or  three  “M”s  [33]:  Miniaturization,  Multiplicity  and 
Microelectronics.  Miniaturization  is  clearly  essential:  Multiplicity  is  the  key  to  sue- 
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ccssful  microsystems  and  Microelectronics  is  critical  for  mechanical  parts  to  cooperate 
with  each  other  and  to  perform  a  given  task.  The  application  areas  of  MEMS  include 
optics,  transportation  and  aerospace,  robotics,  biotechnologies,  medical  engineering, 
etc.  It  is  widely  believed  that  MEMS  technology  will  have  profound  impact  in  the 
future  society  just  like  the  invention  of  IC  chips  forty  years  ago. 

There  are  numerous  opportunities  for  applying  MEMS  to  flow  control.  For  example. 
Ho  et.  al  [34,  35]  demonstrated  the  effectiveness  using  MEMS  to  modify  the  streamwise 
development  for  shear  flow.  MEMS  based  micro  transducers  for  flow  control  can  easily 
increase  the  temporal  response  to  a  large  extent  because  of  the  low  inertia  of  the 
transducer  [36].  Some  other  good  examples  for  MEMS  used  in  the  arena  of  flow  control 
include  MEMS  based  micro  anemometer,  micro  shear  stress  sensor,  micro  pressure 
sensor,  micro  actuator,  etc  [37]. 


2.5.2  MEMS  based  microvalves 


In  this  research,  we  aim  to  integrate  MEMS  technology,  i.e.,  MEMS  based  microvalves, 
into  the  active  TEB  control  systems  as  flow  actuators.  In  the  active  TEB  control 


Jimvei  Feng 


Chapter  2.  Research  Facilities 


34 


Figure  2.8:  The  MEMS-based  microvalve  used  to  control  TEB  rate  (a)  picture  and  (b) 
dimensions. 

systems  developed  in  Chapters  4  and  5.  individual  control  of  the  blowing  rate  of  each 
hole  is  obtain'd  by  the  use  of  MEMS  based  microvalves.  The  microvalves  used  for 
this  research  are  NO-3000  Fluistor,  supplied  by  Redwood  Microsystems  [38.  39].  The 
microvalves  arc  Normally  Open  (NO),  i.e.  if  no  signal  voltage  is  applied  the  microvalves 
are  fully  open  and  operate  at  maximum  flow  rate  as  governed  by  the  source  pressure. 
The  design  specifications  are  given  in  Table  2.1.  Figure  2.8  shows  the  microvalve 
and  the  external  dimensions  and  pin  connections.  Figure  2.9  depicts  the  principle 
of  operation  of  the  microvalves.  The  microvalve  operates  on  the  thermo- pneumatic 
principle.  The  actuator  is  a  sealed  cavity,  which  holds  the  control  liquid  inside  it. 
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Tabic  2.1:  MEMS-based  microvalvc  specifications. 


Fluid  Media 

Non-corrosive  gases 

Outlet  Pressure  Range 

•  0  —  100  psig  (0  —  6.89E5  pa) 

Maximum  Supply  Pressure 

100  psig  (6.89E5  pa) 

Proof  /  Burst  Pressure 

150  /  250  psig  (1.03E6  /  1.72E6  pa) 

Flow  Rate  Range 

0.2  -  3000  seem  (4.1E-12-  G.15E-8  kg/sec) 

Response  Time 

500  ms  typical 

Average  Power  Consumption 

700  m Watts 

Operating  Temperature  Range 

0  -  40  C 

Control  Signal  Voltage  Range 

0  —  15  Volts 

Power  Supply 

15  Volts 

One  wall  of  the  cavity  is  a  diaphragm,  which  flexes  when  the  cavity  pressure  increases 
due  to  evaporation  of  the  liquid.  A  thin  film  resistor  provides  the  heat  required  for 
vaporization.  The  flexing  of  the  diaphragm  wall  is  translated  into  valve  movement. 
Thus,  flow  rate  is  governed  by  varying  the  signal  voltage  that  for  a  Normally  Open 
valve  is  from  0  to  15  volts.  However,  it  is  found  that  the  microvalve  performance 
depends  on  initial  conditions  such  as  fully  open,  fully  close,  or  partially  open.  Figure 
2.10,  which  shows  the  dynamic  pressure  measured  at  the  exit  of  the  microvalve  as  a 
function  of  the  signal  voltage,  is  obtained  when  the  microvalve  starts  from  the  initial 
condition  of  fully-open  [40].  It  is  also  observed  that  the  microvalve  has  a  response  time 
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of  about  10  seconds  if  it  starts  from  an  initial  condition  of  fully  close.  On  the  other 
hand,  the  time  response  is  about  3-4  seconds  if  it  starts  from  partially  or  fully  open 
condition.  This  indicates  that  the  microvales  has  a  very  slow  response.  Hence,  the 
measurements  for  Figure  2.10  were  conducted  when  the  dynamic  pressure  has  reached 
steady  state  condition.  It  can  be  seen  from  Figure  2.10  that  the  microvalve  exhibits  a 
nonlinear  operation  near  the  fully  open  and  fully  closed  regions.  The  figure  also  shows 
that  the  effective  voltage  operating  range  is  supposed  to  be  from  3  to  9  volts. 

Aerodynamic  studies  are  first  conducted  to  investigate  the  effectiveness  of  wake-filling 
obtained  by  controlling  individual  blowing  holes,  before  the  tests  on  the  small  scale 
engine  simulator  are  carried  out  in  the  anechoic  chamber.  This  initial  step  is  important 
because  the  performance  of  MEMS-based  microvalves  as  flow  actuators  needs  to  be 
evaluated.  As  mentioned  in  Chapter  1,  this  is  one  of  the  research  objectives  because 
flow  actuators  play  an  important  role  in  TEB  control  systems. 

The  experiments  are  conducted  on  a  test  bench  by  Rao  [41,  42],  The  aerodynamic 
results  are  presented  in  the  form  of  2-D  contour  plots,  in  Figure  2.11.  The  X-axis 
is  the  pitchwise  direction,  nondimensionalized  by  the  IGY  thickness.  The  Y-axis  is 
the  spanwise  direction,  nondimensionalized  by  the  IGV  height.  The  contours  display 


the  ratio  of  the  axial  velocity  to  the  maximum  free  stream  velocity.  The  choice  of 
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Figure  2.9:  The  MEMS-based  microvalve  working  principle. 


0.0  ’  u  2.:1  j.c  40  so  eo  :.l  an  *.e  icc  n.u  12.U  13c  •*.»; 

Signal  VoHilge  (V) 


Figure  2.10:  The  performance  curve  of  MEMS  based  microvalve 
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maximum  free  stream  velocity  is  made  to  assign  a  common  base  for  all  ratios,  and  also 
to  bring  out  the  variations  that  exist  in  the  flow  field.  Also  seen  in  these  plots  are  two 
vertical  lines  at  0.5,  representing  the  maximum  thickness  of  the  IGY.  The  circles  with 
crosshairs  represent  the  location  of  the  trailing  edge  blowing  holes.  Data  is  acquired 
0.5  IGY  chord  lengths  downstream  of  the  IGV,  with  and  without  blowing. 

Figure  2.11(a)  is  a  velocity  contour  plot  of  the  baseline  case,  i.e.  without  trailing 
edge  blowing.  The  wake  is  centered  about  the  trailing  edge  of  the  strut,  as  seen  by 
the  maximum  velocity  defect  occurring  at  the  X-axis  origin.  Also  noted  from  this 
figure  is  that  the  free  stream  axial  velocity  is  not  uniform  in  the  spanwise  direction. 
Figure  2.11(b)  shows  the  velocity  ratio  contours  with  the  application  of  trailing  edge 
blowing.  The  data  is  acquired  after  the  blowing  rate  of  each  blowing  hole  is  optimized 
by  adjusting  the  control  voltage  to  the  MEMS  based  microvalve.  The  blowing  rate  is 
optimized  when  the  total  pressure  in  the  wake  equals  the  free  stream  total  pressure. 
Figure  2.11(b)  indicates  that  at  most  spanwise  locations,  the  wake  velocity  is  equal  to 
or  varies  in  small  amount  with  the  axial  velocity  in  the  free  stream  at  the  same  span. 
In  addition,  the  spanwise  variation  of  axial  velocity  in  the  wake  is  very  similar  to  that 
in  the  free  stream.  These  results  have  demonstrated  the  advantage  of  individual  hole 
control  by  using  MEMS  based  microvalves  in  the  presence  of  spanwise  variations  in 


free  stream  axial  velocity. 
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Figure  2.11:  Comparison  of  axial  velocity  ratio  (a)  without  TEB,  (b)  with  TEB. 


2.6  Digital  Signal  Processing  Board 


The  control  algorithms  in  this  research  are  implemented  in  software  on  a  digital  sig¬ 
nal  processing  board  from  Spectrum  Signal  Processing  Inc..  The  kernel  part  of  the 
board  is  the  processor  TMSC320C30  manufactured  by  Texas  Instrument  [43].  The 
TMSC320C30  DSP  contains  integer  and  floating-point  arithmetic  units.  2048  x  32  bit 
words  of  on-chip  RAM,  4096  x  32  bit  words  of  on-chip  ROM,  control  unit  with  parallel 
and  serial  interfaces.  Operating  from  a  33.3MHz  clock,  a  performance  of  16.7  million 
instruction  per  second  is  achieved.  The  peak  arithmetic  performance  of  33.3  million 
floating-point  computations  per  second  is  attained  when  the  floating-point  multiplier 
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and  adder  arc  used  in  parallel.  This  board  is  also  equipped  with  dual  channel  1G  bit 
A/D  and  D/A  systems,  a  parallel  expansion  system  referred  to  DSPLIXK  and  a  serial 
expansion  interface.  The  board  is  mounted  in  a  personal  computer  when  it  is  in  use. 
Access  to  memory  passes  through  the  dual  porting  hardware  on  the  board. 


In  our  application  case,  an  external  16-channel  D/A  board  and  an  external  32-channel 
A/D  board  are  utilized  to  work  with  the  DSP  card  through  DSPLINK.  since  the  number 
of  on-board  input/ouput  A/D  and  D/A  channels  is  not  sufficient.  Each  channel  in  the 
A/D  board  consists  of  a  double  buffered  12  bit  DAC.  2nd  order  programmable  analog 
output  filter  and  buffer  amplifier.  The  output  from  the  card  ranges  from  -  8.192  volts 
to  8.192  volts.  The  32  analog  input  channels  in  the  A/D  card  are  multiplexed  to  a  fast 
(3  psec.)  12  bit  ADC.  All  32  channels  have  input  buffering  and  a  first  order  low  pass 
to  reduce  unwanted  high  frequency  noise.  The  channels  are  arranged  in  groups  of  4. 
the  channels  in  each  group  being  sampled  simultaneously  allowing  differential  signal 


acquisition. 
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2.7  Schematic  setup  in  anechoic  chamber 


Figure  2.12  describes  a  schematic  of  the  entire  active  flow  control  setup  for  testing  in 
the  Virginia  Tech  anechoic  chamber.  An  air  tank  with  compressed  air  is  used  to  power 
the  simulator  as  well  as  it  is  connected  to  a  plenum  that  provides  the  air  source  for 
TEB.  The  simulator  and  the  active  flow  control  components  are  placed  in  the  anechoic 
chamber.  The  anechoic  environment  is  required  since  the  acoustic  effect  of  the  TEB  is 
monitored.  The  Virginia  Tech  anechoic  chamber  with  a  dimension  of  13.1  x  8.9  x  6.6 
ft  (3.99  x  2.71  x  2.01meter).  is  constructed  with  Owens  Corning  Type  705  industrial 
fiberglass  wedges.  It  can  be  considered  anechoic  for  frequency  above  200  Hz.  which  is 
far  below  the  blade  passing  frequency  in  the  tests.  The  ambient  noise  pressure  level 
is  about  30  dB  (Ref.  20  x  10~6  Pa).  This  ambient  noise  level  does  not  affect  the 
measured  noise  level  during  testing  since  minimum  in-test  noise  level  is  above  70  dB. 


The  simulator  is  mounted  onto  the  top  of  a  test  stand  which  is  placed  at  the  center  of 
the  chamber.  The  test  stand  is  1.2  meters  (47.2  inches)  in  height  above  the  elevated 
chamber  floor  to  minimizing  the  potential  of  ingestion  of  ground  vortices,  which  would 
add  to  inlet,  distortions.  A  2-inch  (5.08  cm)  diameter  air  line  is  flange  mounted  to  the 
top  plate  to  provide  the  high  pressure  air  to  power  the  simulator.  The  inlet  slides  onto 
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Figure  2.12:  The  schematic  setup  in  Virginia  Tech  anechoic  chamber. 

the  simulator  fan  shroud  and  is  held  in  place  by  a  cowl.  The  cowl  runs  the  length 
of  the  simulator  and  forms  an  exhaust  duct  for  the  turbofan  discharge.  An  orifice 
plate  is  placed  inside  the  cowl  and  located  using  set  screws.  Loading  of  the  turbofan 
is  controlled  by  the  area  of  the  annulus  formed  by  the  orifice  plate  with  the  turbine 
exhaust  duct.  The  turbine  exhaust  is  ducted  separately  from  the  fan  exhaust  so  as  not 
to  affect  the  fan  loading.  Well  downstream  of  the  simulator  exit,  the  exhausts  from  the 
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turbine  and  fan  arc  mixed  in  the  exhaust  pipe.  The  exhaust  pipe  is  12  inches  (30.48 
cm)  in  diameter,  made  from  sheet  metal,  and  exhausts  turbine  and  fan  air  into  the 
atmosphere. 

Compressed  air  is  required  not  only  to  power  the  simulator  but  also  to  supply  blowing 
air.  An  Ingersoll-Rand  four-stage  reciprocating  compressor  meets  the  compressed  air 
requirements.  The  compressed  air  facility  is  equipped  with  an  activated-alumina  drying 
unit  and  a  heat  exchanger.  This  ensures  that  the  compressed  air  that  drives  the  turbine 
is  free  of  water  vapor  and  also  at  room  temperature.  The  compressed  air  is  carried  by 
a  main  line  that  is  2  inches  (5.08  cm)  in  diameter.  It  is  fitted  with  a  solenoid  valve  that 
serves  as  an  emergency  shut-off.  An  1-inch  (2.54  cm)  line  from  the  main  supply  line 
taps  off  air  required  for  blowing.  The  air  is  routed  to  a  supply  plenum  placed  in  the 
chamber.  An  air  regulator  mounted  on  the  feed  line  maintains  plenum  pressure.  Four 
pipes  connect  the  supply  plenum  to  a  manifold.  The  manifold  is  a  block  of  aluminum 
with  four  1/4-inches  (0.635  cm)  holes  drilled  into  it.  Each  hole  has  six  tapped  holes. 
Microvalves  are  mounted  in  these  holes  using  Swagelok  fittings.  Each  microvalve  is 
connected  to  a  single  blowing  hole  on  an  IGY  through  Tygon  tubing.  Care  is  taken  to 
ensure  that  each  joint  is  leak  proof. 

A  Toshiba  process  loop  controller  Model  EC321  [44]  in  conjunction  with  Peek  Mea- 
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Figure  2.13:  The  simulator  speed  control  system. 

surement  STD  6000  I/P  (current  to  pressure)  transducer  [45]  and  ACT-2  A  simulator 
speed  readout  is  used  for  simulator  speed  control.  The  control  configuration  is  shown 
in  Figure  2.13.  Basically  the  process  controller  EC321  is  a  PID  (proportional-integral- 
derivative)  feedback  controller  with  programmable  functions.  The  STD  6000  I/P  trans¬ 
ducer  serves  as  the  actuator  for  EC321  controller  while  the  simulator  speed  readout 
ACT-2.A  n  be  viewed  as  the  control  system  sensor.  This  readout  uses  the  pulses 
from  the  6  magnetic  speed  pickups  that  are  built  in  the  simulator  as  its  input  after 
they  are  amplified  by  a  gain  of  10.  Reference  speed  is  set  by  the  user  using  internal 


programmable  functions. 


Chapter  3 


Turbofan  Noise  Theory 


A  review  of  fan  noise  theory  is  required  for  the  design  and  implementation  of  TEB 
control  systems  using  the  acoustic  sensing  approaches  in  the  latter  part  of  this  disser¬ 
tation.  In  section  3.1,  the  relationship  between  fan  noise  and  blade  vibration  in  terms 
of  the  common  generation  mechanism  is  introduced.  The  basics  of  fan  noise  theory  in¬ 
cludes  the  generation  of  tone  frequencies,  excited  acoustic  duct  modes  and  propagation 
behavior  in  a  circular  duct. 
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3.1  Introduction 


Noise  from  turbofan  and  jet  engines  is  one  of  the  most  significant  noise  sources  of 
military  and  commercial  aircraft.  The  acoustic  signature  of  the  aircraft  engines  is 
characterized  by  a  broadband  noise  level  with  dominant  discrete  tones  at  blade  passing 
frequency  (BPF)  and  its  harmonics.  Among  the  discrete  tones,  the  noise  level  at  BPF 
tone  is  usually  the  loudest.  For  subsonic  fans,  these  discrete  tones  are  usually  10-20 
dB  above  the  broadband  level.  For  transonic  fans,  additional  tones  also  appear  in  the 
inlet  at  frequencies  submultiples  of  the  BPF  due  to  the  shock  waves  at  the  tip  of  the 
fan  blades.  These  tones  are  commonly  referred  as“  multiple  tones”,  or  "combination 
tones”.  However,  these  tones  will  not  be  present  in  this  research  because  the  turbofan 
simulator  is  always  operated  at  subsonic  conditions.  Broadband  noise  is  generated  by 
the  interaction  of  random  disturbance,  e.g.,  those  associated  with  turbulence  boundary- 
layer  flows  over  rotor  blades,  stators,  shroud,  ingested  atmospheric,  and  blade  tip  vortex 
shedding.  Broadband  noise  tends  to  dominate  the  spectrum  at  low  flow  coefficients 
when  the  stage  is  heavily  loaded  and  is  very  near  to  or  in  stall  [46].  The  BPF  and 
harmonic  tones  are  generated  by  the  aerodynamic  interaction  of  wakes  from  either 
blades  or  vanes  with  downstream  vanes  or  blades,  respectively.  Studies  have  shown 
that  one  of  the  most  significant  noise  source  interactions  is  due  to  the  upstream  IGY 
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wakes  and  rotor  blades,  i.c.,  IGY-rotor  interaction.  Other  tonal  noise  sources  may 
include  the  rotor-EGY  interaction  and  rotor-rotor  interaction  as  well  [47.  48]. 

The  noise  generation  mechanism  associated  with  IGY-rotor  interaction  is  the  unsteady 
lift  experienced  by  the  blades  when  they  pass  through  the  upstream  IGY  wakes.  There¬ 
fore,  blade  vibrations  and  tonal  fan  noise  shares  as  a  common  forcing  function  the 
upstream  IGY  wakes.  In  other  words,  the  IGY  wakes  cause  both  blade  vibrations  and 
tonal  fan  noise  as  shown  in  Figure  3.1.  Therefore,  when  the  wake  is  fully  re-energized  by 
TEB,  both  blade  vibration  and  fan  noise  will  be  minimized.  This  implies  that  acoustic 
sensing  approaches  can  be  employed  in  active  TEB  control  systems  to  re-cncrgizc  the 
IGY  wakes  for  the  reduction  of  blade  vibrations  and  HCF.  From  this  point  of  view, 
microphones  used  for  acoustic  sensing  can  be  viewed  as  the  indirect  "wake  sensor” 
and  thus  offers  the  potential  for  the  implementation  of  active  TEB.  As  shown  in  the 
latter  chapters,  the  acoustic  sensing  strategy  will  be  successfully  implemented  on  the 
small  scale  engine  simulator  which  becomes  a  significant  contribution  of  this  research 
endeavor. 

It  is  also  worth  to  point  out  that  TEB  technique  brings  us  an  additional  benefit  which 
is  the  attenuation  of  engine  fan  noise.  In  effect,  trailing  edge  blowing  from  stator  vanes 
or  rotor  blades  has  been  successfully  demonstrated  as  a  new  method  to  reduce  engine 
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Figure  3.1:  The  relationship  between  wake,  tone  noise  and  blade  vibration (HCF). 


fan  noise  in  a  realistic  rotating  environment  [19,  24,  25]. 


3.2  Turbofan  Noise  Theory 


Turbofan  tonal  noise  theory  includes  the  generation  and  propagation  of  acoustic  modes. 
The  first  aspect  is  the  understanding  of  how  the  acoustic  modes  are  generated  while 
the  second  is  how  the  acoustic  modes  propagate  along  the  inlet  duct. 


As  mentioned  above,  the  acoustic  spectra  of  the  fan  noise  is  dominated  by  tones  at 
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BPF  and  its  harmonics.  The  BPF  is  computed  as: 

/bpf  =  ^  0  [Hz]  (3.1) 

where  Q  is  the  speed  of  the  rotor  in  revolution-per-minute  and  Nr  is  the  number  of 
rotor  blades. 


3.2.1  Duct  Acoustic  Modes 


The  homogeneous  wave  equation  for  sound  pressure  in  a  moving  media  is  given  as  [49]: 

<s2» 

where  V2  is  the  Laplacian  operator;  p(r,9,z,t)  is  the  sound  pressure  in  cylindrical  co¬ 
ordinate-  z  and  9  are  the  radial,  axial,  and  circumferential  coordinates,  respectively, 
as  shown  in  Figure  3.2;  c  is  the  sound  speed  in  the  media;  and  V  =  cMk  is  the  flow 
velocity  field  vector  in  the  ^-direction  where  M  is  the  flow  Mach  number. 

By  solving  the  homogeneous  Eqn.  (3.2),  the  general  solution  for  the  sound  field  in  the 
duct  can  be  expressed  as  a  linear  combination  of  eigenfunctions  [50]: 

Mr  Mr 

p(z.  r,  6,  i)  =  £  J2  ‘W 

rn=0  /i=0 


(3.3) 


Jimvei  Feng  Chapter  3 .  Turbofan  Noise  Theory  50 


Figure  3.2:  Sound  field  in  cylindrical  coordinates. 

where  p(z,r,6,t)  denotes  the  tonal  acoustic  pressure  generated  at  the  frequency  of 
n  *  / bpf »  where  n  =  1,2,3- ••  indicates  the  harmonics  of  the  BPF:  and  9)  is 

the  acoustic  mode  with  a  complex  amplitude  Amfl.  The  acoustic  mode  is  defined  by 
two  indices  m  and  p  representing  circumferential  and  radial  variation,  respectively: 
Mc  is  the  number  of  circumferential  modes  while  Mr  is  the  number  of  radial  modes 
included.  The  axial  propagation  characteristic  of  the  (m,p)  mode  is  defined  by  the 
axial  wavenumber  which  will  be  described  in  Section  3.2.2. 

The  acoustic  mode  shape  6)  can  be  found  by  solving  the  wave  equation  (3.2) 

together  with  the  boundary  conditions.  For  example,  in  a  conccntrical  duct  with  rigid 
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walls,  the  following  boundary  conditions  should  be  satisfied: 


dp 

dr¬ 

op 

~dr 


=  0 


=  0 


(3.4) 

(3.5) 


J  r=b 


where  a  is  the  outer  radius  and  b  is  the  inner  radius  of  the  annular  duct  as  shown  in 
Figure  3.2.  Equation  (3.4)  and  (3.5)  basically  state  that  the  radial  particle  velocity 


vanishes  at  r  =  u  and  r  =  b  since  rigid  wall  is  assumed.  From  the  boundary  conditions, 
the  characteristic  equation  for  the  eigenvalue  problem  can  be  found[49.  50]  as: 


^  m  (kynijd)  /o  p\ 

Jm(kmflb)  ~  Ym(km,b) 

where  Jm(-)  is  the  first  kind  m-order  Bessel  function  and  !'„,(•)  is  the  second  kind 
m-ordcr  Bessel  function;  J'm  and  Ym  mean  the  derivative  of  Jm  and  Y,„,  respectively: 
and  km)1  is  the  root  of  Eqn.  (3.6).  It  is  noted  that  for  each  rn  there  is  an  infinite  set 
of  p,  i.e.,  p  =  0. 1, 2  •  •  • .  Since  m  is  associated  with  the  circumferential  coordinate  9 
as  will  be  seen  in  Eqn.  (3.8),  m  is  usually  referred  as  the  circumferential  order  of  the 
acoustic  mode  while  p  is  called  the  radial  order  for  the  similar  reason.  Hence,  one 
circumferential  mode  is  associated  with  a  series  of  radial  modes.  To  this  end.  the  mode 
shape  9)  can  be  expressed  as  [50]: 


$WM)  =  ©(0) 


Jm{f'  hiiiT) 


JjJk’nnO) 


^  rn(km/ia) 


Ym(k 


mu 1 


b  <  r  <  a 


(3.7) 


Because  the  sound  field  in  the  circumferential  direction  must  be  periodic,  the  function 
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0(9)  =  cos(m9  + 


(3.8) 


where  the  initial  mode  phase  <z>m,,  is  a  constant.  Substitution  of  Eqn.  (3.8)  and  Eqn. 
(3.7)  into  Eqn.  (3.3)  leads  to  the  complete  sound  field  expression: 


Mc  mt 

p(z,r,9.t)  =  EE  Am)icos(m0  +  <j>mtx)- 

m— 0  /i~0 


Ym(km,a) 


5  mi^mfi^') 


(3.9) 


It  is  convenient  to  illustrate  the  modal  behavior  by  considering  a  simple  case  of  the 
circular  duct  with  6  =  0.  The  mode  shapes  and  the  characteristic  equation  can  be 
found  from  Eqn.  (3.7)  and  Eqn.  (3.6)  by  setting  6  to  zero.  They  are 

j'jkmila)  =  0  (3  10). 

and 

^mn(^i  9)  ~~  cos(m0  T  ^m/i)  (3.11) 

Thus,  equation  3.10  allows  to  solve  for  the  roots,  kmf, .  As  mentioned  before,  for  each 
value  of  the  circumferential  order  m  there  exists  an  infinite  sequence  of  radial  modes. 
Table  3.1  contains  a  few  of  the  values  for  k„n,a.  It  should  be  noted  that  the  first  radial 
mode  is  p,  —  0  because  the  plane  wave  is  usually  referred  as  mode  (0, 0). 
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Tabic  3.1:  The  values  for  kmfla  for  a  circular  duct  with  a  radius  of  a. 


m\  p 

0 

1 

2 

3 

0 

0 

3.83 

7.02 

10.17 

1 

1.84 

5.33 

8.54 

11.71 

2 

3.05 

6.71 

9.97 

13.17 

3 

4.20 

8.02 

11.35 

14.59 

Figure  3.3  illustrates  several  acoustic  modes  where  the  initial  phase  </>miJl  is  assumed  to 
be  zero.  The  nodal  lines,  as  denoted  by  dotted  lines  in  the  figure,  occur  at  the  zeros 
of  the  mode  shape  $mfl  (r,  6)  in  circumferential  and.  radial  directions.  It  is  apparent 
that  the  integer  m  equals  the  number  of  radial  nodal  lines  while  the  integer  p  equals 
the  number  of  azimuthal  circles.  The  symbol  ;‘+”  in  the  figure  means  that  the  sound 
pressure  in  the  specified  region  is  out  of  phase  as  compared  to  that  in  the  region  denoted 
by 


Tyler  and  Sofrin  [47]  used  the  above  model  to  analyze  the  BPF  and  its  harmonics 
generated  in  a  compressor.  The  model  considered  a  stator  and  a  rotor  with  a  rotating 
speed  of  Q  in  the  positive  0-direction  in  a  circular  duct.  They  came  up  with  two 
important  results  in  their  work.  First,  they  realized  that  the  acoustic  modes  from  the 
turbofan  engine  are  not  stationary  as  discussed  above,  instead  they  are  the  spinning 


Jimvei  Feng 


Chapter  3.  Turbofan  Xoise  Theory 


54 


k  !•  ' _ y±  ■■  ^ 

k-  r  ksr 

J.UKJraj  Mr;atcf»P(Oi  (J  0.ss  il-cuf  j-Vl 

Mack  f  0, 1 )  M  vd*  ■  l  .U  ■  M  nrlc  <X  a « 


Figure  3.3:  The  mode  shapes  for  the  modes  (0.1),  (1,0).  (2,0). 


modes  at  the  speeds  related  to  the  rotor  speed,  i.e,  Eqn.  (3.8)  should  be  changed  to: 

0(0)  =  cos(m9  —  2nnNB^lt/60  +  (3-12) 


where  is  the  constant  initial  phase  of  mode  (m,  //).  It  is  convenient  here  to  calculate 
the  mode  angular  velocity  if  the  above  equation  is  re-written  as: 


0(0)  =  cos 


,n  2ir uNbQ  , 
m{-°  ~  60m  t]  +  0m,‘ 


(3.13) 


The  mode  angular  velocity  is  determined  to  be: 


Q 


rn/x 


‘2'ktiNbQ 

60  m 


(3.14) 


It  indicates  that  the  mode  angular  velocity  depends  on  the  circumferential  order  m. 


harmonic  index  n  and  rotating  speed  of  rotor  Q,  and  it  is  independent  of  the  radial 
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mode  order  p.  Hence,  the  acoustic  modes  with  a  same  m  but  different  harmonic  order 
n  spin  at  different  velocities.  It  is  also  noted  that  the  spinning  modes  usually  rotate  at 
a  different  velocity  compared  to  the  blade  passing  speed  or  its  harmonic  speed  which 
is  2;rnArBft/60  at  the  frequency  of  nfBpF-  In  addition,  positive  value  of  m  denotes  the 
mode  spins  in  the  rotor  rotating  direction,  i.e.  positive  0-direction  while  a  negative 
value  of  m  means  the  mode  spins  in  the  rotor  counter-rotating  direction,  i.e.  negative 
0-direction. 


Based  on  Eqn.  (3.12)  and  Eqn.  (3.9),  the  complete  description  of  the  acoustic  pressure 
at  u/epf  in  the  hard- walled  inlet  of  turbofan  engines  becomes: 


Me  Mr 

p(z,r,0,  t )  =  £X\4m„cos(m0  -  2irnNB^t/60  +  <z>m/,)- 

m= 0  n=0  * 


T  (l-  r\  _  y  ,,  , 

)  \rt(L  \Im\rvmii,J 


£  —  .mp 


z  ^i2nnf  QpFt 


(3.15) 


The  second  discovery  they  made  is  that  the  circumferential  mode  number  rn  in  Eqn. 
(3.15)  is  governed  by: 

m  =  n\B  +  qS  (3.1G) 

where,  q  is  an  integer  index  ±1,  ±2,  ±3.  •  •  •  and  S  is  the  number  of  upstream  or  down¬ 
stream  stator  vanes.  Thus,  the  number  of  rotor  blades  along  with  the  number  of  stator 
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vanes  determines  the  circumferential  order  of  the  modes  present  in  the  inlet. 


3.2.2  Propagation 


After  the  acoustic  mode  is  generated,  we  need  to  address  whether  or  not  the  mode 
can  propagate  down  the  inlet  duct.  As  evident  in  Eqn.  (3.15),  the  propagation  of 
the  modes  depends  on  the  axial  wavenumber  kZt,nfl.  If  k:jnfl  is  a  real  number,  the 
term  e~'k:-,n/,z  in  Eqn.  (3.15)  is  periodic  with  respective  to  2  so  that  it  represents  a 
propagating  mode,  i.e.,  the  mode  is  said  to  be  cut-on.  On  the  other  hand,  if  kz,mu  is 
imaginary  and  negative,  the  term  e~'kz-m>,z  results  in  an  exponentially  decaying  mode, 
i.e.,  the  mode  is  said  to  be  cut-off.  It  is  noted  that  the  above  criteria  holds  for  both 
the  stationary  and  spinning  modes  since  0(0)  (see  Fqn.  (5.8)  for  stationary  modes, 
and  Eqn.  (3.12)  for  spinning  modes  )  does  not  play  a  role  in  the  exponential  term  in 
Eqn.  (3.15). 


To  find  the  axial  wavenumber  kz,mn •  Eqn.  (3.2)  is  written  in  cylindrical  coordinate  and 
assuming  harmonic  motion  results  in: 


d2p  l  dp  1  d2p  d2p  2  dp  .,->d2P 

3 — ~ — t — o  +  •x-rr  =  ~k0  +  2ikf)M  —  +  M  7^ 


dr2  r  dr  r2  dO2  dz2 


(3.17) 
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where  k0  is  the  free  space  wavenumber  as  given  by  k0  =  27r/0/c  and  /o  is  the  excitation 
frequency  which  is  either  BPF  or  its  harmonics  in  this  case. 

Substituting  Eqn.  (3.9)  or  (3.15)  into  Eqn.  (3.17)  gives: 

C  =  H  -  *?,„„(  1  -  M2)  -  2 k0ks.mrM  (3.18) 

Eqn.  (3.18)  allows  to  solve  for  the  value  of  the  wavenumber  The  two  roots  are 

found  to  be: 

K„  =  [-3/io  +  -  «„,(!  -  m  ]  /(l  -  -V2)  (3-19) 

Km  =  [-W*o  -  -  *4,(1  -  m  }  /(l  -  -V2)  (3.20) 

where  and  k~mft  represent  positive  and  negative  traveling  or  decaying  waves, 

respectively. 

From  equations  (3.19)  and  (3.20),  there  are  three  conditions  that  define  the  propagation 
characteristic.  Focusing  only  on  the  modes  in  the  positive  ^-direction  (i.e..  Eqn.  (3.19)), 
they  are: 

(a)  If  ko  >  kmfIy/ 1  -  A/2,  kZm)J  is  a  real  number  and  the  mode  will  propagate  with  a 
wavenumber  of 

Km,  =  [-Uko  +  yfki  -  KJl  -  A P)  ]  /(I  -  A/2) 


(3.21) 
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( b )  If  k0  <  k-mfiV  1  -  -I/2.  is  complex  with  negative  imaginary  component  and  the 

mode  will  exponentially  decay  with  a  wavenumber  of 


k+  = 

nz,mfi 


-Mk0  -  isjk^il  -  IF)  -  hi  ]  /(I  -  M2)  (3.22 


(c)  k(j  =  kmiJ  \/l  —  A/2  defines  the  transition  between  the  previous  two  cases.  This 
condition  is  used  to  define  the  mode  cut-on  frequency  /™,_on  as: 

fcut-°n  =  -  W  (3.23) 

Therefore,  an  acoustic  mode  can  propagate  if  its  excited  frequency  is  greater  than 
fmfi  ~  °n  while  it  decays  otherwise. 


In  addition,  it  is  also  noted  that  the  mode  axial  propagation  velocity  or  mode  phase 
velocity  cz,m and  wavelength  A2,m/1  are  greater  than  those,  i.e.  c  and  A  =  c//0  in  the 
free  field.  They  are 


ym.fi 


X, 


V1  z  (W/o)2 
A 


(3.24) 

(3.25) 


^-(W/o)2 

respectively.  Eqn.  (3.24)  and  Eqn.  (3.25)  indicate  that  the  mode  axial  phase  velocity 


and  mode  axial  wavelength  vary  from  mode  to  mode. 
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3.2.3  Acoustic  Mode  Analysis  of  Simulator 


Based  on  the  physical  configuration  of  the  simulator,  the  excited  acoustic  inodes  in  the 
inlet  can  be  found.  In  addition,  the  cut-on  mode  frequencies  can  be  calculated;  and 
hence,  the  propagating  characteristic  of  the  modes  in  the  inlet  can  be  determined.  Two 
simulator  speeds  of  30, 000  and  40, 000  rpm  are  considered  (the  simulator  is  operated 
from  20,000  rpm  and  40,000  rpm  in  this  research).  The  corresponding  BPF  and  2BPF 
at  30,000  rpm  are  9  kHz  and  18  kHz,  respectively  while  they  are  12  kHz  and  24  kHz, 
respectively  at  40,000  rpm  according  to  Eqn.  (3.1). 

The  simulator  is  equipped  with  4  upstream  IGVs,  18  rotor  blades  and  26  downstream 
EGVs.  Both  IGV-rotor  and  rotor-EGY  interactions  result  in  the  generation  of  discrete 
tones.  Using  Eqn.  (3.16),  the  circumferential  modes  excited  by  IGY-rotor  interaction 
are  listed  in  Table  3.2.  The  modes  excited  by  by  rotor-EGY  interaction  are  presented 
in  Table  3.3.  Note  in  Eqn.  (3.16),  the  integer  index  q  can  range  from  — oc  to  oc.  Only 
a  limited  range  of  q  is  used  in  both  of  the  tables  and  thus  only  a  few  circumferential 
modes  are  computed.  As  indicated  by  Table  3.2.  the  mode  order  generated  by  IGV- 
rotor  interaction  increases  by  an  interval  of  4.  On  the  other  hand,  the  mode  order 
generated  by  rotor-EGY  interaction  increases  by  an  interval  of  26.  Note  that  these  two 
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Table  3.2:  The  circumferential  order  of  some  acoustic  modes  excited  by  1G\  -rotor 
interaction  in  the  simulator  inlet. 


BPF  (n=l) 

2BPF  (n=2) 

3BPF  (n=3) 

Q 

m  =  18  +  4q 

Q 

-10 

BO 

-12 

-6 

BO 

-8 

-6 

B 

-2 

-4 

m 

-2 

B 

2 

D 

Q| 

2 

-3 

6 

-8 

4 

-12 

6 

-2 

10 

-7 

8 

-11 

10 

-1 

14 

-6 

12 

14 

Table  3.3:  The  circumferential  order  of  some  acoustic  modes  excited  by  rotor-EGY 
interaction  in  the  simulator  inlet. 


BPF  (n=l) 

2BPF  (n=2) 

3fcSPF  (n=3) 

Q 

m  =  18  4-  2 6q 

Q 

m  =  36  +  26g 

Q 

m  =  54  4-  26q 

-2 

-34 

-2 

-16 

-3 

-24 

-i 

-8 

-1 

10 

-2 

2 

0 

18 

0 

36 

-i 

28 

1 

44 

1 

62 

0 

54 

2 

70 

2 

88 

1 

80 
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tables  only  provide  information  on  the  circumferential  order  of  the  acoustic  modes.  As 
stated  before,  a  specific  m  consists  of  a  sequence  of  radial  order,  denoted  by  p. 


The  calculation  of  the  mode  cut-on  frequency.  f™tl~on  can  be  carried  out  according 
to  Eqn.  (3.23)  where  kmil  can  be  found  by  solving  Eqn.  (3.6).  The  outer  radius  of 
a  =  2.05  inches,  the  hub-to-tip  ratio  of  b/a  =  0.439,  the  sound  speed  at  25°C  of 
c  =  346.6  meter/second,  and  the  estimated  Mach  number  of  0.3  at  30.000  rpm  and 
0.35  at  40,000  rpm  are  used  in  the  calculation.  The  results  are  listed  in  the  tables  3.4 
and  3.5  for  the  two  simulator  speeds.  As  can  be  seen,  the  cut-on  frequency  increases 
with  the  circumferential  order  m  of  the  acoustic  modes  with  same  radial  order  //. 
Similarly  for  the  acoustic  modes  with  same  circumferential  order,  the  cut-on  frequency 
increases  with  p.  Therefore,  the  relative  smaller  cut-on  frequencies  are  positioned  in 
the  upperleft  corner  of  the  tables. 

Based  on  Table  3.4  and  Table  3.5  and  the  wave  propagation  criteria  discussed  in  section 
3.2.2  the  propagating  acoustic  modes  in  the  simulator  inlet  duct  can  be  determined. 
Table  3.6  lists  the  propagating  acoustic  modes  at  BPF  and  2BPF  from  both  IGY- 
rotor  and  rotor-EGV  interactions  when  the  simulator  runs  at  30, 000  rpm  and  40, 000 
rpm.  Note  that  in  the  table.  (2,0  -  2)  means  the  range  of  modes  (2.0)  through  (2,2). 
For  example,  the  modes  (2,0)  and  (2,1)  are  present  in  the  inlet  at  30,000  rpin,  i.e., 
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Table  3.4:  The  mode  cut-on  frequencies  at  the  speed  of  30,000  rprri  (M=0.3)  [Hz]. 


fi\  in 

0 

±2 

±4 

±6 

±8 

±10 

±12 

0 

0 

2,811 

5,308 

7,  570 

9.748 

11,896 

14.026 

1 

5,796 

6,673 

8,874 

11,552 

14,157 

16.594 

18,938 

2 

11,394 

11,818 

13.046 

14,984 

17,461 

20,147 

22.756 

3 

17,029 

17,307 

18.126 

19,450 

21,244 

23,480 

26,054 

4 

22,676 

22,883 

23,498 

24,501 

25,868 

27,582 

29.639 

Table  3.5:  The  mode  cut-on  frequencies  at  the  speed  of  40.000  rpm  (M=0.35)  [Hz]. 


fi\  m 

0 

±2 

±4 

±6 

±8 

±10 

±12 

0 

0 

2,761 

5,212 

7,433 

9,573 

11.681 

13,773 

1 

5,692 

6,553 

8,715 

11,344 

13,902 

16,295 

18,596 

2 

11,189 

11.605 

12;811 

14,714 

17,146 

19,784 

22.346 

3 

16,722 

16,995 

17.799 

19.099 

20,861 

23,057 

25.585 

4 

22,267 

22,471 

23,075 

24,059 

25,402 

27,085 

29,105 

Jimvei  Feng 


Chapter  3.  Turhofan  Noise  Theory 


63 


Table  3.6:  The  propagating  modes  at  BPF  and  2BPF. 


Speed 

IGV  -  rotor  Interaction 

rotor  -  EG! 

'  Interaction 

BPF  (9  kHz) 

2BPF  (18  kHz) 

BPF  (9  kHz) 

2BPF  (18  kHz) 

30,000  rpm 

i — l 

1 

O 

oi 

HH, 

(0,0-3) 

none 

(10,0-  1) 

(±6,0) 

(±4, 0  -  2) 

(M=0.3) 

(±8,0-2) 

(±12,0) 

BPF  (12  kHz) 

2BPF  (24  kHz) 

BPF  (12  kHz) 

2BPF  (24  kHz) 

40,000  rpm 

(±2,0-2) 

(0,0-4) 

(-8.  0) 

(10,0-3) 

1— ( 

1 

o 

CO 

;H 

(±4,0-4) 

(M=0.35) 

(±10,0) 

(±8,0-3) 

(±12,0-2) 

they  are  cut-on.  because  the  circumferential  order  rn  =  2  is  excited  by  the  IGY-rotor 
interaction  according  to  Table  3.2.  In  addition,  the  BPF  of  9,000  Hz  at  30,000  rpm  is 
greater  than  the  cut-on  frequencies  of  modes  (2.0)  and  (2,1),  (2,811  Hz  and  6,673  Hz, 
respectively)  as  seen  from  Table  3.4.  However,  mode  (2.2)  cannot  propagate  since  its 
cut-on  frequency  (i.e.,  11,818  Hz)  exceeds  the  BPF  (i.e.,  9,000  Hz). 

It  should  be  noted  that  the  magnitude  of  the  acoustic  modes  listed  in  the  table  depends 
on  the  distribution  of  loading  on  the  blades.  It  implies  that  some  of  the  modes  will 
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dominate  the  inlet  sound  field  while  the  other  modes  may  contribute  insignificantly. 

As  found  from  the  table,  the  number  of  the  propagating  modes  at  BPF  is  limited 
although  IGY-rotor  and  rotor-EGV  interaction  excite  an  infinite  set  of  the  circumfer¬ 
ential  order  m.  This  is  because  most  of  the  excited  modes,  especially  the  modes  due 
to  rotor-EGV  interaction,  are  cut-off.  In  addition,  more  modes  are  cut-on  at  40.000 
rpm  as  compared  to  those  at  30,000  rpm  since  the  corresponding  BPF  at  40.000  rpm  is 
higher  than  that  at  30,000  rpm.  For  the  same  reason,  the  sound  field  at  2BPF  contains 
much  more  modes  than  that  at  BPF. 

Figure  3.4  shows  several  propagating  mode  shapes.  The  solid  circle  denotes  the  inlet 
centerbody  which  has  a  radius  of  0.9  inches.  As  compared  to  Figure  3.3,  the  acoustic 
modes  here  are  spinning.  A  positi'v  value  for  m  defines  the  mode  spins  in  the  same 
direction  as  the  rotor  while  a  negative  value  is  for  the  opposite  direction.  Since  mode 
(0,1)  has  no  variation  in  ^-direction,  there  is  no  spinning  direction  shown  in  the  figure. 
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Figure  3.4:  The  mode  shapes  of  the  spinning  modes  (0,1),  (-2,0)  and  (2,0). 


3.2.4  Radiation 


When  the  propagating  acoustic  modes  reach  the  opening  of  the  inlet,  they  will  radiate 
into  the  far  field.  It  is  interesting  to  investigate  the  directional  distribution  of  the 
propagating  modes.  However,  the  radiation  process  is  complicated  due  to  the  diffrac¬ 
tion  effects  and  the  acoustic  impedance  change  at  the  inlet  opening.  The  first  work 
on  the  sound  radiation  from  an  unflanged  cylindrical  duct  may  be  from  Levine  and 
Schwinger[51]  and  Weinstein[52]  using  Wiener-Hopf  technique.  Later  work  included 
Homicz  and  Lordi[53]  and  Lansing[54].  The  earlier  work  of  Levine  and  Schwinger  fo- 
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cuscd  on  non-spinning  modes.  Later  work  by  Homicz  and  Lordi  included  the  directivity 
plot  for  spinning  modes.  According  to  their  work,  each  mode  will  produce  a  radiation 
pattern  in  the  far  field.  The  sound  pressure  level  directivity  is  then  a  composite  of  all 
the  acoustic  mode.  The  location  of  the  principal  lobe  for  the  mode  (m,  p)  in  a  circular 
duct  with  a  radius  of  a  can  be  predicted  as  [53]: 


Tn  =  co.s  1  \/l  —  M2 


l-( 


fcpfl 


)2 


(3.26) 


where  Tp  is  the  angle  of  the  principal  lobe  of  the  radiation  directivity  of  the  acoustic- 
mode  measured  from  the  inlet  axis.  Equation  3.26  is  applicable  for  aircraft  engine 
inlets  with  a  hub  to  tip  ratio  of  less  than  0.5. 


3.3  Summary 


The  basics  of  fan  noise  theory  were  reviewed  in  this  Chapter.  Some  important  concepts 
such  as  spinning  acoustic  mode,  mode  spinning  speed,  axial  mode  phase  velocity  and 
wavelength  were  introduced.  Propagation  criteria  in  the  duct  with  flow  was  described. 
The  cut-on  frequencies  for  the  small  scale  simulator  using  a  concentrical  duct  model 
were  calculated.  The  propagation  modes  at  BPF  and  2BPF  for  the  simulator  running 
at  the  speeds  of  30, 000  and  40, 000  rpm  were  tabulated.  These  knowledge  is  neces- 
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sary  when  we  design  and  implement  TEB  control  systems  using  the  acoustic 
approaches  described  in  the  Chapter  5. 


Chapter  4 


A  Pitot-probe  Based  Flow  Control 
System 


This  chapter  describes  a  TEB  control  system  with  Pitot  probes  as  wake  sensors  on  the 
small  scale  turbofan  engine  simulator.  As  mentioned  early  in  Chapter  1,  one  of  the 
research  objectives  in  the  dissertation  is  to  investigate  the  feasibility  of  implementing 
active  flow  control  systems  in  a  realistic  environment.  Though  the  Pitot  probe  in 
this  chapter  is  an  intrusive  flow  sensor  which  are  not  probably  practical  in  realistic 
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situations  of  turbofan  engines,  it  is  appropriate  as  a  preliminary  step  toward  a  more 
in  depth  investigation  of  advanced  wake  sensing  approaches.  In  addition,  an  upper 
performance  limit  in  terms  of  wake-filling  can  be  obtained  since  the  Pitot  probes  provide 
the  direct  measurement  of  the  wake  which  is  sought  to  be  minimized.  Consequently, 
the  probe-based  control  system  developed  in  this  chapter  can  serve  as  the  baseline  case 
in  evaluating  the  performance  of  other  control  systems  with  indirect  wake  sensors. 

In  this  chapter,  the  working  principle  of  Pitot  static  probes  is  first  described.  The 
experimental  setup  in  the  anechoic  chamber  follows  and  the  configuration  of  the  con¬ 
trol  system  is  then  presented.  The  wake-filling  performance  in  terms  of  aerodynamic 
and  acoustic  measurements  are  investigated.  Finally,  the  main  conclusions  are  then 
discussed. 


4.1  Experimental  setup  in  the  simulator 


The  sensing  approach  used  in  this  part  is  simple  Pitot-Static  probes  to  measure  the 
total  pressure  in  the  wakes  behind  each  IGV  as  illustrated  in  Figure  4.1.  The  probe 
has  a  head  diameter  of  1/16”  with  a  sensing  hole  of  diameter  1/32”  and  length  of  0.5”. 
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These  probes  are  insensitive  to  incidence  angle  variations  of  5°.  Four  threaded  holes 
are  provided  in  the  inlet  for  mounting  the  four  flow  sensors.  The  holes  are  positioned 
such  that  each  probe  measures  the  maximum  total  pressure  defect  in  a  wake.  The  inlet 
flow  field  is  found  to  be  radially  uniform,  and  hence  the  use  of  a  single  flow  sensor  for 
each  wake  is  adequate.  This  will  be  further  described  in  the  next  section.  These  wake 
probes  are  placed  0.5  Cs  downstream  of  the  IGV  where  Cs  represents  the  IGY  chord. 
This  spacing  is  determined  based  on  a  bench  test  result  which  shows  that  a  distance  of 
0.5  Cs  is  required  to  allow  the  jet  to  mix  into  the  wake.  Measurements  made  closer  to 
the  IGV  will  sense  the  potential  core  of  the  jet  issuing  out  of  the  blowing  hole  rather 
than  the  mixed  out  flow,  and  hence  could  generate  an  erroneous  wake  signal. 


VVitki  Vt  ;  jIk  o  ’  sj 


Figure  4.1:  The  Pitot-probe-based  experimental  setup. 
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Another  probe  is  placed  in  the  free  stream  to  measure  the  free  stream  total  pres¬ 
sure  (not  shown  in  Figure  4.1).  The  difference  between  the  output  of  the  wake  and 
free  stream  probes  forms  the  error  signal  for  a  traditional  feedback  PID  controller 
(proportional-integral-derivative),  which  is  to  be  introduced  in  a  later  section.  Basi¬ 
cally,  this  algorithm  aims  to  minimizing  the  error  signal.  Therefore,  the  wake  is  going 
to  be  filled  up  to  the  free  stream  when  the  controller  has  converged.  The  blowing  rate 
is  adjusted  by  the  MEMS-based  microvalves  according  to  the  control  signal  from  the 
controller. 


4.2  Working  principle  of  the  Pitot-static  probe 


In  order  to  employ  Pitot-static  probes  as  flow  sensors,  grasp  of  some  basic  fluid  knowl¬ 
edge  is  needed.  Therefore,  it  is  appropriate  to  introduce  a  simple  review  of  fluid 
dynamics.  The  famous  Bernoulli’s  equation  is  to  be  first  introduced.  The  working 
principle  of  the  velocity  measurement  device.  Pitot-static  probe  is  then  described. 


Bernoulli’s  equation  can  be  derived  from  the  first  law  of  thermodynamics  (conservation 
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of  energy)  to  give: 

l-pV2  +  P  =  P,  (4.1) 

where  P0  is  usually  called  total  pressure  which  remains  constant  streamwise.  p  is  the 
fluid  density,  P  is  pressure  and  V  denotes  flow  velocity.  The  term  ^pV'2  is  referred  as 
dynamic  pressure  while  P  is  called  static  pressure. 


It  is  important  to  review  the  assumptions  in  Bernoulli's  equation.  First,  the  steady  flow 
assumption  means  that  the  relationship  cannot  be  applied  to  a  flow  which  has  time 
dependent  properties  or  behavior.  There  will  be  some  exceptions  to  this  when  flow 
changes  are  slow  enough  to  consider  the  flow  “quasi-steady".  Secondly,  incompressible 
flow  was  assumed  meaning  that  density  was  assumed  constant  throughout  the  flow. 
For  most  problems  encountered  in  aerodynamics,  the  incompressibility  assumption 
can  be  used  for  flow  velocity  of  up  to  about  0.65  Mach  number.  Thirdly,  flow  along  a 
streamline  was  assumed.  Thus,  P0  may  be  considered  constant  throughout  the  entire 
flow  if  it  is  known  that  all  streamlines  emanated  from  a  common  source  such  as  a 
uniform  flow.  If  the  flow  is  known  to  have  a  velocity  or  pressure  gradient,  Bernoulli’s 
equation  must  be  used  with  care  and  only  applied  along  streamlines.  Finally,  negligible 
buoyancy  was  assumed  in  the  equation.  The  buoyancy  term  pgz  is  always  negligible 
when  Eqn.  (4.1)  is  used  in  air  where  c  denotes  the  altitude  from  sea  level  and  g  is  the 
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acceleration  of  gravity. 


In  the  simulator  setup  for  the  speed  range  of  20.000  to  40.000  rpm.  all  the  above 
assumptions  arc  satisfied.  The  axial  Mach  number  of  steady  flow  in  the  inlet  duct 
for  this  speed  range  is  only  around  0.3.  The  changes  in  simulator  speed  and  thus 
flow  conditions  can  also  be  considered  “quasi-steady” .  Finally,  the  total  pressure  P0  is 
constant  and  equal  to  the  atmospheric  pressure,  i.e.,  10,1000  Pa  or  14.3  Psi.  because 
the  simulator  is  drawing  ambient  stationary  air  into  the  inlet.  An  schematic  of  stream 
lines  for  the  simulator  case  is  shown  in  Figure  4.2. 


''  /•  -  P 
V  ~  0 


Figure  4.2:  Example  of  streamline  in  the  simulator  when  Bernoulli's  equat  ion  is  applied. 
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Bernoulli's  equation  provides  an  easy  method  for  determining  the  speed  of  a  fluid  flow. 
Rearranging  the  equation,  the  flow  velocity  is  given  as: 


V 


(2(P0  -  P) 


(4.2) 


The  total  pressure  P0  in  a  flow  is  easy  to  measure  if  the  flow  direction  is  known.  One 
needs  only  place  a  simple,  open  ended  tube  in  the  flow  such  that  the  tube  is  aligned 
with  the  local  flow  direction  and  the  open  end  of  the  tube  is  pointing  upstream.  The 
flow  then  stagnates  or  comes  to  rest  at  the  tube’s  open  end.  giving  a  local  pressure  equal 
to  the  total  pressure.  The  tube  can  then  be  connected  to  any  pressure  sensors  to  allow 
a  determination  of  that  pressure  as  shown  in  Figure  4.3(a).  This  type  of  tube  is  called 
a  Pitot  probe.  The  static  pressure  P  is  equally  easy  to  measure  using  another  type 
of  tube  with  a  closed  end  and  pressure  ‘taps”  or  openings  around  its  circumference. 
This  tube  must  also  be  aligned  with  the  flow  direction.  When  connected  to  a  sensor, 
the  flow’s  static  pressure  will  oe  indicated.  It  should  be  noted  that  for  a  static  probe 
to  obtain  good  readings  the  upstream  end  should  be  properly  shaped  to  avoid  flow 
separation  and  pressure  taps  placed  far  enough  from  the  end  to  avoid  errors  that 
might  be  produced  by  the  local  acceleration  of  the  flow  over  the  front  of  the  probe. 
The  static  probe  is  shown  in  Figure  4.3(b).  A  Pitot-static  probe,  as  shown  in  Figure 
4.3(c),  is  simply  a  static  tube  wrapped  around  a  Pitot  probe. 
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Figure  4.3:  The  operation  of  (a)  Pitot  probe,  (b)  Static  probe,  and  (c)  Pitot-static 
probe. 

The  Pitot  probes  are  used  in  the  experiments  which  are  shown  in  Figure  4.4.  One 
probe  is  placed  in  the  free  stream  while  another  probe  in  the  wake.  Therefore,  the 
total  pressure  measured  in  the  free  stream  probe  is 

Poj  =  Pf  +  \pVj  (4.3) 

where  Pj  is  the  free  stream  static  pressure  and  V)  is  the  free  stream  velocity.  On  the 
other  hand,  the  total  pressure  in  the  wake  is  measured  as 

Po,„,  =  Pu,  +  \pV*  (4-4) 

where  Pw  is  the  wake  static  pressure  and  Vw  is  the  wake  flow  velocity.  As  described  in 
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Chapter  1.  the  velocity  Vw  is  less  than  Vf  due  to  the  viscous  effect  as  flow  passes  the 
surface  of  IGY.  However,  the  static  pressure  in  the  wake  is  the  same  as  that  in  the  free 
stream,  i.e..  Pw  —  P )  since  streamline  flow  field  is  assumed  in  the  inlet  [55],  Therefore. 
The  wake  total  pressure  P0,u>  is  less  than  the  free  stream  total  P0j  because  of  the  losses 
in  the  flow  velocity.  The  difference  between  these  two  total  pressures  is  used  as  the 
error  signal  for  the  control  system,  which  will  be  described  in  the  following  section. 


Figure  4.4:  The  use  of  Pitot  probes  in  the  experiments. 

The  error  signal  P0,/  ~  Po,w  is  minimized  in  the  control  system  by  using  trailing  edge 
blowing  technique  adding  mass  flow  into  the  wake,  thus  augmenting  the  wake  flow 
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velocity.  Especially  as  the  digital  controller  converges,  the  error  signal  will  vanish.  Le¬ 
the  total  pressure  of  the  wake  with  TEB,  PB,f,B  will  be  equal  to  Pq.j.  That  is 

Po,uB  =  Po,f  (4.5) 


Note  that  PBBB  is  defined  as: 

Pq,wB  =  Pw  +  ^P^w.TEB  (4-6) 

It  is  apparent  from  Eqn.  (4.6)  and  Eqn.  (4.3)  that  the  augmented  wake  flow  velocity 
Vw,teb  will  be  equal  to  the  free  stream  velocity  V).  In  consequence,  complete  wake 
re-energizing  is  achieved. 


4.3  Proportional  integral  derivative  feedback  con¬ 
troller 


The  control  scheme  selected  for  this  application  is  the  feedback  PID  algorithm  [56.  57]. 
Figure  4.5  shows  a  schematic  of  a  PID  controller  in  a  typical  feedback  control  block 
diagram  which  includes  the  plant,  the  sensor,  and  the  reference  input  r(t).  The  error 
signal  e(t)  is  defined  as  the  difference  between  the  reference  signal  r(t)  and  the  sensor 
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output.  The  output  of  the  PID  controller  is  usually  called  control  signal  v(t)  which 
drives  the  plant. 


reference,  n't) 


error,  eitl  ;  p jjyj 

r' s  • - * 


controller 


Control  signal. 

inti 

Plant 


Sensor  i-*~ 


Figure  4.5:  A  typical  feedback  control  system. 


Feedback  PID  controllers  are  widely  used  in  the  process  control  industries.  The  control 
output  u(t)  of  a  PID  controller  consists  of  the  weighted  contribution  of  three  terms. 
These  three  control  signal  components  are  proportional  to  the  error  signal  e(t)  and  its 
integral  and  derivative  as: 

u(t)  =  kpe(t)  +  kt  jf  e{t)dn  +  k/-^  (4.7) 

where  t  denotes  time:  kp,  kt  and  kd  represent  the  proportional,  integral  and  derivative 
gains,  respectively.  Proportional  control  can  lead  to  reduced  errors  to  disturbance  but 
it  results  in  small  steady-state  errors.  It  can  also  increase  the  speed  of  response  but 
typically  at  the  cost  of  a  larger  transient  overshoot.  If  the  controller  also  includes  a 
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term  proportional  to  the  integral  of  the  error,  the  steady-state  error  can  be  eliminated. 
However,  there  tends  to  be  a  further  deterioration  of  the  dynamic  response  called 
windup  that  requires  a  long  period  of  time  after  the  error  changes  its  sign  to  bring  the 
system  back  to  normal.  This  is  the  result  of  the  integral  term  becoming  extremely  large, 
thus  requiring  a  large  number  of  samples  (of  opposite  sign)  to  reduce  the  magnitude 
of  the  integral  term.  Finally,  addition  of  a  term  proportional  to  the  error  derivative 
can  add  damping  to  the  dynamic  response  to  improve  stability.  A  disadvantage  of 
the  derivative  term  is  that  it  can  amplify  the  high-frequency  noise  in  the  system.  In 
practice,  values  of  the  PID  parameters  are  mainly  determined  empirically.  These  three 
terms  combined  form  the  classical  PID  controller.  Therefore,  design  of  a  PID  controller 
involves  the  determination  of  the  three  gains.  Figure  4.6a  shows  the  structure  of  the 
analog  PID  controller  while  Figure  4.6b  illustrates  it  in  the  s-domain  obtained  by  the 
Laplace  transform.  The  s  domain  expression  for  Eqn.  (4.7)  is  given  as: 

U{s)  =  (kp  +  ^  +  kds^j  E(s)  (4.8) 

where  U(s)  and  E(s)  represent  in  the  s-domain  the  signals  in  the  time  domain.  u(t) 
and  e(t),  respectively. 

In  recent  years,  digital  signal  processors  (DSP)  arc  usually  used  to  implement  PID 


controllers  since  DSP^  become  more  affordable.  This  requires  that  the  analog  propor- 
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Figure  4.6:  Block  diagram  of  analog  PID  controller  in  (a)  time  domain  and  (b)  s-domain 
froms. 


tional,  integral  and  differential  terms  in  Eqn.  (4.7)  be  approximated  numerically  in 
a  discrete  system.  If  the  numerical  approximation  is  performed  accurately,  there  is 
little  difference  in  the  system  response  from  the  analog  implementation.  The  Euler 
approximation  method,  or  the  right-side  rectangular  rule,  which  is  the  simplest  widely 
used  approach,  results  in  the  controller  difference  eqinuon  [57]: 

u(k  +  1)  =  kpe(k  +  1)  +  ktu,{k  +  1)  +  ±11  ~  (4.9) 

where  k  is  the  time  index,  i.e.,  tk;  T  is  the  sampling  period:  and  uj[k  +  1)  is  given  by: 

ui{k  +  1)  =  Uj(k )  +  Te{k  +  1)  (4.10) 

Alike  the  analog  PID  algor:' hm  which  has  the  expression  in  the  .s-domain,  the  digital 
PID  controller  can  be  expressed  in  the  discrete  2-domain  through  the  2-transform  of 
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Eqn.  (4.9).  The  2-domain  expression  is  useful  in  the  analysis  of  the  control  system 
performance  such  as  stability,  response  time.  etc.  It  is  given  as 


U(z) 


kp  + 


kjTz 
2  —  1 


.  Uz  -  1) 
Tz 


E(z) 


(4.11) 


Figure  4.7:  Digital  PID  controller. 


A  block  diagram  of  the  application  of  the  PID  controller  to  the  TEB  is  shown  in  Figure 
4.8.  The  control  loop  starts  from  the  two  Pitot  probes.  The  pressure  sensors  are  needed 
to  convert  the  probe  total  pressures  to  voltage  signals.  Differential  voltage  amplifiers 
with  a  gain  of  10  arc  implemented  at  the  output  of  the  pressure  sensors  (the  probe 
sensitivity  is  about  20  mVolts/psi).  Differential  amplification  is  used  since  the  signal 
from  a  Pitot  probe  sensor  needs  to  be  converted  to  a  single  ended  ground  referenced 
signal.  This  is  necessary,  since  all  other  devices  in  the  system  are  ground  referenced. 
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Low  pass  filters  (LPF)  with  a  cut-off  of  10  Hz  arc  used  to  serve  as  anti-aliasing  filters 
to  remove  the  high  frequency  noise.  The  digital  PID  algorithm  given  in  Eqn.  (4.9) 
is  implemented  on  a  Texas  Instrument  DSP  TMS320C30.  The  sampling  frequency  is 
chosen  to  be  30  Hz.  This  is  the  lowest  sampling  frequency  which  can  be  obtained 
using  the  on-board  timer.  However,  this  frequency  is  high  enough  as  compared  to  the 
slow  response  of  the  mechanical  components  in  the  system  such  as  the  MEMS  based 
microvalves,  blowing  air  travelling  in  the  tubing  and  the  Pitot  probes.  Since  the  output 
voltage  range  of  the  DSP  board,  i.e.,  -8.2  —  8.2  volts  does  not  match  the  operating 
range  of  the  MEMS  based  microvalves  as  shown  in  Figure  2.10.  a  voltage  offset  circuit 
with  4.5  Volts  are  used  for  each  microvalve  to  shift  the  range  of  -8.2  —  8.2  volts  to 
-3.7  —  12.7  Volts.  In  practice,  this  range  are  further  limited  so  that  it  only  covers  the 
effective  working  range  of  the  microvavles.  which  is  3  —  9  Volts.  Since  six  microvalves 
are  used  in  the  experiments  to  control  independent lv  the  six  blowing  holes  in  each  of 
the  four  IGVs,  24  channels  of  the  offset  circuits  are  required. 

The  blowing  air  controlled  by  the  microvalvcs  mixes  with  the  flow  in  the  wake  which  is 
sensed  again  by  the  flow  sensor.  Therefore,  a  single-input-single-output  (SISO)  control 
loop  is  formed.  In  the  experiments,  all  six  microvalves  connected  to  a  particular  IGV 
are  driven  by  the  same  SISO  controller.  Since  there  are  four  IGVs  in  the  simulator  inlet, 
four  identical  SISO  PID  controllers  are  required  which  hence  forms  a  4140  controller. 
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Figure  4.8:  Block  diagram  of  the  PID  controller  in  anechoic  chamber  setup. 

Before  testing  the  4140  controller  on  the  simulator  in  the  anechoic  chamber,  prelim¬ 
inary  tests  of  the  SISO  system  are  conducted  in  a  bench  setup.  The  purpose  of  the 
bench  tests  is  to  determine  the  appropriate  values  for  kp,  ki  and  kd  of  the  controller. 
These  values  are  difficult  to  predict  for  two  reasons.  The  first  reason  is  that  the  exact 
model  for  MEMS  based  microvalve  in  terms  of  blowing  air  total  pressure  versus  con¬ 
trol  voltage  is  not  available.  As  discussed  in  section  2.5.2.  the  pneumatic  microvalve 
exhibits  a  nonlinear  and  time-varying  characteristics.  The  second  difficulty  lies  in  that 


the  transfer  function  from  the  blowing  air  total  pressure  to  the  total  pressure  sensed 
by  the  flow  sensor  is  also  difficult  to  obtain.  This  is  because  jet  mixing  in  the  wake  is 
a  very  complicated  process  [58].  As  in  many  applications,  these  three  parameters  are 
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selected  experimentally  [57]. 
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Figure  4.9:  Bench  setup  for  BID  TEB  control. 


The  bench  setup  is  shown  in  Figure  4.9.  It  consists  of  two  parallel  Plexiglas  plates  with 
the  IGV  mounted  between  the  plates.  This  assembly  is  completed  by  side  supports 
that  are  fastened  onto  the  plates  by  bolts.  The  bench  is  then  mounted  on  a  wooden 
test  stand,  and  the  assembly  is  placed  at  the  exit  of  a  4”  x  4”  nozzle  that  is  connected 
to  the  discharge  of  a  centrifugal  blower.  The  bench  is  aligned  so  that  the  midspan  of 
the  IGV  lines  up  with  the  midspan  of  the  nozzle  cross  section.  One  Pitot  static  probe 
is  placed  in  the  free  stream  to  serve  as  the  reference  signal.  Another  probe  is  placed 


so  that  its  center  collinear  with  the  center  of  a  blowing  hole.  In  addition,  the  probe  is 
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located  0.5  chord  downstream  IGY  so  that  it  measures  the  full}'  mixed  flow  as  shown 
in  Figure  2.7.  Free  stream  velocity  can  be  varied  by  putting  a  blockage  at  the  inlet  of 
the  blower. 

Figure  4.10  and  Figure  4.11  present  the  performance  of  the  proportional  part  of  the 
PID  controller,  i.e.,  P  controller  with  a  gain  of  kp  =  0.6  and  2,  respectively.  The  free 
stream  flow  velocity  is  changed  abruptly  to  test  the  step  response  of  the  system.  The 
long  convergence  time  is  due  to  the  response  time  of  the  MEMS  based  microvalve  and 
the  traveling  time  of  blowing  air  in  the  tube.  The  steady  state  error  is  nonzero  which 
decreases  as  kp  increases  at  the  price  of  the  greater  control  signal  variation. 

Figure  4.12  and  Figure  4.13  demonstrate  the  system  performance  when  PI  controllers 
are  used.  With  the  introduction  of  the  integral  term,  die  steady  state  error  is  removed. 
As  ki  increases,  the  control  signal  reaches  the  steady  state  more  quickly,  however, 
usually 'at  the  cost  of  overshooting.  The  overshoot  is  not  seen  in  the  figures  because 
the  pressure  of  the  blowing  supply  air  is  not  high  enough.  Similarly,  the  effect  of  the 
derivative  term  is  also  tested  which  however  does  not  cause  obvious  difference  in  the 
system  step  response  since  the  actuation  part  has  a  very  long  response  time. 
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Figure  4.10:  PID  controller  bench  test  with  Figure  4.11:  PID  controller  bench  test  with 
kp  =  0.6,  hi  =  0,  k,:  =  0.  kp  =  2.  k{  =  0,  kj  =  0. 
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Figure  4.12:  PID  controller  bench  test  with 
JL  =  2.  h  =  0.01.  k *  =  0. 


Figure  4.13:  .PID  controller  bench  test  with 
JL  =  2.  h  =  0.02,  kd  =  ( 
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As  a  final  note,  unstability  can  be  caused  in  the  controller  system  if  the  gains  are  not 
chosen  properly.  Figure  4.14  shows  that  the  system  is  at  the  edge  of  going  unstable  for 
the  gains  of  kp  =  1.0,  kj  =  0.2  and  kd  =  0.0.  As  can  be  seen  in  the  figure,  the  control 
signal  oscillates  with  a  large  amplitude  which  leads  to  the  oscillation  in  the  wake  signal 
as  well.  Therefore,  such  values  for  gains  should  be  avoided. 


Figure  4.14:  PID  controller  bench  test  with  kp  —  1,  kt  =  0.2,  kd  =  0. 


Through  extensive  testing  on  the  bench,  the  gain  values  in  the  PID  controller  of  kp  = 
2.0,  kv  —  0.02  and  kd  —  0.0  are  selected  for  the  simulator  testing  in  the  anechoic 
chamber  to  achieve  “satisfactory"  performance  in  terms  of  both  transient  response  and 
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convergence  rate.  It  is  noted  that  the  chamber  test  uses  an  4140  controller  instead  of 
the  SISO  used  in  the  bench  test.  Another  difference  is  that  the  free  stream  probe  in 
the  simulator  would  always  measure  a  constant  as  of  the  ambient  pressure  according  to 
the  Bernoulli’s  law  described  in  section  4.2.  even  though  the  inlet  flow  Mach  number 
may  change  with  time. 


4.4  Active  TEB  system  performance 


The  active  flow  control  experimental  results  presented  in  this  section  are  obtained  on 
the  simulator  as  shown  in  Figure  4.1.  using  the  PID  controller  shown  in  Figure  4.8.  As 
mentioned  early,  the  purpose  of  the  PID  feedback  controller  is  to  use  the  error  signal  to 
control  the  microvalves  so  that,  as  the  flow  condition  changes  in  an  engine  the  amount 
of  wake-filling  can  be  optimized  by  adjusting  the  microvalves. 

The  first  test  performed  is  to  demonstrate  the  effectiveness  of  the  active  control  system 
to  adjust  to  varying  flow  conditions  as  the  simulator  speed  is  varied.  The  total  pressure 
of  the  four  wake  probes,  free  stream  probe  and  four  control  signals  are  recorded.  Figure 
4.15  shows  the  time  histories  of  these  signals  for  IGV  A  (IGVs  A,  B,  C,  D  are  defined 
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in  Figure  4.1).  The  lower  two  traces  in  this  figure  are  the  pressures  measured  in  the 
free  stream  and  in  the  wake  of  IGV  A.  The  difference  between  these  signals  is  the  error 
signal  to  be  minimized.  Figure  4.15  clearly  shows  that  the  controller  is  successful  in 
adaptively  driving  the  microvalves  so  that  the  wake  and  free  stream  essentially  have 
the  same  velocity,  thus  ensuring  optimal  wake-filling.  The  dynamics  of  the  controlled 
system  is  revealed  only  if  a  drastic  change  in  the  simulator  speed  is  produced,  e.g.  speed 
changes  at  t=150  and  200  seconds.  For  example,  the  simulator  speed  is  decreased  from 
40,000  rpm  down  to  20,000  rpm  in  about  9  seconds  at  t= 150  seconds.  From  the  free 
stream  and  wake  traces,  the  control  system  requires  about  12  seconds  to  optimize  the 
wake  filling,  i.e.  minimize  the  error  signal.  Similar  behavior  is  observed  at  t=200 
seconds  where  the  simulator  speed  is  ramped  up  from  20.000  to  40.000  rpm.  These 
results  demonstrate  that  the  active  flow  control  of  wakes  is  feasible  under  realistic 
changes  in  flow  conditions. 

The  total  pressure  signals  of  the  four  wake  probes  are  compared  in  Figure  4.1G  for  the 
time  range  of  140  -  220  seconds.  It  indicates  that  the  wake  dynamics  behave  differently 
when  the  simulator  speed  changes  at  t=150  and  200  seconds.  The  wake  of  IGV  A  shows 
the  largest  peak  among  the  four  while  the  wake  of  IGV  C  is  the  smallest.  The  trace  for 
IGV  B  has  the  longest  re-convergence  time  while  the  one  for  IGV  D  has  the  shortest. 
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Figure  4.15:  The  time  histories  of  the  signals  on  IGY  A. 

The  control  signals  for  the  four  IGYs  are  illustrated  in  Figure  4.17.  The  control  signal 
for  IGY  A  has  the  largest  variation  at  t=150  seconds  which  corresponds  to  the  fully 
open  condition  in  MEMS  based  microvalves  according  to  the  performance  curve  shown 
in  Figure  2.10.  The  fullv-open  condition  “holds  ”  there  since  the  system  response  is  fast 
enough  when  the  simulator  speed  decreases  abruptly,  which  leads  to  the  largest  peak 
at  t,=150  seconds  in  IGV  A  wake  total  pressure  as  seen  in  Figure  4.16.  Similarly,  the 
control  signal  for  IGY  C  has  the  smallest  variation  which  corresponds  to  the  smallest 
peak  in  the  wake  total  pressure  of  IGY  C.  In  addition,  Figure  4.17  indicates  that  the 
control  signals  for  the  four  IGYs  are  not  identical  when  the  four  wakes  are  completely 
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Figure  4.16:  The  total  pressure  of  the  four  wake  probes. 


EU 


*  i  rji 

Time  (j'ctmnls) 


Figure  4.17:  The  control  signals  for  the  four  IGVs. 
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filled. 

These  variability  in  the  results  between  IGYs  is  probably  due  to  the  physical  geometry 
in  the  IGYs  and  in  the  response  time  of  the  MEMS  based  microvalvcs  and  the  Pitot 
probes.  Thus,  these  results  suggest  that  an  independent  control  system  for  each  IGY 
could  be  important  to  compensate  for  the  variabilities  between  them. 


4.5  Acoustic  Results 


The  successful  wake-filling  leads  to  the  reduction  of  the  tonal  noise  components  due 
to  the  suppression  of  the  interaction  of  the  IGYs  and  the  rotor  blades.  Thus,  the 
attenuation  of  the  noise  at  BPF  and  its  harmonics  is  an  indication  of  the  reduction 
of  blade  vibration  and  the  HCF  damage  [1,  3].  Therefore,  the  sound  field  with  and 
without  TEB  is  monitored  at  far-held  locations. 

The  monitored  far-held  microphone  locations  are  depicted  in  Figure  4.18.  In  total, 
twelve  positions  are  measured  in  the  horizontal  plane,  from  0°  to  110°,  every  10°.  Since 
the  inlet  used  is  symmetric,  the  measured  directivity  holds  for  corresponding  positions 
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in  the  vertical  plane,  and  also  in  the  mirrored  horizontal  plane.  The  0°  position  lines 
up  with  the  axis  of  the  simulator.  Since  acoustic  data  are  inherently  unsteady  due  to 
the  unsteady  IGV-rotcr  interaction,  the  sound  pressure  level  results  at  each  position 
are  obtained  by  taking  an  average  of  50  spectra. 


Simulator 


Side  View 


Front  View 


Figure  4.18:  Position  of  far-held  monitoring  microphones. 


Figure  4.19  shows  the  sample  far-held  acoustic  spectra  acquired  at  the  20°  microphone 
location  at  a  simulator  speed  of  30.000  rpm.  The  corresponding  BPF  tone  is  at  9.000 
Hz.  These  spectra  demonstrate  a  reduction  of  about  8.2  dB  at  the  BPF  tone  and 


reduction  between  6  to  3  dB  at  the  first  four  harmonics.  The  wake  filling  does  not  lead 
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to  reduction  of  the  broadband  noise.  This  is  because  the  broadband  noise  is  generated 
by  the  interaction  of  random  disturbance  while  the  tonal  noise  is  caused  by  the  IGY 
mean  wakes  [24].  These  tonal  noise  reductions  again  clearly  indicate  the  successful 
wake  filling  as  well  as  the  reduction  of  the  blade  forcing  function  and  vibration.  It  is 
interesting  to  point  out  that  the  most  objectionable  noise  is  from  the  tones  and  the 
tonal  reductions  are  clearly  perceived  by  the  human  ear. 


Figure  4.19:  Far-field  acoustic  spectrum  at  20°  with  and  without  active  control. 
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Table  4.1:  Confidence  interval  analysis  at  20°  (ju=mean.  r=standard  deviation). 


SPL  (dB) 

BPF  tone 

2BPF  tone 

Without  TEB 

With  TEB 

Without  TEB 

With  TEB 

p  +  2r 

95.6 

87.6 

94.3 

88.4 

h 

94.8 

86.6 

93.7 

87.6 

p  -  2r 

94.0 

85.5 

93.0 

86.8 

An  uncertainty  study  was  performed  to  determine  the  confidence  in  the  experimental 
results.  The  mean  value  p  and  standard  deviation  r  for  the  measurements  at  the  BPF 
and  2BPF  tones  were  computed  without  and  with  active  flow  control.  Using  these 
values,  the  confidence  analysis  is  presented  in  Table  4.1  as  well  as  shown  in  Figure 
4.19.  It  is  noted  that  the  confidence  level  of  this  range  from  p  —  2 r  to  p  +  2 t  is  94 if 
the  probability  density  of  measurement  data  follows  the  Gaussian  distribution,  which 
is  the  case  for  many  engineering  applications  [59].  Therefore,  it  can  be  concluded  that 
the  amount  of  noise  reductions  at  BPF  and  harmonics  is  very  reliable. 
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4.5.1  Sound  radiation  pattern  at  30,000  rpm 

The  radiating  directivities  are  also  recorded  for  the  purpose  of  estimating  the  sound 
power  level.  Figure  4.20  compares  the  directivity  plots  at  a  simulator  speed  of  30.000 
rpm  without  and  with  trailing  edge  blowing.  The  data  presented  consists  of  the  sound 
pressure  levels  at  the  blade  passing  frequency  and  the  first  harmonic,  and  overall  sound 
pressure  level.  Figure  4.20(a)  shows  the  directivity  pattern  at  BPF.  The  confidence 
interval,  i.e.,  mean  ±  standard  deviation,  is  also  shown  as  dotted  lines.  As  observed, 
the  BPF  tone  is  reduced  at  all  recorded  positions  except  at  0°,  which  lines  up  with  the 
axis  of  the  simulator.  Significant  reductions,  greater  than  3  dB,  are  obtained  at  most 
positions.  The  maximum  reduction  is  8.2  dB  at  the  20°  position  which  agrees  well  with 
the  7  dB  reduction  obtained  in  the  work  by  Leitcli  [25.  30].  Note  that  aerodynamic 
measurements  by  Leitch  showed  nearly  full  filling  of  the  wake  in  his  experiments.  Since 
the  same  setup  is  used,  it  is  quite  reasonable  to  conclude  that  nearly  complete  filling 
of  the  wake  is  also  achieved  here.  At  this  speed,  the  principal  radiation  angles  of  the 
modes  cut-on  at  the  BPF  due  to  IGY-rotor  interaction  are  also  calculated  according 
to  Eqn.  (3.26).  The  rotor-EGY  interaction  does  not  have  any  cut-on  modes  at  BPF. 
From  Table  3.6,  the  cut-on  modes  due  to  IGY-rotor  interaction  at  the  speed  of  30,000 
rpm  are  (±2,0),  (±2, 1)  and  (±6,0)  whose  principal  radiation  angles  are  computed  to 
be  ±19°.  ±48°  and  ±61°,  respectively.  The  maximum  reduction  angle  observed  in  the 
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Figure  4.20:  Directivity  of  BPF  tone  at  30,000  rpm. 
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experiments  is  20°  which  agrees  very  well  with  the  predicted  19°  for  mode  (2.0).  This 
implies  that  mode  (2.0)  is  greatly  reduced  due  to  the  trailing  edge  blowing. 


Figure  4.20(b)  shows  the  radiation  directivity  at  the  first  harmonic,  i.c..  2BPF.  In 
this  case  the  reductions  are  confined  to  the  sector  from  20°  to  50°,  inclusive.  The 
maximum  reduction  of  6.1  dB  also  occurs  at  the  20°  position.  The  pressure  levels 
observed  outside  of  this  sector  are  increased.  As  predicted  in  Table  3.6.  there  are  much 
more  propagating  modes  present  in  the  inlet.  The  far- field  directivity  depends  on  the 
combined  radiation  of  all  the  generated  modes,  including  interference  effects.  Trailing 
edge  blowing  reduces  the  amplitude  of  modes  propagating  due  to  IGV-rotor  interaction. 
The  destructive  interference  of  these  modes  with  modes  generated  by  other  mechanisms 
is  affected  by  this  amplitude  reduction.  This  would  cause  the  sound  pressure  level  to 
increase  along  certain  directions. 


4.5.2  Sound  radiation  pattern  at  40,000  rpm 


Figure  4.21  presents  the  acoustic  data  recorded  at  a  simulator  speed  of  40.000  rpm.  As 
shown  in  Figure  4.21(a),  the  BPF  tone  shows  reduction  in  the  sector  0°  to  50°  which 
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agree  with  the  calculated  principal  lobe  angles  due  IGY-rotor  interaction  of  14°.  35°. 
40°,  51°  and  78°.  It  implies  that  the  far-field  sound  level  at  BPF  is  mainly  contributed 
by  IGV- interaction  as  compared  to  other  interactions.  The  average  reduction  over  all 
measured  positions  is  2.7  dB,  with  a  maximum  of  7.3  dB  at  20°.  In  the  sector  60°  to 
110"  it  can  be  observed  that  the  sound  pressure  level  with  and  without  trailing  edge 
blowing  remains  the  same. 


The  first  harmonic  levels  exhibit  reductions  at  most  positions,  as  shown  in  Figure 
4.21(b).  The  sound  pressure  level  at  20°  increases  by  3.2  dB  after  trailing  edge  blowing. 
This  could  be  due  to  unmasking  of  modes  due  to  other  interaction  mechanisms  such 
as  core  flow  distortion  in  the  inlet.  The  peak  reduction  is  3.1  dB  at  30°. 

Trailing  edge  blowing  to  re-energize  the  wake"-  from  the  IGVs  results  in  the  suppression 
of  the  IGY-rotor  interaction  effect.  However,  the  far-field  radiation  is  due  to  the 
contribution  of  several  noise  sources,  in  particular  the  IGY-rotor,  rotor-rotor,  rotor- 
EGY  interactions  and  inlet  flow  distortions  as  a  result  of  the  lack  of  an  inlet  flow 
control  device.  Thus,  the  maximum  allowed  reduction  in  the  radiation  directivity 
from  re-energizing  the  IGV  wakes  is  limited  by  the  presence  of  the  rotor-rotor,  rotor- 
EGY  interaction  and  inlet  flow  distortion.  In  the  experiments  performed  to  obtain 


the  radiation  directivities,  the  wakes  were  effectively  re-energized,  i.e.,  the  wake  total 
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Figure  4.21:  Directivity  of  BPF  tone  at  40,000  rpm. 
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pressure  was  equal  to  the  free  stream  total  pressure  (see  Figure  4.15).  Because  of 
this,  it  is  believed  that  the  radiation  directivities  with  trailing  edge  blowing  shown  in 
Figure  4.21  are  due  from  both  the  rotor-rotor  and  rotor-EG\  interactions,  and  inlet 
flow  distortion.  In  addition,  since  all  the  modes  excited  by  the  rotor-EG\  interaction 
are  cut-off  at  a  simulator  speed  of  30.000  rpm  as  indicated  in  Table  3.6,  the  radiation 
directivities  with  trailing  edge  blowing  shown  in  Figure  4.20  are  due  only  from  the 
rotor-rotor  interaction  and  inlet  flow  distortion. 


4.5.3  Sound  power  level  reductions  at  30,000  rpm  and  40,000 
rpm 


Determining  the  sound  power  has  the  advantage  that  it  describes  the  sound  energy 
output  of  the  source,  without  regard  to  the  directivity.  Therefore,  the  sound  power 
reduction  may  provide  a  better  indication  of  the  effect  of  trailing  edge  blowing  on 
reducing  the  sound  source  strength.  This  analysis  provides  an  estimate  of  the  change 
in  sound  power  level  with  trailing  edge  blowing  helps  better  understand  the  effect 
trailing  edge  blowing  has  on  the  unsteady  interaction.  The  radiated  sound  power, 
W,  is  computed  from  the  radiation  directivities  by  integrating  the  far-field  acoustic 
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w  =  Y*™A  (4.12) 

U  pc 

where  p,  rms  is  the  effective  sound  pressure  measured  at  the  ith  microphone  location: 
Ai  a  surface  area  where  the  sound  pressure  is  assumed  to  be  uniform  and  represented 
by  the  measurement  at  the  zt,(  microphone  location.  In  the  calculation,  a  surface  area 
.4;  =  2nR2sin6  ■  vr/18  is  used  where  the  radius  R  is  shown  in  Figure  4.18  [40]. 


The  BPF  power  reduction  is  significant  at  both  test  speeds  as  shown  in  Table  4.2.  At 
30,000  rpm,  the  power  in  the  BPF  tone  was  reduced  by  4.4  dB.  At  40.000  rpm  the 
reduction  in  the  BPF  tone  power  was  2.9  dB.  A  reduction  in  far-field  acoustic  power 
can  be  directly  interpreted  as  a  reduction  in  the  source  power.  Trailing  edge  blowing 
reduces  the  magnitude  of  the  velocity  defect  seen  by  the  rotor  and  also  smooths  the 
velocity  profile.  In  doing  so,  this  decreases  the  unsteady  loading  experienced  by  the 
blade  as  it  passes  through  the  wake,  and  hence  the  reduction  in  source  strength. 


At  the  same  speed  the  blade  passing  tone  power  is  reduced  more  than  the  first  harmonic 
tone.  This  is  probably  because  the  wake-filling  is  not  perfect,  i.e..  "jet-like"  or  "double¬ 
humped”  which  will  be  discussed  in  next  chapter.  The  “jet-like”  wake-filling  does  not 
lead  to  much  reduction  in  the  harmonics  because  the  harmonic  loading  on  rotor  blades. 
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Tabic  4.2:  Sound  power  level  reduction. 


30.000  rpm 

40,000  rpm 

BPF 

First  Harmonic 

BPF 

First  harmonic 

Before  active  control 

101.4  dB 

96.2  dB 

111.6  dB 

108.0  dB 

After  active  control 

97.0  dB 

94.2  dB 

108.7  dB 

107.0  dB 

Power  level  reduction 

4.4  dB 

2.0  dB 

2.9  dB 

1.0  dB 

obtained  through  transforming  the  jet-like  forcing  function  into  the  frequency  domain 
using  FFT,  do  not  alter  much  as  compared  to  pure  wake  forcing  function.  The  lesser 
power  reduction  at  the  higher  speed  is  probably  because  with  an  increase  in  speed  the 
acoustic  power  generated  by  the  rotor-EGV  interaction  increases  as  well  as  inlet  flow 
distortions. 


4.6  Summary 


An  active  flow  control  system  with  the  Pitot  probes  as  the  wake  sensors  is  success¬ 
fully  developed  to  obtain  effective  wake  management  through  trailing  edge  blowing  in 
a  rotating  environment.  Although  Pitot  probes  are  intrusive  sensors  which  are  not 
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appropriate  for  realistic  engine  systems,  use  of  pitot  probes  is  a  preliminary  step  to¬ 
ward  more  in-depth  for  TEB  control  research.  More  importantly,  results  obtained  from 
this  system  can  serve  as  the  baseline  in  evaluations  of  other  possible  control  systems 
indirect  wake  sensing. 

A  4-channel  PID  controller  generates  the  signal  voltage  that  controls  the  flow  rate 
through  each  microvalve.  The  inputs  to  the  PID  controller  are  the  total  pressures 
in  the  free  stream  and  wakes  behind  each  IGV.  The  error  signal  is  computed  as  the 
difference  between  these  signals,  with  the  free  stream  total  pressure  as  the  reference. 
Pitot  probes  immersed  in  the  flow  measure  the  wake  total  pressures  that  are  used  as 
inputs  to  the  controller.  The  advantages  of  this  Pitot-probe  based  system  are  that 
effective  wake-filling  is  achieved  and  changes  in  operating  condition  can  be  tracked  for 
effective  wake  management  with  changing  conditions. 


From  the  results  presented  it  can  be  concluded  that  the  Pitot-probe  based  active  TEB 
control  system  produces  a  very  uniform  flow  field  downstream  of  the  IGV.  This  is 
because  a  Pitot  probe  measures  the  direct  wake  signal  which  is  exactly  what  the  con¬ 
troller  seeks  to  minimize.  The  ability  of  the  system  to  achieve  optimum  wake  filling 
when  subjected  to  a  change  in  inlet  flow  conditions  demonstrates  the  feasibility  and 
advantage  of  active  flow  control.  It  was  also  demonstrated  that  far-field  sound  pressure 
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levels  at  the  relevant  discrete  frequencies  arc  reduced  by  trailing  edge  blowing.  Fur¬ 
thermore.  the  reduction  in  sound  source  strength  supports  the  effectiveness  of  trailing 
edge  blowing  in  reducing  unsteady  IGY-rotor  interaction.  The  maximum  sound  level 
reduction  is  observed  at  the  blade  passing  frequency  of  each  tested  speed.  The  maxi¬ 
mum  tone  reductions  obtained  are  around  8.2  dB  at  30.000  rpm  and  7.3  dB  at  40.000 
rpm.  Sound  power  level  at  the  BPF  tone,  calculated  from  the  measured  directivity,  is 
reduced  by  4.4  dB  at  30,000  rpm  and  2.9  dB  at  40,000  rpm. 


However,  as  mentioned  early  in  this  chapter,  the  drawbacks  of  using  Pitot  probes  are 
obvious.  This  is  mainly  because  Pitot-probe  based  TEB  control  is  not  appropriate  for 
the  real  engine  environments.  Therefore,  there  exists  a  demand  to  develop  an  advanced 
TEB  control  system  which  does  not  require  Pitot  probes.  In  the  next  chapter,  a  non- 
intrusive  control  system  based  on  microphones  is  introduced. 


Chapter  5 


Microphone-based  Flow  Control 
Systems 


New  flow  control  systems  based  on  the  microphone  sensors  will  be  described  in  this 
chapter.  The  non-intrusive  wake  sensing  approach  enables  the  active  flow  control 
system  to  eliminate  the  need  of  intrusive  Pitot  probes  as  presented  in  Chapter  4.  The 
significance  of  using  acoustic  sensing  lies  in  that  it  provides  a  more  feasible  approach 
toward  the  realistic  engine  TEB  control  applications. 
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In  section  5.1,  the  motivation  of  the  acoustic  sensing  approach  will  bo  briefly  described. 
The  experimental  setup  using  the  microphone  sensors  follows.  A  unique  signal  pro¬ 
cessing  device  will  be  discussed  in  section  5.3.  Two  control  schemes  will  be  developed 
in  section  5.4  and  5.5,  respectively  as  well  as  the  results  obtained  from  both  control 
schemes. 


5.1  Introduction 


The  experimental  tests  in  Chapter  4  demonstrated  the  potential  of  using  Pitot  probes 
as  the  wake  sensors  for  optimum  re-energizing  of  the  wakes.  However,  there  exist 
several  drawbacks  of  using  Pitot  probes.  First,  such  probe  is  not  practical  to  be  used 
in  the  true  engine  environments.  Extreme  care  needs  to  be  taken  for  the  sensor  to 
be  mounted  into  the  right  position  and  angle.  Secondly,  Pitot  probes  are  intrusive 
sensors  which  themselves  can  create  the  wakes.  Lastly,  the  Pitot  probes  can  only  be 
used  to  sense  the  wake  of  non-turning  IGVs.  This  is  because  Pitot  probe  can  only 
sense  correctly  the  flow  velocity  at  a  fixed  angle  (i.e.,  pointing  along  a  streamline). 
However,  the  wake  region  of  a  curved  IGV  turns  pitchwise  as  the  engine  operating 
condition  changes,  i.e.  when  inlet  Mach  number  changes.  Curved  IGVs  are  widely 


Jimvei  Feng 


Chapter  5.  Microphone-based  Flow  Control  Systems 


108 


used  in  aircraft  engines.  Thus,  there  is  a  need  to  develop  non-intrusivc  wake  sensing 
strategies  to  infer  the  filling  of  the  wakes. 

In  this  chapter,  an  acoustic  sensing  approach  using  microphones  is  proposed  to  indi¬ 
rectly  sense  the  wake-filling.  The  basic  idea  can  be  obtained  from  Figure  3.1.  The 
figure  show's  that  the  upstream  IGV  wake  is  a  common  forcing  function  to  tonal  fan 
noise  and  blade  vibrations.  Therefore,  when  the  wake  is  fully  re-energized  by  TEB, 
both  blade  vibration  and  tonal  fan  noise  are  minimized.  It  follows  logically  that  a 
reduction  of  tonal  fan  noise  by  TEB  implies  a  reduction  of  fan  blade  vibrations  as 
well.  From  this  point  of  view,  microphones  used  for  acoustic  sensing  can  be  viewed  as 
the  indirect  “wake  sensors”  or  “vibration  sensors”  and  thus  offers  the  potential  to  be 
integrated  with  the  other  elements  of  TEB  control  system. 

Two  TEB  control  schemes  in  conjunction  with  the  microphone  sensing  approach  will 
be  presented  in  this  chapter.  In  the  first  control  scheme  only  the  magnitude  of  the 
BPF  tone  is  used  to  generate  the  error  signal  for  the  controller.  In  this  case,  an 
optimal  controller  is  implemented.  The  second  control  scheme  consisted  of  using  both 
the  magnitude  and  phase  of  the  BPF  tone  sensed  by  the  microphones.  In  this  second 
approach,  a  PID  feedback  controller  is  implemented.  These  control  approaches  arc 
next  described  and  their  performance  w'ere  tested. 
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5.2  Experimental  Setup 


Figure  5.1  shows  a  schematic  of  the  experimental  setup  with  the  microphones  as  the 
indirect  wake  sensors.  Compared  to  Figure  4.1,  its  shows  that  the  Pitot  probes  are 
removed  while  the  microphones  are  flush  mounted  on  the  inlet  wall.  The  basic  control 
idea  is  that  the  microphones  sense  the  noise  induced  by  the  IGY-rotor  interaction.  This 
information  is  then  provided  to  a  digital  controller  that  generates  the  control  signals  for 
the  MEMS  based  microvalves  to  reduce  the  noise  from  this  interaction.  The  reduction 
of  the  noise  signature  results  in  effective  wake-filling. 


Fan  rotor 


Figure  5.1:  Microphone-based  experimental  setup. 
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The  microphones  used  in  the  experiments  are  high-intensity  microphones  (model  10GB) 
from  PCB  Electronics,  as  shown  in  Figure  5.2(a).  This  microphone  is  specially  designed 
to  measure  pressure  perturbations  in  air  or  in  fluids  in  severe  environments.  They 
can  also  be  used  to  measure  very  small  pressure  disturbances  on  a  much  higher  static 
head,  with  certain  precautions.  Such  microphone  features  high-sensitivity,  acceleration- 
compensated  quartz  pressure  elements  coupled  to  built-in  integrated  circuit  impedance 
converting  amplifier.  The  quartz  elements  in  model  10GB  microphone  utilizes  a 
special  cut  in  quartz  to  produce  a  proportionately  higher  output  voltage  than  the 
standard  X-cut  compression  crystal  normally  used.  A  built-in  seismic  mass  acting  on 
another  quartz  crystal  effectively  cancels  the  spurious  signal  produced  by  the  mass  of 
the  diagram  and  end  piece  acting  upon  the  very  sensitive  crystals  in  the  presence  of 
axial  vibration  inputs.  This  design  produces  an  extremely  high  level  output  signal  with 
good  resolution,  relatively  free  from  unwanted  vibration  effects.  Figure  5.2(b)  shows 
the  cross  section  of  106B  microphone  elements.  As  with  all  quartz  sensors,  the  high 
rigidity  results  in  negligible  diaphragm  producing  excellent  linearity  from  the  threshold 
pressure  to  full-scale  pressure  as  seen  in  Figure  5.2(c).  The  high  rigidity  of  quartz  also 
results  in  sensors  with  high  natural  frequency,  giving  a  very  wide  useful  frequency 
range.  The  built-in  electronics  consists  of  a  voltage  follower  with  unity  gain.  A  single 
wire  feeds  constant  current  power  to  the  voltage  follower  and  also  carries  the  dynamic 
signal.  Caution  needs  to  be  taken  never  to  apply  power  directly  to  the  microphone 
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Figure  5.2:  The  model  106B  microphone  from  PCB  Electronics,  (a)  microphone  view; 
(b)  microphone  cross-section:  and  (c)  performance  curve. 


Jimvei  Feng 


Chapter  5.  Microphone- based  Flow  Control  Systems 


112 


without,  the  current-limiting  protection  (2  mA  maximum),  as  to  do  so  will  destroy  the 
built-in  amplifier. 


In  the  experimental  setup,  four  high  intensity  microphones  are  flush  mounted  on  the 
inlet  wall  right  behind  the  four  IGVs,  named  IGY  A,  B,  C  and  D,  respectively  as 
shown  in  Figure  5.3.  Accordingly,  the  four  microphones  behind  IGV  .4,  B.  C,  and 
D  are  referred  to  as  microphone  1,  2,  3,  and  4,  respectively  for  convenience.  The 
microphone  3  can  be  seen  to  be  flush  mounted  behind  IGY  C.  As  stated  before,  the 
four  wakes  behind  the  IGYs  become  the  four  noise  sources  at  the  fan  face  when  they 
arc  chopped  through  the  rotor.  Placing  the  microphones  in  between  the  IGYs  and  the 
rotor  makes  them  to  sense  these  noise  sources  effectively  at  BPF  and  its  harmonics.  In 
the  experiments  here,  the  simulator  is  operated  at  a  speed  of  30.000  rpm.  which  leads 
to  a  BPF  tone  frequency  of  9,000  Hz. 


5.3  Adaptive  Extraction  of  BPF  Tone 


The  implementation  of  the  acoustic  sensing  strategy  requires  some  unique  signal  pro¬ 
cessing.  The  simulator  radiated  noise  spectrum  consists  of  a  set  of  dominant  tones  and 
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Flush  mounted  microphone 


Figure  5.3:  Simulator  inlet  with  four  casing  microphones  behind  four  IGYs. 

a  broadband  component  as  can  be  seen  in  Figure  4.19.  One  of  the  sources  of  the  tonal 
components  at  the  BPF  and  harmonics  are  due  to  the  IGV-rotor  interaction,  i.e..  noise 
due  to  the  wakes  Thus,  the  tonal  component  information  instead  of  broadband  noise 
needs  to  be  extracted  out  for  use  in  the  active  flow  control  system. 


Furthermore,  the  BPF  tone  is  selected  to  construct  the  error  signal  for  the  control  sys¬ 
tem.  This  is  primarily  because  BPF  tone  is  expected  to  dominate  the  tonal  component 
according  to  the  fan  noise  generation  mechanism.  (The  fan  noise  generation  will  be 
discussed  in  detail  in  Appendix  C).  The  basic  concept  of  the  mechanism  is  given  briefly 
here  as  follows.  Since  the  fan  noise  is  generated  as  the  rotor  blades  chop  through  the 
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upstream  IGV  wakes,  the  upstream  wake  can  be  viewed  as  the  forcing  function  ap¬ 
plied  to  the  rotor.  A  possible  shape  of  the  forcing  function.  f(t).  can  be  illustrated 
in  Figure  5.4  with  a  period  of  1  j  /bpf  since  the  blades  chop  through  the  wake  at  the 
frequency  of  BPF.  The  variation  of  the  trace  is  due  to  the  flow  variation  in  the  wake. 
The  generated  sound  energy  distribution  at  BPF  and  its  harm'onics  is  determined  by 
the  frequency  distribution  of  the  forcing  function,  which  can  be  computed  by  using 
Fourier  scries  expansion  with  the  fundamental  frequency  of  Jeff-  The  regular  Fourier 
series  expansion  is  expressed  as: 

'  OC 

■/(f)  =  Y.  F^bpf)  ■  einfBPFt  (5.1) 

n=~o o 

where  F(nfBPF)  is  the  complex  Fourier  coefficient;  n  denotes  the  harmonics  index. 
According  to  the  periodical  shape  shown  in  Figure  5.4,  the  largest  Fourier  coefficient 
would  occur  at  the  fundamental  frequency,  Jbpf-  Therefore,  the  generated  noise  at 
BPF  is  dominant.  Another  reason  of  choosing  BPF  tone,  instead  of  its  harmonics,  to 
construct  the  error  signal  is  because  high  order  harmonics  would  not  be  expected  to 
be  reduced  if  a  momcntumlcss  wake  is  filled  “jet-like"  or  called  double-hump  as  shown 
in  Figure  5.5.  “Jet-like”  wake-filling  can  occur  easily  in  practice  since  the  TEB  jet  can 
not  be  spread  out  adequately  downstream  IGV. 


The  most  challenging  part  in  obtaining  BPF  tone  from  the  broadband  spectrum  is  the 
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Figure  5.4:  The  shape  of  the  forcing  function  applied  to  the  rotor. 

/<0 


T  ime  (seconds) 

Figure  5.5:  The  shape  of  the  forcing  function  with  TEB. 

requirement  of  adaptivity.  This  is  because  as  the  engine  speed  changes,  the  BPF  will 
also  vary,  which  makes  a  conventional  bandpass  filter  with  a  fixed  center  frequency 
insufficient.  Other  requirements  includes  the  reliability,  insensitivity  to  temperature, 
precision  and  so  forth.  The  conventional  approach  is  using  DSP-based  algorithm,  i.e. 
the  adaptive  filter  with  LMS  (least  mean  square)  algorithm  [GO].  This  method  requires 
a  complicated  hardware  system  which  includes  ADCs  to  convert  the  continuous  signals 
into  digital  signals,  a  DSP  with  adaptive  filters,  and  DACs  to  convert  the  digital  signals 


Jimvei  Feng  Chapter  5.  Microphone-based  Flow  Control  Systems  116 

back  to  the  continuous  domain. 

To  this  end,  a  novel  adaptive  bandpass  filter  device  is  developed  and  tested  to  capture 
the  BPF  tone  component  from  the  acoustic  signatures  of  the  simulator.  Compared 
to  the  conventional  DSP-based  methods  [60],  this  extractor  pertains  some  attractive 
properties  which  includes  the  high  performance  in  terms  of  reliability,  steady-state 
error,  respond  time,  precision  and  sensitivity  to  temperature.  In  addition,  it  is  easy  to 
implement  since  there  is  no  need  for  DSP  and  their  programming. 


f  :  * 

I  uni?  irarimri*} 


Figure  5.6:  The  block  diagram  for  the  adaptive  BPF  tone  extractor. 
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Figure  5.6  shows  the  block  diagram  of  the  device  which  consists  of  two  important 
components.  The  BPF  tone  is  obtained  by  using  an  analog.  high-Q  bandpass  switched- 
capacitor  filter  (SGF)  whose  center  frequency  adjusted  using  the  simulator  shaft  signal, 
i.e.  referred  here  as  the  reference  signal  or  RPM  signal.  Another  component  is  the 
phase-locked-loop  (PLL)  which  is  used  to  condition  the  rpm  signal.  The  idea  is  that 
the  reference  signal,  after  conditioning,  is  input  into  a  PLL  circuit.  This  PLL  circuit 
is  designed  to  generate  the  clock  signal  for  the  SCF.  The  clock  signal  changes  with  the 
reference  signal,  i.e.  the  clock  signal  tracks  the  reference  signal.  The  center  frequency 
of  the  bandpass  SCF  is  determined  by  the  clock  signal  which  in  turn  is  controlled  by 
the  reference  signal  through  the  PLL  circuit.  Thus  tlie  filter  center  frequency  “adapts” 
or  “tracks”  the  reference  signal.  The  RMS  (root  mean  square)  circuit  follows  the  SCF 
to  compute  the  BPF  tone  power  or  magnitude.  The  various  time  domain  signal  in  the 
device  arc  also  shown  in  Figure  5.6.  As  can  be  seen,  the  BPF  tone  embedded  in  the 
microphone  output  signal  is  converted  by  the  device  to  a  DC  value  which  will  be  input 
to  a  digital  controller. 

Figure  5.7  shows  its  excellent  static  performance  in  frequency  domain  (Figure  5.7a) 
and  in  time  domain  (Figure  5.7a)  using  the  experimental  data  obtained  from  simulator 
testing.  The  unfiltered  signal  which  consists  of  tonal  component  and  broadband  back¬ 
ground  noise  is  from  one  of  the  casing  microphones  when  the  simulator  is  running  at 
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Figure  5.7:  The  static  perfc 
domain;  (b)  in  time  domain. 
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around  30,000  rpm.  Figure  5.7a  indicates  that  only  the  BPF  signal  is  filtered  out  while 
the  higher  harmonics  and  broadband  background  noise  are  rejected.  This  comparison 
can  also  be  seen  in  time  domain  in  Figure  5.7b  which  shows  a  clean  sinusoidal  wave  and 
original  microphone  output  signal.  It  is  worth  to  point  out  that  no  low  pass  filters  are 
needed  for  the  original  microphone  signal  in  order  to  prevent  frequency  aliasing  before 
it  is  fed  to  the  SCF.  This  is  because  the  clock  frequency  for  the  SCF  is  much  higher 
than  the  upper-limit  frequency  in  the  frequency  response  of  a  casing  microphone.  Low- 
pass  filters  are  often  required  in  a  conventional  data  acquisition  system.  In  addition, 
higher  Q-factor  for  this  adaptive  filter  can  be  obtained  through  a  serial  connection 
of  multiple  SCFs.  Four  SCFs  are  cascaded  in  the  device  for  each  microphone  output 
signal. 

The  dynamic  performance  of  the  device  mainly  depends  on  the  dynamic  performance 
of  the  PLL.  The  major  performance  parameters  of  a  PLL  includes  the  tracking  speed, 
hold-on  range,  lock-in  range  and  pull-in  range.  Appendix  A  has  more  details  about 
these  parameters  and  how  to  obtain  them.  In  this  case,  the  two  most  important 
parameters  are  hold-on  range  and  tracking  speed.  The  hold-on  range  Q,  which  is  the 
range  of  frequency  over  a  loop  will  remain  locked,  is  evaluated  for  this  extractor  to  be 
20,000  rpm  —  81,000  rpm  in  terms  of  simulator  speed.  Note  this  range  is  broad  enough 
to  cover  the  operating  range  of  the  simulator  which  is  20,000  rpm  —  40,000  rpm.  The 
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tracking  speed,  when  evaluated  in  terms  of  setting  time  ts  in  the  step  response,  is  about 
1  //Second.  This  means  the  center  frequency  of  the  adaptive  filter  can  adapt  to  the 
input  signal  in  1  //Second  in  response  to  a  unit  step  in  the  simulator  speed.  It  is  noted 
that  this  setting  time  is  short  enough  for  our  application  case. 


5.4  Acoustic  Sensing  Approach  Using  BPF  Magni¬ 
tude 


5.4.1  Control  algorithm 


As  mentioned  earlier,  the  IGY  wake  -  rotor  interaction  effect  results  in  strong  tones  at 
BPF  and  harmonics.  The  re-energizing  of  the  wake  by  TEB  leads  to  noise  attenuation 
due  to  the  reduction  of  IGV-rotor  interaction.  In  addition,  the  optimum  wake-filling 
minimizes  the  BPF  tone  acoustic  signal.  Thus,  the  power  of  the  BPF  tone  sensed 
by  the  casing  microphones  as  a  function  of  the  trailing  edge  blowing  rate  will  have  a 
minimum  at  the  optimum  blowing  as  shown  in  Figure  5.8.  It  is  important  to  reiterate 
that  the  acoustic  signal  will  not  be  driven  to  zero  because  of  two  factors.  The  first 
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factor  is  that  the  IGY-rotor  interaction  is  not  the  only  noise  source  at  the  BPF  tone 
and  harmonics.  The  other  noise  sources  are  the  rotor-rotor.  rotor-ECA  interactions  as 
well  as  flow  distortions  ingested  into  the  inlet  due  to  the  lack  of  a  flow  inlet  control 
device  (ICD).  Though  the  acoustic  modes  induced  by  the  rotor-rotor  and  the  rotor- 
EGY  interactions  are  cut-off  at  the  simulator  speeds  tested,  ingested  distortions  will 
not  be  eliminated  by  the  TEB.  The  second  factor  that  the  BPF  tone  is  not  completely 
eliminated  is  simply  because  the  wake  re-energizing  process  is  not  perfect,  i.e.  there  is 
some  wake  deficit  left  after  TEB. 


Figure  5.8:  The  BPF  power  in  casing  microphones  as  a  function  of  TEB  rate. 


An  experiment  is  performed  on  the  simulator  in  the  anechoic  chamber  to  verify  this 
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behavior  in  the  inlet  microphones  before  the  active  control  scheme  is  developed  and 
tested.  This  is  conducted  by  measuring  the  SPL  of  the  BPF  tone  in  both  far-field  and 
inlet  microphones  as  a  function  of  TEB  rate.  The  results  for  the  simulator  speed  of 
30,000  rpm  are  shown  in  Figure  5.9  where  the  blowing  rate  sweeps  from  no-blowing 
to  over-blowing.  The  dashline  represents  SPL  of  the  BPF  tone  sensed  in  the  far-field 
microphone  at  20°  of  the  simulator  axis.  The  solid  line  is  the  SPL  of  the  BPF  tone 
measured  by  the  casing  microphone  behind  IGV  .4  (see  Figure  5.3.  Page  113).  It  can 
be  found  from  the  figure  that  the  maximum  noise  reduction  in  both  of  the  casing 
and  far-field  microphones  occurs  at  the  same  blowing  rate.  This  indicates  that  the 
maximum  BPF  tone  reduction  in  the  casing  microphones,  which  will  be  achieved  after 
the  non-intrusive  controller  has  converged,  can  also  infer  the  best  wake-filling,  just  as 
the  maximum  BPF  tone  reduction  in  the  far-field  microphones  infers  the  best  wake¬ 
filling.  The  significance  of  this  conclusion  lies  in  that  it  provides  the  justification  for 
the  TEB  control  using  acoustic  sensing  approach. 

Therefore,  the  first  control  scheme  is  to  implement  a  simple  optimal  controller  that 
determines  the  optimum  blowing  rate  by  searching  for  the  minimum  of  the  BPF  power 
curve  in  Figure  5.8.  Figure  5.10  shows  the  block  diagram  for  the  control  system. 
The  feedback  loop  starts  from  the  four  microphones  signals,  which  are  pre-conditioned 
before  they  are  used  by  the  adaptive  BPF  extractor.  The  simulator  RPM  signal,  which 
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Figure  5.9:  Sound  pressure  level  at  BPF  tone  in  the  far-field  and  casing  microphones. 


is  obtained  from  a  magnetic  pickup  built  in  the  simulator,  is  used  to  synchronize  the 
BPF  extractor  so  that  its  center  frequency  can  track  the  BPF  signal  with  simulator 
rotating  speed.  The  output  of  BPF  extractor,  which  is  a  DC  value,  serves  as  the 
error  signal  for  die  digital  controller.  The  controller  generates  the  control  signal  for 
the  MEMS  based  microvalves.  It  is  important  to  point  out  here  that  different  offset 
voltages  should  be  used  in  the  microvalve  offset  circuits  to  compensate  for  the  variations 
between  the  four  IGVs.  According  to  the  results  shown  in  Figure  4.17  where  Pitot 
probes  were  used  as  flow  sensors,  the  offset  voltages  of  4.5,  4.0,  3.5,  3.0  Volts  are 
applied  to  IGV  C,  D,  B  and  A,  respectively.  Hence,  the  same  control  signal  can  be 
used  to  drive  all  the  microvalves  of  the  four  IGVs,  i.e.,  SISO  controller.  The  control 
algorithm  used  here  aims  to  locating  the  minimum  of  the  curve  shown  in  Figure  5.8.  In 
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the  context  of  control  theory,  this  is  a  typical  one-dimensional  optimal  control  problem: 
A  gradient  method  is  utilized  to  search  for  the  minimum  moving  along  the  opposite  of 
the  function  gradient.  That  is, 

£4+ 1  =  £4  +  Dk  ■  S  (5.2) 

where  £4  is  the  control  signal  at  time  tk.  S  is  the  step  size,  and  D *.  is  the  search 
direction  with  the  initial  value  of  1  which  is  updated  as: 

{Dk  if  e*+i  >  ek 

(5.3) 

-Dk  if  ek+i  <  ek 

where  ek  denotes  the  time-averaged  BPF  power  from  all  the  four  casing  microphones, 
at  time  tk. 

Time-averaging  is  needed  because  the  BPF  tone  signal  fluctuates  severely  due  to  the 
unsteady  behavior  of  the  rotor  interacting  with  the  distorted  incoming  flow.  The 
upstream  incoming  flow  distortion  is  probably  mainly  due  to  the  lack  of  an  ICD  which 
is  commonly  used  to  make  the  inlet  flow  uniform  [61].  Another  source  is  from  the 
unsteady  nature  of  the  IGV  wakes  such  as  vortex  shedding,  etc.  A  typical  BPF  tone 
waveform  is  shown  in  Figure  5.11  when  the  simulator  runs  at  a  constant  speed.  It  can 
be  found  from  the  figure  that  the  maximum  amplitude  of  BPF  is  about  9  times  the 
minimum  amplitude,  i.e.  the  maximum  variance  of  BPF  signal  is  about  19.1  dB.  In 
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Figure  5.10:  The  block  diagram  for  the  control  system  using  casing  microphones. 

the  experiments,  an  averaging  time  of  4  seconds  is  required  to  get  a  steady  value  for 
BPF  power.  The  long  averaging  time  of  4  seconds  would  lead  to  a  slow  response  in  the 
TEB  control  system  which  is  a  clear  drawback.  However,  tnis  control  approaches  is 
still  pursued  since  one  of  the  objectives  in  this  chapter  is  to  demonstrate  the  concept 
of  TEB  control  using  acoustic  signature. 

Another  note  on  Eqn.  (5.2)  is  that  the  step  size.  S,  needs  to  be  relatively  large, 
which  is  again  required  because  of  the  fluctuation  of  BPF  tone.  This  is  because  it 
would  be  difficult  to  detect  the  trivial  change  in  the  BPF  power  if  the  blowing  rate 
varies  insignificantly  as  a  result  of  the  use  of  a  small  step  size.  In  technical  words, 
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Time  (seconds) 

Figure  5.11:  A  typical  time  history  of  the  BPF  tone  in  casing  microphone. 

the  fluctuation  of  BPF  tone  may  lead  to  some  local  minima  as  well  on  the  curve  in 
Figure  5.8  in  addition  to  the  global  minimum,  which  may  prevent  the  controller  from 
converging  to  the  global  minimum  if  a  small  step  size  is  used.  In  the  experiments,  a 


step  size  of  0.2  volts  is  used. 
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5.4.2  Results  and  discussions 


This  section  describes  the  experimental  results  obtained  from  the  tests  conducted  in 
the  anechoic  chamber.  These  results  include  the  TEB  controller  performance,  and  the 
noise  reduction  in  the  far-field  and  inlet  microphones. 

Controller  performance 

As  mentioned  earlier,  the  purpose  of  the  optimal  feedback  control  system  is  to  achieve 
the  optimum  blowing  rate  by  searching  for  the  minimum  point  of  the  curve  shown  in 
Figure  5.8  according  to  the  control  algorithm  described. 

Figure  5.12  the  RMS  signal  time  history  of  the  BPF  tone,  i.e.  error  signal,  during  the 
adaptation  process  as  well  as  the  control  signal  driving  the  microvalves.  The  simulator 
is  running  at  the  speed  of  30,000  rpm.  The  air  source  pressure  for  TEB  is  80  psi  in 
the  plenum.  The  initial  conditions  of  the  MEMS  based  microvalves  are  fully  open.  It 
is  clear  that  the  control  system  is  effective  in  reducing  the  BPF  tone  in  the  casing 
microphones  to  the  minimum.  The  control  signal,  which  was  updated  every  4  seconds, 
is  able  to  converge  to  the  optimum  solution,  i.e.  the  optimum  blowing  rate  in  Figure  5.8. 
in  about  40  seconds.  It  is  interesting  to  note  that  the  control  signal  oscillates  around 
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the  optimum  steady-state  value  that  also  results  in  a  very  low  frequency  oscillation  of 
the  error  signal,  i.e.  the  control  system  meanders  around  the  optimum  blowing. 
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Figure  5.12:  The  convergence  process  of  the  optimal  controller. 


Acoustic  results 


Microphone  Signal  (Volts) 
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Similarly  as  in  Chapter  4,  the  sound  field  is  monitored  at  far-field  locations  in  the 
anechoic  chamber.  It  is  known  that  the  successful  wake-filling  leads  to  the  reduction  of 
the  tonal  noise  components  due  to  the  suppression  of  the  interaction  between  the  IGY 
wakes  and  the  fan  blades.  Thus,  the  amount  of  the  noise  reduction  at  the  tonal  noise 
can  serve  as  an  indication  of  the  successful  re-energizing  of  the  wakes.  Figure  5.13  shows 
the  noise  spectra  in  the  far-field  at  20°  of  the  simulator  axis  with  and  without  active 
TEB  control.  For  the  purpose  of  comparison  with  the  results  obtained  in  last  chapter, 
the  BPF  tone  occurs  at  9.000  Hz  and  corresponds  to  a  simulator  speed  of  30.000  rpm. 
These  spectra  demonstrate  a  reduction  of  7.9  dB  at  the  BPF  tone  and  reduction  of 
5.2  dB  at  the  first  harmonic  (2BPF).  These  results  agree  very  well  with  those  obtained 
using  Pitot  probes  as  sensors  presented  in  Chapter  4  with  reductions  of  8.2  and  6.1  dB 
at  the  BPF  and  2BPF  tones,  respectively.  Figure  5.14  shows  the  radiating  directivities 
at  the  BPF  tone  of  9,000  Hz  without  and  with  TEB  control  system.  As  observed,  the 
BPF  tone  is  reduced  over  the  10°  to  80°  sector.  The  sound  power  level  is  reduced  by 
4.9  dB.  Again,  the  reduction  is  very  close  to  the  4.4  dB  obtained  using  Pitot  probe  as 
error  sensors  in  last  chapter.  It  is  worth  to  point  out  that  the  difference  in  the  far-field 
acoustic  results  obtained  from  the  Pitot-probe-based  and  microphone-based  sensing 
approaches  might  be  from  the  intrusive  Pitot-probe  flow  sensors  themselves.  That 
is,  the  Pitot  probes  can  create  their  own  wakes.  These  acoustic  results  demonstrate 
the  effectiveness  and  potential  of  this  noil-intrusive  acoustic  sensing  strategy  for  active 
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Figure  5.14:  Directivity  plot  with  and  without  the  TEB  controller. 
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TEB. 


The  casing  microphone  output  with  and  without  active  TEB  is  also  recorded  after  the 
controller  has  converged.  The  results  in  frequency  domain  and  time  domain  for  the 
casing  microphone  behind  IGV  .4  are  shown  in  Figures  5.16a  and  5.16b,  respectively. 
As  seen  in  Figure  5.16a,  a  reduction  6  dB  is  achieved  at  the  BPF  tone.  This  result 
can  also  be  seen  in  Figure  5.16b  where  the  root-mean-square  (RMS)  signal  of  the  BPF 
tone  is  reduced  by  a  factor  of  2.  As  explained  above,  the  significant  variability  of  the 
power  of  the  tone  is  due  to  the  upstream  flow  distortion.  In  addition,  similarly  as  in 
the  far-field  acoustic  spectrum,  the  wake-filling  does  not  lead  to  the  reduction  in  the 
broadband  noise. 

Discussions  and  conclusions 

Based  on  the  comparison  of  the  results  using  Pitot  probes  and  the  microphone  sensors, 
it  demonstrates  the  feasibility  of  the  acoustic  sensing  approach  to  eliminate  the  need 
for  intrusive  flow  sensors  to  successfully  fill  the  wakes  and,  thus,  reduce  the  unsteady 
loading  experienced  by  the  blades.  However,  the  approach  has  a  main  drawback,  which 
is  the  slow  system  response  as  shown  in  Figure  5.12.  This  is  because  the  BPF  power 
signal  fluctuates  severely  due  to  the  unsteady  behavior  of  the  IGV-rotor  interaction. 
Thus,  the  controller  has  to  take  4  seconds  to  obtain  a  steady  time-averaged  value 
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of  the  BPF  tone  power.  It  is  important  to  remark  that  the  optimal  controller  had 
difficulties  converging  to  the  minimum  if  the  starting  point  was  not  relatively  close  to 
the  minimum.  This  slow  response  (  ~  40  seconds)  and  the  poor  reliability  makes  the 
control  system  not  practical  for  real  applications  where  the  system  has  to  track  sudden 
changes  in  engine  speeds.  This  limitation  led  to  the  development  of  the  second  control 
scheme  investigated  in  this  work. 


5.5  Acoustic  Sensing  Approach  Using  BPF  Magni¬ 
tude  and  Phase 


5.5.1  Control  algorithm 


The  slow  time  response  of  the  microphone-based  control  system  can  be  accomplished 
by  once  again  implementing  a  PID  controller.  To  this  end.  the  error  signal  needs  to 
be  able  to  provide  information  that  allows  to  distinguish  the  under  and  over  blowing 
conditions  through  a  change  in  the  sign  of  the  error  signal.  This  can  be  accomplished 
by  using  both  the  magnitude  and  phase  information  from  the  BPF  tone,  instead  of 
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Figure  5.17:  The  wake  profiles  in  undcrblowing  and  overblowing  cases. 


only  the  BPF  magnitude  information  used  in  the  first  control  scheme. 


Plotting  the  circumferential  flow  velocity  profile  on  the  face  of  the  fan  for  both  under 
and  over  blowing,  as  shown  in  Figm~  5.17,  clearly  indicates  that  the  two  conditions  are 
characterized  by  being  out-of-phase.  It  follows  that  the  phase  of  the  casing  microphone 
signals  will  also  show  a  180  phase  change  between  the  two  blowing  cases.  Thus,  the 
phase  change  information  in  the  casing  microphone  signals  in  conjunction  with  the 
BPF  tone  magnitude  can  be  used  by  a  PID  controller  as  error  signal. 


Preliminary  experimental  investigation  at  the  simulator  speed  of  30,000  rpm  verifies 
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the  above  idea.  The  rpm  signal  from  the  magnetic  pickup  of  the  simulator  serves  as  the 
phase  reference  for  measuring  the  casing  microphone  phases.  Figure  5.18  illustrates 
that  the  time-averaged  BPF  tone  phase  in  casing  microphone  2  changes  by  180  degree 
as  we  keep  increasing  the  TEB  air  source  pressure  (0  degree  is  defined  at  no  TEB). 
The  measurement  results  for  all  the  four  microphones  are  tabulated  in  Table  5.1  when 
the  simulator  runs  at  30,000  rpm.  It  indicates  that  the  BPF  tone  phase  of  all  the 
microphone  signals  change  by  about  180  degree  from  undcrblowing  to  overblowing 
case.  A  perfect  180°  change  is  not  achieved  because  the  BPF  tone  is  generated  not 
only  by  IGV-rotor  interaction  but  rotor-rotor,  rotor-EGV  and  inlet  flow  distortion  as 
well.  Note  that  the  phase  measurement  exhibits  a  large  variation  at  the  blowing  rate 
of  40  psi  and  50  psi  as  compared  to  the  other  blowing  rates.  This  is  probably  because 
the  actual  supply  air  pressure  has  some  variation  and  40  psi  or  50  psi  is  very  close  to 
the  optimum  blowing.  Therefore,  the  actual  pressure  at  the  measurement  40  psi  or  50 
psi  meanders  around  the  optimum  blowing,  i.e.,  alternating  between  overblowing  and 
underblowing  which  thus  leads  to  the  large  phase  variation.  Finally,  it  should  be  noted 
that  Table  5.1  shows  the  result  of  180°  phase  change  only  at  one  simulator  speed.  The 
results  for  other  speeds  will  be  addressed  in  the  latter  part  of  this  chapter. 


The  formation  of  the  error  signal  is  depicted  in  the  schematic  of  Figure  5.19.  The 
magnitude  of  the  BPF  tone  as  a  function  of  the  blowing  rate  is  shown  in  Figure  5.19a 
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Figure  5.18:  The  relationship  between  BPF  phase  and  TEB  air  source  pressure. 


Table  5.1:  BPF  tone  phase  as  a  function  of  blowing  rate. 


Mic.  \  TEB 

0  psi 

20  psi 

40  psi 

50  psi 

60  psi 

80  psi 

No.  1 

206° 

211° 

209°  -  224° 

44°  -  -13° 

3°  -  -36° 

-12° 

No.  2 

23° 

25° 

20°  -  -30° 

147°  -  177° 

162°  -  183° 

168° 

No.  3 

182° 

165° 

102°  -  133° 

20°  -  0° 

-2°  -  -19° 

2° 

No.  4 

-5° 

3° 

29°  -  81° 

130°  -  160° 

183°  -  191° 

180° 
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while  the  phase  is  plotted  in  Figure  5.19b.  It  is  important  to  note  that  the  phase  of 

a  signal  is  always  measured  relative  to  another  coherent  reference  signal.  The  RPM 

signal  from  the  magnetic  speed  pickup  sensor  was  used  as  the  reference  signal  for 

measuring  the  casing  microphone  phases.  The  error  signal  for  the  PID  controller  can 

be  constructed  by  first  assigning  a  value  from  -1  to  1  for  the  BPF  phase  of  0  to  jt. 

The  phase,  now  in  the  form  of  a  function  ranging  from  -1  to  1,  and  the  magnitude  of 

the  BPF  tone  arc  multiplied  to  yield  the  error  signal  in  terms  of  the  blowing  rate  as 

shown  in  Figure  5.19c.  Thus,  a  negative  error  signal  will  be  obtained  for  under  blowing 

) 

while  a  positive  error  signal  will  be  obtained  for  the  over  blowing  case.  In  addition, 
the  error  signal  will  vanish  at  the  optimum  blowing.  Hence,  a  typical  PID  algorithm 
can  be  used. 


In  order  to  implement  this  control  scheme,  a  real  time  phase  detector  is  designed.  A 
phase  detector  converts  the  relative  phase  between  two  signals,  the  microphone  and 
RPM  signals,  to  a  voltage  signal  so  that  it  is  then  input  into  the  digital  controller. 
A  phase  detector  with  the  characteristics  shown  in  Figure  5.20  is  proposed  for  our 
application  [62].  For  the  phase  detector  to  be  insensitive  to  the  amplitude  of  the  input 
signals,  the  BPF  tone  from  the  casing  microphones  and  the  RPM  signals  are  first 
converted  to  TTL  (transistor  to  transistor  logic)  signals.  This  is  important  since  the 
BPF  tone  from  the  turbofan  engine  is  fluctuating  severely  in  amplitude.  The  kernel 
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Figure  5.19:  Schematic  of  (a)  BPF  tone 


error  signal. 
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Figure  5.20:  The  characteristics  of  the  phase  detector. 

part  in  this  phase  detector  is  an  Exclusive-OR  gate  with  the  two  TTL  signals  as  inputs. 
It  is  known  that  the  output  of  an  Exclusive-OR  gate  goes  high  when  only  one  input 
goes  high,  and  the  gate  output  remains  at  a  logic  low  for  all  other  inputs.  Thus,  it  can 
be  shown  that  the  mean  value  of  the  gate  output,  which  can  be  obtained  by  a  low  pass 
filter,  exhibits  a  triangular  function  with  respect  to  the  phase  difference  between  the 
two  TTL  input  signals.  Appendix  B  has  the  details  for  the  implementation  of  a  phase 
detector. 

Both  the  magnitude  and  phase  information  for  the  BPF  tone  for  each  microphone  is 
input  into  the  digital  controller.  The  error  signal  for  the  PID  controller  is  then  obtained 
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ek  =  '%2bi(k)  * p(i)(k){-iy  (5.4) 

i= 1 

where  bt(k)  and  pi(k)  are  the  magnitude  (RMS  signal  from  the  adaptive  tone  extractor 
device)  and  the  relative  phase  (from  the  phase  detector  circuit)  of  the  ith  microphone 
(<  =  1,2, 3, 4  in  the  azimuth  order  behind  IGV  A,  B,  C .  D,  respectively  )  at  time 
tk.  respectively.  It  is  important  to  mention  that,  in  order  for  the  control  system  to 
operate  in  the  linear  region  of  the  curve  shown  in  Figure  5.20,  a  phase  shifter  should 
be  used  for  RPM  signal  so  that  the  relative  BPF  tone  phase  in  microphone  1  is  set 
to  zero  for  the  no-blowing  case.  The  term  (—1)*  in  Eqn.  (5.4)  is  included  to  account 
for  the  relative  phase  between  microphone  signals  due  to  the  acoustic  modes  generated 
by  the  IGV-rotor  interaction.  Since  the  4  IGYs  and  the  18  fan  blades  excites  the 
m  =  2,  —2,  6  and  10  circumferential  propagating  acoustic  modes  as  shown  in  Table 
3.6,  adjacent  microphones  are  out-of-phase.  Thus  \ without  the  term  (-1)'  in  Eqn. 

(5.4) ,  the  error  signal  will  vanish.  Note  that  the  error  signal  as  obtained  in  Eqn. 

(5.4)  serves  the  purpose  of  averaging  the  effect  of  the  four  wakes.  This  is  convenient  to 
account  for  variability  in  the  trailing  edge  blowing  between  the  IGVs.  Another  possible 
advantage  of  combining  the  error  signals  from  the  casing  microphone  as  in  Eqn.  (5.4) 
is  that  the  IGV-rotor  interaction  can  be  separated  from  the  rotor-rotor  and  rotor-EGY 
interactions  if  they  generate  different  m  order  circumferential  acoustic  modes.  This 


will  be  further  discussed  in  future  work. 
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Figure  5.21:  The  block  diagram  for  the  control  system  using  both  BPF  power  and 
phase  information. 

Figure  5.21  shows  the  block  diagram  of  the  control  system.  Compared  to  the  setup  in 
the  previous  scheme,  the  phase  detector  has  been  added  and  the  control  algorithm  is 
a  feedback  PID  system.  The  inlet  setup  with  four  microphones  flush  mounted  on  the 
inlet  casing  remains  unchanged. 


Controller  performance 


This  second  control  scheme  was  also  tested  in  the  anechoic  chamber  setup.  The  time 
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Figure  5.22:  The  time  histories  of  (a)  control,  (b)  BPF  tone  RMS,  and  (c)  BPF  tone 


phase  signals. 
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histories  of  the  control,  and  the  magnitude  and  phase  of  the  error  signals  are  recorded 
and  they  are  shown  in  Figure  5.22.  This  figure  shows  the  convergence  process  after 
the  system  is  turned  on  at  about  5  second.  It  can  be  found  that  the  control  signal 
converges  to  the  steady  state  of  5.5  Volts  in  about  4-5  seconds.  The  BPF  power  is 
significantly  reduced  in  about  4  seconds.  Completed  minimization  is  achieved  in  about 
10  seconds.  This  result  should  be  compared  to  about  40  seconds  shown  in  Figure  5.12. 
It  is  also  interesting  to  observe  the  real  time  error  phase  information  from  the  output 
of  the  phase  detector  circuit  shown  in  Figure  5.22c.  Before  the  controller  is  turned 
on,  the  phase  time  history  is  stable  with  minimum  variability.  After  the  controller  is 
turned  on,  as  the  error  signal  is  minimized  the  phase  goes,  in  an  average  sense,  through 
a  change  of  nearly  90°  as  expected  from  Figure  5.19b.  Moreover,  the  phase  also  shows 
a  significant  increase  in  the  variability.  This  can  be  explained  using  Figure  5.18  which 
show  an  abrupt  phase  change  near  the  optimum  blowing  rate.  This  indicates  that 
a  small  deviation  from  the  optimum  blowing  rate  can  cause  the  large  change  in  the 
phase.  In  the  real  system,  small  deviations  from  the  optimum  blowing  rate  can  always 
exist  even  though  the  control  system  has  converged.  As  a  matter  of  fact,  the  converged 
value  just  wanders  around  the  “true”  optimum  value.  In  addition  to  that,  inlet  flow 
distortion  due  to  the  lack  of  an  ICD  could  also  be  a  cause  since  it  is  an  unsteady  process. 
These  results  clearly  show  that  the  main  drawback,  i.e.  the  slow  time  response,  of  the 
first  control  scheme  was  overcome  in  this  new  implementation. 
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The  converged  value  of  the  control  signal  is  compared  to  that  obtained  using  the  first 
controller  developed  above.  Figure  5.12a  shows  the  steady  state  value  of  4.5  \olts  of 
the  control  signal  obtained  by  the  first  controller  while  the  converged  control  signal 
from  the  second  controller  is  about  5.5  Volts  in  Figure  5.22a.  The  1-volt  difference 
should  be  due  to  the  flow  condition  changes  as  a  result  of  slight  physical  difference  in 
mounting  the  experimental  setup  for  the  two  testing.  The  difference  may  come  from 
the  slight  changes  in  the  inlet  duct  axial  and  exhaust  duct  axial  direction.  Hence,  for 
comparison  purpose,  the  first  controller  is  also  tested  on  the  same  experimental  setup 
here.  The  result  for  the  converged  control  signal  of  about  5.5  Volts  is  illustrated  in 
Figure  5.23.  Therefore,  it  can  be  confirmed  that  the  second  controller  does  converge  to 
the  steady  state  value  as  the  first  controller  did,  leading  to  the  optimum  TEB  blowing 
rate  as  well.  The  control  signal  in  Figure  5.23  is  not  so  “noisy"  as  compared  to  that 
in  Figure  5.22a.  This  is  because  the  control  signal  from  the  first  controller  is  the  time- 
averaged  value  with  a  period  of  4  seconds.  It  can  be  expected  that  the  control  signal  in 
Figure  5.22a  would  be  similar  to  that  in  Figure  5.23  if  it  was  also  time-averaged  with 
the  same  period  of  4  seconds. 

Acoustic  Results 

The  acoustic  data  in  the  far-field  and  in  the  casing  microphones  are  also  recorded.  The 
figures  are  not  presented  here  since  they  are  very  similar  to  those  obtained  using  the 
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Time  (seconds) 

Figure  5.23:  The  steady-state  control  signal  from  the  controller  only  using  BPF  mag¬ 
nitude. 

first  control  scheme.  The  measured  far  field  pressure  at  20°  after  convergence  shows 
a  reduction  of  the  BPF  and  2BPF  tones  of  7.5  and  3  6  dB.  respectively.  In  addition, 
about  7.0  dB  noise  reduction  is  observed  at  BPF  u:.  all  the  four  casing  microphones. 
These  results  again  agree  very  well  with  the  attenuation  obtained  using  the  Pitot 
probes. 

Discussions  and  conclusions 

The  control  scheme  has  demonstrated  the  much  faster  convergence  rate  as  compared 
to  the  first  controller  due  to  the  use  of  the  PID  controller.  The  error  signal  for  the  PID 
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controller  can  be  formed  because  the  additional  information  of  BPF  phase  is  utilized 
besides  the  BPF  magnitude.  As  a  matter  of  fact,  signal  phase  in  many  engineering 
cases  contains  abundant  information  which  can  be  employed  in  one  way  or  another. 
This  is  another  excellent  application  case  of  making  use  of  phase  for  a  TEB  control 
system. 


In  addition  to  the  much  faster  convergence  rate,  the  control  system  can  always  converge 
to  the  global  optimum  in  the  cost  function  shown  in  Figure  5.8.  As  mentioned  earlier, 
the  BPF  tone  is  usually  so  fluctuating  that  it  can  lead  to  some  local  minima  in  the  cost, 
function.  The  worst  case  may  cause  the  first  controller  converge  to  a  local  minimum 
rather  than  the  global  minimum  if  the  time-average  period  and  step  size  in  Eqn.  (5.2) 
are  not  big  enough.  However,  the  local  minimum  in  the  cost  function  can  does  not 
affect  the  convergence  of  the  PID  controller  since  the  use  of  the  BPF  phase  allows  to 
distinguish  underblowing  and  overblowing  conditions. 


It  is  important  to  remark  that  this  control  approach  works  only  at  one  engine  speed. 
The  reason  is  because  the  phase  reference  signal,  which  is  the  rpm  signal  from  the 
magnetic  pickup,  varies  in  phase  with  the  simulator  speed.  It  follows  that  the  BPF 
tone  phase  relative  to  the  rpm  signal  cannot  be  set  to  zero  for  all  speeds  in  the  noblow¬ 
ing  case  once  it  is  initialized  to  zero  at  one  specific  speed  by  the  phase  shifter.  In 
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consequence,  the  phase  detector  cannot  work  in  the  linear  region  shown  in  Figure  5.20. 
One  possible  way  to  solve  this  problem  is  to  conduct  a  phase  survey  in  advance  for  all 
speeds  providing  the  phase  detector  with  a  prior  information.  Apparently,  this  method 
is  not  practical.  In  the  following,  a  new  scheme  is  proposed  and  some  preliminary 
experimental  testing  are  conducted  to  attempt  to  resolve  this  problem. 

The  novelty  of  the  proposed  control  scheme  is  to  use  the  casing  microphone  signal, 
instead  of  the  rpm  signal,  to  serve  as  the  phase  reference,  i.e.,  phase  difference  between 
adjacent  casing  microphones  as  depicted  in  Figure  5.24.  The  phase  difference  between 


Microphone  1 


Figure  5.24:  The  casing  microphones  mounted  on  the  inlet. 
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two  adjacent  microphones,  A0  =  0(f)  -0(0)  can  be  calculated  out  according  to  Eqn. 
(3.12)  to  be: 

A0  =  [m9  —  2nnNBFlt/Q0  +  (pm^]e=  \  ~  [rn&  ~  27mAsQf/60  +  o„4(]«=o  (5.5) 
=  7r  for  m  =  ±2,  ±6.  ±10  (5.6) 

Note  that  the  phase  difference  of  n  is  obtained  for  all  the  circumferential  modes  m  = 
±2,  ±6  and  ±10.  Equation  (5.6)  clearly  indicates  that  two  adjacent  microphone  signals 
are  out-of-phase.  It  is  important  to  note  that  Eqn.  (5.  ')  is  independent  of  the  simulator 
speed.  Ft.  It  is  this  conclusion  that  provides  the  possibility  to  develop  a  new  acoustic 
sensing  approach  which  can  work  for  a  broad  range  of  simulator  speeds. 

It  is  known  from  Figure  5.25a  (the  re-illustration  of  Figure  5.18).  the  BPF  tone  phase 
in  microphone  2  changes  from  180°  to  0°  as  the  TEB  rate  increase  from  underblowing  to 
overblowing.  Since  microphone  1  is  adjacent  to  microphone  2.  the  BPF  Tone  phase  in 
microphone  1,  as  illustrated  in  Figure  5.25b,  would  change  from  0°  to  180°  as  TEB  rate 
increases.  This  was  also  verified  by  the  measurement  results  for  microphone  1  as  found 
in  Table  5.1.  Figure  5.25c  is  the  differential  of  Figure  5.25a  and  Figure  5.25b.  It  is 
important  to  note  that  Figure  5.25c  is  independent  of  the  simulator  speed  Ft,  as  evident 
in  Eqn.  (5.6).  (The  phase  values  in  Figure  5.25a  and  Figure  5.25b  were  obtained  using 
magnetic  pickup  signal  as  phase  reference  which  would  change  with  speed.)  It  follows 
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Figure  5.25:  The  BPF  tone  phase  in  (a)  Microphone  1  referenced  by  rpm  signal,  (b) 
Microphone  2  reference  by  rpm  signal  and  (c)  Microphone  1  referenced  by  Microphone 
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that  Figure  5.25c  can  be  used  to  replace  Figure  5.19c  in  constructing  a  new  error  signal 
for  the  PID  controller  which  now,  however,  could  work  for  all  simulator  speeds. 

Preliminary  testing  are  performed  to  attempt  to  confirm  the  proposed  idea.  Tables 
5.2  and  5.3  list  the  10  measurements  on  the  BPF  tone  phase  of  microphone  1  with 
microphone  2  signal  as  the  reference  without  blowing  and  with  overblowing,  respec¬ 
tively.  The  overblowing  is  obtained  by  setting  the  blowing  air  pressure  to  80  psi.  The 
simulator  speeds  ranging  from  20,000  to  35,000  rpm  are  tested.  According  to  Eqn. 
(5.6)  or  Figure  5.25c,  all  the  measurements  are  supposed  to  be  around  180".  However, 
as  shown  in  Tables  5.2  and  5.3,  only  the  measurement  results  in  the  case  of  30,  000  rpm 
are  expected  since  they  have  a  mean  value  of  ft  =  166°  with  a  small  standard  deviation 
of  t  =  8  at  no  blowing,  and  ft  =  174  with  r  =  11  at  overblowing.  The  measure¬ 
ment  results  at  the  other  speeds  are  not  consistent.  The  phase  differences  between  the 
adjacent  microphones  are  not  out-of-phase  or  /  and  exhibiting  very  large  variances. 
The  measurements  on  other  two  microphones,  i.e.,  microphone  3  and  microphone  4 
demonstrate  similar  results  as  well. 

These  unexpected  results  are  mainly  due  to  the  flow  distortions  ingested  to  the  inlet 
because  of  the  lack  of  an  ICD.  When  interacted  with  the  rotor  blades,  the  flow  distortion 
can  also  generate  other  circumferential  acoustic  modes  in  the  inlet  in  addition  to  the 
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Without  Blowing 

RPM 

10  measurements 

F 

r 

20,000 

44 

45 

58 

65 

94 

111 

131 

122 

139 

147 

96 

38 

23,000 

36 

2 

14 

25 

2 

-73 

103 

-46 

-77 

-151 

-17 

68 

25,000 

29 

15 

-18 

18 

10 

12 

6 

14 

22 

4 

11 

12 

27,500 

132 

117 

117 

119 

131 

122 

125 

108 

136 

115 

122 

8 

30,000 

165 

166 

165 

147 

171 

1 

152 

137 

141 

155 

159 

166 

i . _  - 

8 

32,500 

-42 

40 

-4 

-15 

-43 

-81 

-49 

-52 

-51 

-51 

-35 

32 

35,000 

-16 

-35 

-5 

-35 

1.8 

-15 

0 

-32 

-35 

-33 

-20 

15 
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Table  5.3:  Phase  difference  between  two  adjacent  microphones  with  overblowing 
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modes  induced  by  IGY-rotor  interaction.  This  explains  the  deviations  of  the  mean 
values  from  180°  in  Tables  5.2  and  5.3.  In  addition,  the  flow  distortion  is  an  unsteady 
process  which  thus  results  in  the  large  variance  in  the  phase  measurements.  Another 
possible  source  for  the  unexpected  results  is  IGY-rotor  and  rotor-EG\  interactions. 
Although  they  are  cut-off  at  the  simulator  speed  tested,  they  can  contribute  to  the 
output  of  near-field  microphones  as  well.  In  Chapter  7,  acoustic  mode  separation 
techniques  will  be  proposed  to  filter  out  these  disturbance  modes  so  that  only  the 
modes  induced  by  IGV-rotor  interaction  are  extracted  out.  In  addition,  it  is  noted 
that  the  anticipated  measurements  obtained  at  30.000  rpm  in  Tables  5.2  and  5.3  is 
probably  because  the  IGY-rotor  interaction  happens  to  dominate  among  the  other 
noise  sources  at  this  particular  speed. 


5.6  Summary 


Two  active  TEB  control  systems  with  the  non-intrusive  acoustic  sensing  approach  were 
successfully  developed  and  tested  on  a  turbofan  simulator  in  the  anechoic  chamber.  The 
non-intrusive  casing  microphones,  which  are  compact  and  acceleration-compensated, 
proved  to  be  performing  very  well  in  this  application.  The  significance  of  using  the 
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non-intrusive  sensing  technique  in  the  control  system  lied  in  that  it  provided  a  more 
feasible  way  of  wake-filling  in  the  realistic  engine  environments. 


The  first  control  system  used  only  the  magnitude  of  the  BPF  tone  to  generate  the  error 
signal  for  the  optimal  controller.  The  system  did  converge  to  the  optimum  blowing, 
however,  with  a  slow  response  of  about  40  seconds.  The  slow  response  made  the  control 
system  not  practical  for  real  applications.  The  second  control  system  was  developed  to 
solve  this  limitation.  The  novelty  of  the  second  control  system  was  the  formation  of  the 
error  signal  for  a  PID  controller.  The  phase  signal  of  BPF  can  be  used  to  distinguish 
between  the  undcrblowing  or  overblowing  cases.  The  magnitude  information  of  BPF 
tone  was  used  as  a  weight  which  represents  how  far  away  the  current  blowing  rate  is 
from  the  optimum  blowing  rate.  In  doing  so,  a  PID  control  algorithm  could  be  used 
again.  The  convergence  rate  was  accelerated  bv  about  10  times  as  compared  to  the 
first  controller. 

A  novel  analog  adaptive  BPF  power  extractor  was  proposed  in  this  chapter  and  em¬ 
ployed  successfully  in  the  control  systems.  This  extractor,  the  key  technology  to  suc¬ 
cess,  had  the  capability  to  extract  the  BPF  tone  or  its  harmonics  adaptively  from  the 
microphone  output  when  the  simulator  speed  changed. 
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The  acoustic  results  obtained  from  both  systems,  including  far-ficld  directivity  and 
sound  power  reduction,  agreed  with  those  obtained  in  the  probe-based  experiments.  In 
the  first  control  system,  the  noise  reduction  of  7.9  dB  was  achieved  in  the  far  field,  which 
was  comparable  with  the  8.2  dB  reduction  with  Pitot-probe  based  control  system.  The 
sound  power  level  was  reduced  by  4.9  dB,  which  was  close  to  the  power  reduction  of 
4.4  dB  using  Pitot  probe  as  error  sensors.  Similar  results  were  also  obtained  in  the 
second  control  system.  These  results  demonstrated  the  effectiveness  and  potential  of 
this  noil-intrusive  acoustic  sensing  strategy  for  active  TEB. 


The  limitation  of  the  second  controller  was  that  it  could  only  work  at  one  engine 
speed.  This  was  because  the  phase  reference,  i.e.,  rpm  signal  from  the  magnetic  pickup 
changed  with  the  simulator  speed.  The  new  concept,  which  involved  using  adjacent 
microphone  signal  as  phase  reference,  was  attempted  to  solve  this  problem.  However, 
it  turned  out  that  further  work  is  needed  to  accomplish  the  goal. 


Chapter  6 


Analytical  Prediction  of  Inlet 
Sound  Field 


In  Chapter  5  and  Chapter  4,  experimental  studies  were  conducted  to  investigate  TEB 
control  using  a  small  scale  engine  simulator.  This  chapter  aims  to  conducting  analytical 
studies  on  the  active  TEB  using  an  existing  aeroacoustic  code.  Results  obtained  from 
the  studies  are  expected  to  provide  the  theoretical  guidance  for  the  development  of 
more  advanced  experimental  TEB  control  systems.  In  section  6.1.  the  objective  of  this 


Jinwei  Feng  Chapter  6.  Analytical  Prediction  of  Inlet  Sound  Field  157 

analytical  work  is  introduced.  A  TEB  model  is  proposed  in  section  C.2.  Using  the  TEB 
model,  the  sound  field  prediction  in  the  simulator  inlet  is  conducted  in  section  6.3  In- 
adapting  the  existing  aeroacoustic  CUP2D  code. 


6.1  Introduction 


In  Chapter  5.  a  microphone-based  TEB  control  system  was  presented  that  made  use 
of  the  BPF  tone  magnitude  and  phase  to  infer  the  wake-filling.  However,  the  system 
has  some  drawbacks  in  terms  of  the  response  time  or  the  effective  working  range  in 
terms  of  simulator  speed.  More  advanced  acoustic  sensing  approaches,  which  could 
involve  using  the  amplitude  and  phase  of  specific  acoustic  modes  contributing  to  the 
BPF  tone,  may  solve  this  problem.  To  this  end,  a  clear  understanding  of  inlet  sound 
field  with  respect  to  TEB  management  becomes  critical. 


Such  understanding  will  include  the  quantitative  relationship  of  inlet  sound  field  with 
respect  to  wake  management  approaches.  More  specifically,  the  objective  of  this  chap¬ 


ter  is  to  establish  a  toi.l  to  answer  questions  such  as: 
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-  How  many  acoustic  modes  in  the  inlet  are  significantly  excited  by  IGY-rotor  in¬ 
teraction?  Mode  separation  techniques  may  be  used  in  advanced  acoustic  sensing 
approaches  to  extract  out  these  acoustic  modes  from  other  noise  sources  such  as 
the  rotor-EGV  interaction.  This  may  be  required  since  a  TEB  control  system 
only  seeks  to  reduce  IGY-rotor  interaction,  thus  other  noise  sources  are  viewed 
as  “disturbances”  by  the  sensors  of  a  control  system. 

-  What  is  the  effect  of  the  jet-like  wake-filling,  commonly  encountered  in  practice, 
in  the  inlet  sound  field? 

-  What  is  the  effect  of  axial  spacing  between  IGYs  and  rotors  on  inlet  sound  field 
with  and  without  TEB? 


In  the  following  sections,  two  steps  are  taken  to  achieve  the  above  objectives.  The  first 
step  is  to  develop  an  analytical  TEB  model.  The  second  step  is  to  introduce  this  model 
into  an  existing  aeroacoustics  code  to  conduct  inlet  sound  prediction  studies  of  active 
TEB  approaches  in  turbofan  engines. 
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6.2  TEB  Model  Development 


To  the  author’s  knowledge,  this  is  the  first  research  on  the  modeling  of  TEB  model. 
Thus,  the  model  development  in  this  section  may  appear  crude  at  first.  However,  it  is 
the  key  first  step  to  provide  preliminary  results  and  serve  as  a  stepping  stone  towards 
more  in-depth  research  in  this  field. 


Before  the  development  of  the  TEB  model,  several  IGY  pure  wake  models  are  studied 
in  terms  of  their  mass  flow  deficit  and  momentum  flux  deficit  downstream  the  IGY.  As 
the  name  suggests,  mass  flow  rate  is  defined  to  be  the  time  rate  of  mass  flow,  i.e..  the 
product  of  flow  density,  velocity  and  cross-section  area.  Momentum  flux  is  defined  as 
the  time  rate  of  flow  momentum  per  cross-se  ,cion  area,  i.e..  mass  flow  rate  times  flow 
velocity  [55].  It  is  known  that  complete  wake-filling  implies  that  the  wake  momentum 
flux  deficit  and  wake  mass  flow  rate  deficit  are  eliminated  (i.e..  momentumless  wake). 
Therefore,  these  pure  wake  deficits  should  be  first  defined  so  that  the  required  amount 
of  blowing  air  can  be  determined  to  achieve  a  momentumless  wake. 


The  physical  features  of  the  time  averaged  non-turning  IGY  wakes  are  characterized 
by  measuring  the  semi-wake  width  <5,  and  the  wake  centerline  velocity  defect  B as 
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Figure  6.1:  The  parameters  characterizing  a  pure  wake  velocity  profile. 


shown  in  Figure  6.1.  Numerous  empirical  correlations  have  been  developed  for  these 
two  quantities  which  quantify  the  physical  nature  of  the  wakes.  Three  models  are 
investigated  here  developed  by  Majjigi  and  Glicbe  [63].  Lakshminarayana  and  Davino 
[64].  and  Silverstein.  Katzoff  and  Bullivant  [65],  respectively.  A  pure  wake  model  can 
be  expressed  as: 

W(y,s)  =  \Vdc(s)e-°m:il^2  (6.1) 

where  \V(y.  s )  denotes  the  flow  velocity  in  the  IGV  wake  which  is  a  function  of  pitchwise 
distance  y  and  streamwise  distance  from  the  IGV  trailing  edge,  s. 


The  three  models  are  now  described  in  terms  of  the  wake  parameters  5(s)  and  lFdc($). 
(a).  Majjigi-Gliebe  pure  wake  model: 
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For  the  Majjigi-Gliebe  model,  the  wake  parameters  are  defined  by 


SM 

C, 

0.2375  (s/C,)C°Dm  +  0.034125 
0.375(s/C/)Cp125 

(6.2) 

Wdc(s) 

Wo 

ro.2o°-3675(s/G/)  +  1.95 

D  7.65  (s/C,)  +  1.0 

(6.3) 

where  C/  is  the  IGY  chord:  CD  is  the  IGV  drag  coefficient;  IT'0  is  the  free  stream 
velocity. 


(b).  Lakshminarayana-Davino  pure  wake  model: 

For  the  Lakshminarayana-Davino  pure  wake  model,  the  wake  parameters  are  given  for 
two  regions.  For  the  near  wake  (0.06  <  s/C,  <  0.4),  they  are  given  by: 


(6.4) 


\Vdc{s)  _  0.168  C°p5  0.353 C°J 

W0  ~  (s/C,  +  0.12)0-22  +  (s/C,  +  0.12)0-22 

and  for  the  far  wake  ( s/C ,  >  0.4),  they  are 

6(s)  _  0.m(s/C,)C°j 

C,  -  s/C,  +  0.688 

Wdc(s)  __  0.9792 C°D5 

Ho  “  s/C,  +  0.688 


(c).  Silverstein-Katzoff-Bullivant  pure  wake  model 
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For  the  Silvcrstein-Katzoff-Bullivant  model,  the  wake  parameters  are: 


m 

C, 


^(«/C,  +  0.3) 


n-rfc(g) 

Wo 


0.68C/  \fC~D^sJC~i  +  0.15) 


(0.8) 

(6.9) 


It  is  important  to  note  that  the  three  pure  wake  models  above  are  developed  in  a 
2-dimensional  easeade  model  (i.e.,  y  and  s  directions  only).  The  normalized  mass 
flow  rate  deficit,  md,  and  normalized  momentum  flux  deficit,  Md,  are  computed  and 
used  to  compare  these  models.  The  normalized  mass  flow  rate  deficit  and  normalized 
momentum  flux  deficit  are  defined  as 


fOC 

fhd  = 

J  —  oo  1 

/>00  r 

Md  = 

L\ 

W0-  W(y.s) 
Wo 

W’n  —  W*2(y.  s) 


U’2 

u  0 


dy 

dy 


(6.10) 

(6.11) 


The  computed  mass  flow  rate  and  momentum  flux  deficits  are  presented  in  Figure 
6.2  and  6.3,  respectively.  These  results  indicate  that  only  the  Lakshminaravana- 
Davino  model  conserve  both  the  mass  flow  rate  and  momentum-flux  deficits  stream- 
wise  in  the  far  wake  downstream  the  IGY.  The  deficits  for  the  Majjigi-Glicbe  and 


Silverstein-Katzoff-Bullivant  models  are  decaying  streamwise.  Therefore,  the  far  wake 
Lakshminarayana-Davino  model  is  used  in  developing  the  active  TEB  model  in  the 
following  since  it  states  unequivocally  the  required  amount  of  blowing  air  to  achieve  a 
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momentumless  wake. 

A  model  for  the  TEB  is  now  proposed.  It  is  assumed  here  that  the  TEB  velocity  profile 
can  be  represented  as 

WT(y.  s)  =  -Wjc(s)  •  e"0'693  (6.12) 

where  the  parameters  ST{s)  and  define  the  jet  width  and  depth,  respectively. 

It  is  further  proposed  that  these  parameters  can  be  expressed  as  proportional  to  the 
pure  wake  width  5(s)  and  depth  B'dc(s).  That  is: 


*T(s) 

=  a0-<5(s) 

(6.13) 

=  b0-\]dc(s) 

(6.14) 

where  <S(s)  and  Wdc(a)  can  be  found  from  Eqn.  (6.6)  and  (6.7).  The  proportionality 
coefficients  a0  and  o0  are  the  width  and  depth  scaling  factors.  It  is  noted  that  the 
TEB  jet  width  and  depth  have  the  same  decaying  rate  along  ^-direction  as  those  of  the 
Lakshminarayana-Davino  far  wake.  Therefore,  as  in  this  Lakshminarayana-Davino  far 
wake  model,  mass  flow  rate  in  the  TEB  jet  model  is  also  conserved  downstream  the 
IGY.  This  conservation  of  the  model  is  important  since  the  mass  flow  rate  downstream 
should  be  constant  in  real  cases.  In  practice,  the  value  for  the  width  scaling  factor.  a0 
should  be  determined  by  the  TEB  physical  geometry  such  as  hole  diameter,  number, 
location  etc.  The  depth  scaling  factor.  b0  can  be  determined  by  the  given  TEB  mass 


Jimrei  Feng 


Chapter  6.  Analytical  Prediction  of  Inlet  Sound  Field 


1G4 


Figure  6.2:  Normalized  mass  flow  rate  deficit  of  the  pure  wake  models. 


Figure  6.3:  Normalized  momentum  flux  deficit  of  the  pure  wake  models. 
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flow  rate,  fnT ,  once  a0  is  known.  This  is  because  60  along  with  a0  determines  mT  as 


where  p  is  the  fluid  density.  Equation  (6.16)  is  obtained  by  substituting  Eqn.  (6.13) 
and  (6.14)  into  Eqn.  (6.12). 

The  mixed  flow  velocity,  i.e.,  wake  and  TEB  velocities,  TV™x(y,  .s)  is  then  obtained 
through  a  linear  summation  of  the  pure  wake  flow  and  TEB  jet  velocities  as 

WmiAy-s)  =  W(y,S)  +  \VT(y,s)  (6.17) 

Although  the  actual  jet  mixing  process  is  much  more  complicated  than  a  linear  sum¬ 
mation  [58],  Eqn.  (6.17)  provides  a  simple  model  for  analytical  studies  which,  however, 
does  not  violate  any  fluid  laws.  It  should  be  noted  that  the  linear  summation  in  Equa¬ 
tion  6.17  satisfies  mass  flow  rate  conservation  if  incompressible  flow  (usually  Mach 
number  less  than  0.65)  is  assumed.  This  condition  holds  in  this  research  since  the  inlet 
flow  Mach  number  in  the  tests  is  about  0.3.  This  is  important  since  the  aeroacoustics 
analysis  to  be  conducted  in  the  latter  part  of  this  chapter  will  focus  on  the  studies  of 
inlet  sound  pressure  with  relation  to  the  TEB  mass  flow  rate. 


P  f  WT{y,s)dy 

J -OC 

/OC  _  0.6931/2 

bQWdc(s)e~^dy 

-OC 


(6.15) 

(6.16) 
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A  comparison  is  made  with  reported  experimental  results  to  validate  this  model.  Ac¬ 
cording  to  the  experimental  work  by  Park,  etc  [22,  23],  the  mixed  flow  velocity  in 
the  momentumless  wake  case  was  found  to  be  “self-similar"  downstream  the  IGY.  i.e.. 
normalized  velocity  profiles  for  a  momentumless  wake  at  different  streamwise  distances 
collapsed  into  a  single  curve.  The  normalization  factors,  i.e.,  the  overshoot  H  and 
the  spread  L  which  usually  quantify  a  “jet-like”  wake-filling  profile,  are  described  in 
Figure  G.4  which  is  an  example  of  a  momentumless  wake  velocity  profile  at  0.75  IGY 
chord  downstream  with  ao  =  0.707  and  bo  =  1.414  used  in  Eqn.  (6.17).  The  pitehwise 
distance  is  normalized  by  the  IGY  thickness  while  the  flow  velocity  is  normalized  by 
the  free  stream  velocity.  For  comparison  purposes,  the  pure  wake  velocity  profile  of 
the  Lakshminaravana-Davino  model  is  also  shown  in  Figure  6.4.  Using  the  normaliza¬ 
tion  factors  H's  and  L’s  as  indicated  in  Figure  6.4,  obtained  from  the  momentumless 
wake  profiles  at  each  of  the  downstream  locations  of  s/C/  =  0.5,  0.75,  1.0  and  1.5.  the 
normalized  rnomentumless  wake  profiles  at  these  axial  locations  are  shown  in  Figure 
6.5.  As  seen,  the  mixed  flow  velocity  in  the  momentumless  wake  case  is  exactly  “self¬ 
similar”  as  expected.  In  addition,  the  decaying  rate  of  H  and  L  along  streamwise  is 
also  computed  to  be  (s/C/)-0-56  and  (s/C/)0-56,  respectively,  which  should  be  compared 
to  (s/C/)-0-92  and  (s/C/)03  reported  in  references  [22,  23].  These  results  show  that 
the  model  for  TEB  proposed  here  agrees  well  with  the  experimental  results  presented 
in  the  literature  [22,  23].  Thus,  this  comparison  provides  some  degree  of  confidence  in 
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the  model  to  conduct  sound  prediction  in  the  next  section. 


6.3  Sound  Field  Prediction  in  Simulator  Inlet 


6.3.1  CUP2D  with  active  TEB 


The  aeroacoustics  code  used  in  inlet  sound  pressure  prediction  is  the  code  referred  to  as 
CUP2D  developed  by  D.  Hanson  [48].  This  code  considers  a  typical  2D  cascade  model 
for  stator/rotor  interaction  as  sketched  in  Figure  6.6.  Either  blade  row  can  rotate  in 
either  direction.  As  the  free  stream  passes  through  the  first  row.  the  wake  is  generated 
which  imposes  on  the  second  row  blades  leading  to  the  unsteady  loading.  The  unsteady 
loading  on  the  blades,  when  coupled  with  the  fluid,  results  in  the  sound  generation. 
The  CUP2D  code  is  capable  of  predicting  inlet  sound  field  upstream  the  first,  row,  in 
between  two  rows,  and  downstream  row.  Compared  to  the  previously  existing  codes 
[66],  CUP2D  presents  a  more  complete  2D  cascade  analysis  in  which  flow  tangency 
boundary  conditions  are  satisfied  on  the  two  cascades  in  relative  motion  for  several 
harmonics  simultaneously.  The  new  features  of  the  code  include  unsteady  coupling, 


Flow  velocity  normalized  by  the  overshoot  H 
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Figure  6.4:  A  momentumless  wake  velocity  profile  (ao=0.707  and  6q = 1  - 4 1 4 ) . 


Figure  6.5:  Normalized  momentumless  wake  profile  (ao=0.707,  &o=l-414). 
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Figure  6.6:  The  2D  flat  plate  model. 


frequency  scattering,  and  flow  turning  between  first  and  second  row. 


The  physical  parameters  characterizing  the  2D  cascade  model,  which  are  required  to 
input  to  CUP2D,  are  now  described.  The  effective  fan  radius  R,  which  is  the  important 
parameter  in  converting  a  real  3D  stator-rotor  configuration  into  a  2D  cascade  model, 
is  assumed  to  be  85%  the  actual  fan  radius;  B,  is  the  blade  number  in  the  first  row 
with  a  rotating  speed  of  while  the  second  row  has  B2  blades  rotating  at  speed  Q2; 
the  blade  chords  of  the  first  and  second  row  are  defined  by  C \  and  C2,  respectively; 
y\  and  y2  are  the  translational  moving  direction  for  the  two  rows;  Si  and  s2  are  the 
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strcamwisc  coordinates  along  the  blades;  L  is  the  axial  flow  velocity;  1 !  and  1  ■>  are  the 
swirl  flow  velocities  while  Il'i  and  IT2  represent  the  velocities  along  s i  and  s-j  direction, 
respectively;  xs  is  the  axial  spacing  between  the  rows;  $i  and  On  are  the  blade  angles 
of  the  first  and  second  row;  and  G\  and  G<i  are  the  blade  gaps. 

It  is  important  to  note  that  the  code  uses  the  Silverstein-Katzoff-Bullivant  pure  wake 
model  for  the  first  row  [65],  In  addition,  the  code  utilizes  the  Fourier  integral  coefficients 
of  the  pure  wake  model  rather  than  the  model  equation  itself. 


The  existing  CUP2D  should  be  adapted  so  that  it  can  be  used  in  this  research. 
First,  the  Silverstein-Katzoff-Bullivant  pure  wake  model  should  be  replaced  by  the 
Lakshminarayana-Davino  model  since  the  latter  presents  a  more  reasonable  wake  pro¬ 
file  for  TEB  studies,  i.e.,  mass  flow  rate  conserved.  It  follows  that  the  Fourier  co¬ 
efficients  of  the  Lakshminarayana-Davino  model  should  be  the  derived.  The  second 
adaptation  is  that  the  TEB  model  should  be  introduced  into  the  CUP2D  code.  Sim¬ 
ilarly,  the  Fourier  coefficients  of  TEB  should  be  found.  Appendix  C  has  more  details 
regarding  the  code  adaptation.  As  a  brief  summary,  the  Fourier  coefficients  of  the 
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Lakshminarayana-Davino  pure  wake  model  arc  found  to  be: 

' inB\ [— $2sin[di  -  02)  +  xssin0\] 


ly/F  Wdc(Sl) 

F’l(si)  =  ~q~  ’  iri  p 


-(-)2 

V. 


exp 


cos9\ 


(G.18) 


where. 


Q  = 


\/FG\ cos9\ 


2.13  <5(s,) 

In  addition,  the  Fourier  coefficients  associated  with  the  TEB  jet  arc  obtained  as: 

'  inB\{—S2Sin(9y  —  +  orssin9]] 


(C.19) 


Cj(Sl)  =  ^-^eXp 


where. 


—(IlH)2 
QT  . 


exp 


cos9\ 


QT  = 


'■%( 


(6.20) 

y/HG^COsO] 

2.13a05(s,) 

(6.21) 

^0  *  W  dc  (  S 1 ) 

(6.22) 

These  Fourier  coefficients  in  Eqn.  (6.18)  and  Eqn.  (6.20)  need  to  be  input  to  the 
subroutine  GTWAKE  in  CUP2D. 


6.3.2  Prediction  results  in  simulator  inlet 


This  section  uses  the  adapted  CUP2D  code  to  analyze  the  simulator  running  at  30,000 
rpm  which  was  one  of  the  test  speeds  in  Chapters  4  and  5.  As  mentioned  earlier,  the 
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prediction  results  are  expected  to  provide  an  insight  on  the  quantitative  relationship 
between  the  inlet  sound  field  and  TEB  in  terms  of  mass  flow  rate  and  wake-filling  shape 
(i.e.,  jet-like  shape),  and  on  the  effect  of  the  axial  IGV-rotor  spacings.  The  acoustic 
information  obtained  would  be  helpful  in  developing  more  advanced  microphone-based 
controllers.  In  addition,  the  results  are  also  expected  to  provide  some  guidance  to  the 
design  of  TEB  configurations  like  the  number  of  holes,  hole  size,  and  so  forth. 

Figure  6.7  shows  the  schematic  diagram  of  the  2D  cascade  model  for  the  simulator.  As 
compared  to  Figure  6.6,  the  first  row  here  is  the  IGYs  and  the  second  row  is  the  rotor. 
In  addition,  the  swirl  flow  velocity  in  between  IGVs  and  rotor  is  zero  (i.e.,  \\  =  0 
and  6i  =  0  in  Figure  6.6  or  M\\  =  0  in  Figure  6.7)  since  the  non-turning  IGYs  are 
used;  Mxa .  Aa  and  pa  represent  the  flow  axial  Mach  number,  sound  speed  and  fluid 
density  in  the  region  upstream  of  the  IGV;  the  subscripts  a.  b  and  c  denote  the  region 
upstream  IGY,  in  between  IGV  and  rotor,  and  downstream  rotor,  respectively;  M\\ 
and  M\  ■,  are  the  swirl  Mach  numbers.  The  values  of  these  parameters  at  the  simulator 
speed  of  30.000  rpm  arc  listed  in  Table  6.1.  The  detailed  derivation  of  these  values 
can  be  found  in  Appendix  D.  Using  the  information  provided  by  the  table,  the  inlet 
sound  pressure  level  in  the  location  of  inlet  microphones  with  and  without  TEB  is  then 


predicted  using  CUP2D. 
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Tabic  6.1:  The  input  parameters  for  CUP2D  when  the  simulator  runs  at  30k  rpm. 


number  of  IGV  ( B\ )  ;  ^ 

number  of  rotor  blades  (B2)  '■  18 

Axial  Mach  number  upstream  IGV  ( Mxa )  :  0.3 

Axial  Mach  number  between  IGV  and  rotor  (MXb)  '■  0.3 

Axial  Mach  number  downstream  rotor  (A/xc)  :  0.6 

Swirl  Mach  number  between  IGV  and  rotor  {My )  :  0.0 

Swirl  Mach  number  downstream  rotor  {My2)  :  0.6 


Translational  Mach  number  of  the  1st  row  (Myi) 
Translational  Mach  number  of  the  2nd  row  {My2) 

IGV  gap/chord  ratio  (G\/C\) 
rotor  gap/chord  ratio  {G2/C2) 

normalized  a  nal  spacing  between  IGV  and  rotor  (A's / R) 

sound  speed  upstream  IGV  (A„) 

sound  speed  between  IGV  and  rotor  (Ab) 

sound  speed  downstream  rotor  (Ac) 

density  upstream  IGV  (pa) 

density  between  IGV  and  rotor  (pb) 

density  downstream  rotor  (pc) 

normalized  axial  location  for  SPL  printout  Xa/R 


0.0 

0.405 

1.31 

0.74 

2.395 

1127.6  feet/second 
1127.6  feet /second 
1141.1  feet/ second 
0.0702  Ibm/feet 3 
0.0702  Ibm/feet 3 
0.0674  Ibm/feet 3 
1.79 


drag  coefficient  Cp 


0.01 
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Figure  6.7:  Schematic  digram  of  the  2D  cascade  model  for  the  simulator. 

Baseline  results 

The  baseline  case  (i.c.  without  TEB)  is  first  calculated.  The  CUP2D  output  at  BPF 
and  2BPF  at  the  inlet  microphone  locations  is  shown  in  Table  6.2.  Note  that  the  code 
output  does  not  contain  any  information  of  mode  radial  order  since  radial  uniformity 
is  assumed  in  the  2D-cascade  model.  The  first  observation  from  the  table  is  that  BPF 
tone  dominates  the  inlet  sound  field.  The  BPF  tone  is  about  50  dB  louder  than  that 
at  2BPF  tone.  Furthermore,  the  circumferential  mode  m  =  2  dominates  the  BPF 
tone.  The  second  most  significant  mode  at  BPF  is  m  =  6  which  is  27  dB  quieter  as 
compared  to  that  of  m  —  2.  Therefore,  the  code  predicts  that  IGV-rotor  interaction 
primarily  generates  m  =  2  modes.  This  information  is  very  valuable  when  designing 
advanced  microphone-based  TEB  control  systems.  As  mentioned  before,  BPF  tone 
sensed  by  inlet  near-field  microphones  are  due  to  IGV-rotor,  rotor-rotor  and  rotor- 
EGV  interactions,  and  inlet  flow  distortion.  The  sound  prediction  result  here  suggests 
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Table  6.2:  The  CUP2D  output  at  BPF  and  2BPF  tones. 


BPF 

2BPF 

circumferential  mode 

SPL 

circumferential  mode 

SPL 

m=2 

132.2  dB 

m=8 

82.6  dB 

m=6 

105.2  dB 

m=4 

43.3  dB 

m=-10 

82.2  dB 

m=-14 

43.4  dB 

Total:  132.2  dB 

Total:  82.0  dB 

that  IGV-rotor  interaction  can  be  separated  from  other  noise  sources  by  extracting 
out  only  m  =  2  modes,  (assume  that  the  other  two  noise  sources  generate  different 
circumferential  modes).  Design  of  advanced  microphone- based  systems  using  mode 
separation  technique  will  be  discussed  in  the  next  chapter  as  <x  future  work. 


As  an  additional  note,  it  is  worth  to  point  out  that  the  absolute  SPL  (i.e.,  133.2  dB) 
from  the  prediction  is  comparable  to  that  obtained  from  the  experimental  simulator 
testing  (i.e.,  128  dB  in  Figure  5.15).  which  gives  some  degree  of  confidence  of  using 
CUP2D  in  the  analytical  studies  of  this  chapter.  However,  it  is  important  to  note  that, 
the  code  greatly  under-predicts  the  2BPF  tone  level. 


Jimvei  Feng 


Chapter  6.  Analytical  Prediction  of  Inlet  Sound  Field 


176 


Figure  6.8:  The  inlet  sound  pressure  level  with  respect  to  IGY-rotor  axial  spacing. 

The  effect  of  the  axial  spacing  Xs  on  inlet  sound  field  is  also  investigated  here  without 
TEB.  The  axial  spacing  is  one  of  the  major  design  specifications  for  modern  turbofan 
engines.  Therefore,  it  is  important  to  take  into  account  the  axial  spacing  effect  when 
designing  acoustic  sensing  approach  for  TEB  control.  Figure  6.8  depicts  the  results 
Vote  that  the  SPL  value  for  mode  m  =  2  is  used  to  represent  the  pressure  level  at  BPF 
tone  since  m  —  2  mode  dominates  all  other  modes.  The  y-axis  is  the  pressure  level  at 
BPF  in  the  inlet  microphone  locations:  x-axis  is  the  axial  spacing  between  IGV  and 
rotor  normalized  by  IGV  chord.  The  figure  shows  that  the  pressure  level  decreases  as 
the  axial  spacing  increases.  This  result  agrees  with  those  reported  in  the  literature 

[15]- 
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BPF  tone  reduction  versus  TEB  mass  flow  rate  and  wake-filling  shape 

The  wake-filling  shapes  studied  here  are  obtained  by  varying  «0  in  the  active  TEB 
model.  According  to  Eqn.  (6.13),  a  small  value  for  a0  would  lead  to  a  narrow  TEB 
jet  while  a  large  value  for  a0  means  a  wide  jet.  Figure  6.9  shows  several  wake-filling 
shapes  at  simulator  fan  face  (i.e.,  one  IGY  chord  downstream)  in  the  cases  of  o()  =  0.3. 
0.5,  0.707  and  0.9. 

Figure  6.10  illustrates  the  BPF  tone  reduction  resulted  from  these  wake-filling  shapes 
shown  in  Figure  6.9.  The  horizontal  axis  is  the  TEB  mass  flow  rate  normalized  by 
the  pure  wake  mass  flow  rate  deficit.  So  a  value  of  1.0  refers  to  the  amount  of  mass 
flow  rate  which  leads  to  the  moinentuinless  wake.  As  indicated,  the  maximum  noise 
reductions  (i.e.,  the  peaks  in  Figure  6.10)  do  not  occur  at  the  momentumless  wake  if 
the  wake-filling  shape  is  “jet-like".  The  closer  to  the  pure  wake  width  the  TEB  jet 
width  is,  the  closer  to  the  momentumless  wake  the  peak  approaches.  The  importance 
of  this  conclusion  lies  in  that  it  warns  that  the  TEB  hole  size  should  be  designed  such 
that  the  jet  width  should  be  as  close  as  possible  to  the  pure  wake  width.  In  doing  so, 
the  maximum  sound  pressure  reduction  can  serve  as  the  exact  acoustic  indicator  to 
a  momentumless  wake.  Such  design  would  also  leads  to  the  largest  amount  of  sound 
pressure  reduction  (see  the  peaks  in  Figure  6.10,  around  28  dB  for  the  case  of  uq  =  0.9 
while  only  15  dB  for  ao  =  0.3).  From  another  viewpoint,  the  conclusion  warns  that 
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more  mass  flow  rate  is  actually  needed  to  achieve  a  momcntuinlcss  wake  in  a  converged 
microphone-based  TEB  control  system  where  the  maximum  sound  pressure  reduction 
is  used  to  infer  the  jet-like  wake-filling. 

The  effect  of  axial  spacing  on  TEB 

The  effect  of  the  IGY-rotor  axial  spacing  on  active  TEB  in  terms  of  the  inlet  sound 
pressure  is  also  examined.  Three  cases  of  Xs  =  0.5,  0.75  and  1.0  IGV  chord  are  studied. 
Figure  6.11  shows  the  results  when  the  jet  width  is  that  of  oo  =  0.5.  It  indicates  that 
sound  pressure  reduction  increases  with  Xs,  given  the  same  amount  of  TEB  mass  flow 
rate(  when  <  0.6).  For  instance,  if  a  limited  amount  of  TEB  mass  flow  rate  is  available, 
e.g.  0.5.  a  reduction  of  15  dB  is  achieved  in  the  case  of  A's  =  l.O'lGY  chord  while  only 
about  8  dB  is  obtained  in  Xs  —  0.5  IGY  chord.  The  second  observation  from  the  figure 
is  that  the  peak  gets  closer  to  the  solid  vertical  line  (i.e.,  the  momentumless  wake)  as 
the  axial  spacing  decreases.  It  implies  that  if  the  maximum  sound  reduction  in  the 
inlet  is  employed  as  an  indicator  to  infer  a  momentumless  wake  as  did  in  the  optimal 
controller  in  Chapter  5,  it  would  work  better  in  the  smaller  spacing  Xs. 
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Figure  6.9:  The  four  mon.enturaless  wake-filling  shapes. 
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Figure  6.10:  BPF  tone  reduction  in  the  four  wake-filling  cases. 
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Figure  6.11:  Noise  Reduction  with  TEB  mass  flow  rate. 


6.4  Summary 


This  chapter  conducted  analytical  studies  on  the  active  TEB  and  sound  prediction 
using  CUP2D.  In  the  first  part  of  the  chapter,  an  analytical  model  for  TEB  was  pre¬ 
sented  which  was  based  on  Lakshminaravana-Davino  pure  wake  model.  The  active 
TEB  model  was  mass  flow  rate  conserved.  In  addition,  the  momentumless  wake  pro¬ 
files  were  self-similar  downstream  of  the  IGV.  This  result  agreed  with  the  experimental 
results  obtained  from  Park[22,  23].  The  second  part  of  the  chapter  was  devoted  to 
the  analytical  predictions  of  inlet  sound  pressure  using  an  adapted  CUP2D  code.  The 
adaptation  of  the  CUP2D  code  consisted  of  two  steps.  First,  FFT  coefficients  for  the 
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Lakshminarayana-Davino  pure  wake  model  were  computed  to  replace  the  Sil verst ein 
model  which  was  originally  used  in  CUP2D.  Second,  the  active  TEB  was  introduced 
into  the  CUP2D  code. 

The  simulator  at  a  rotating  speed  of  30,000  rpin  was  analyzed.  The  effects  of  "jet¬ 
like”  wake  profiles  on  inlet  sound  field  and  IGV-rotor  axial  spacings  on  active  TEB 
were  examined.  It  was  found  that  maximum  sound  pressure  reductions  do  not  occur 
at  the  momentumless  wakes  if  the  wake  profiles  are  “jet-like”.  More  blowing  air  than 
the  amount  which  caused  maximum  sound  reduction  was  needed  to  achieve  a  momen¬ 
tumless  wake.  The  analysis  also  warned  that  TEB  holes  should  be  such  designed  that 
the  jet  width  approaches  the  pure  wake  width  as  possible.  The  investigation  on  the 
IGV-rotor  axial  spacing  effects  on  active  TEB  and  inlet  sound  field  demonstrated  that 
a  larger  spacing  would  lead  to  a  larger  inlet  sound  reduction  using  a  limited  amount 


of  blowing  air. 


Chapter  7 


Conclusions  and  Future  Work 


The  TEB  control  systems  using  different  sensing  approaches  are  summarized  and  main 
conclusions  of  the  research  are  drawn  in  this  chapter.  Future  work  is  also  discussed. 
It  includes  improvements  of  the  acoustic  sensing  approaches  and  development  of  other 
practical  sensing  techniques. 
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7.1  Conclusions 


To  the  author’s  knowledge,  this  was  the  first  research  effort  to  conduct  studies  on 
the  feasibility  of  employing  active  TEB  control  to  reduce  stator-rotor  interaction  and 
the  associated  HCF  damage  in  blades.  Active  flow  control  systems  with  different 
wake  sensing  approaches  had  been  successfully  implemented  in  the  simulated  engine 
environment. 

The  Pitot-probe  based  flow  control  system  was  developed  and  implemented  as  a  first 
step  toward  more  in  depth  investigations  of  active  TEB.  A  PID  controller  was  used  to 
generate  control  signal  for  MEMS-based  microvalves  to  adjust  the  TEB  rate.  Results 
had  shown  that  the  control  system  produced  a  very  uniform  flow  field  downstream 
of  the  IGV  in  a  rotating  environment.  The  ability  of  the  system  to  achieve  effective 
wake  filling  when  subjected  to  a  change  in  inlet  flow  conditions  demonstrated  the 
feasibility  and  advantage  of  active  flow  control.  It  was  also  demonstrated  that  far-field 
sound  pressure  levels  at  BPF  tone  were  reduced  by  TEB.  The  maximum  tone  reduction 
obtained  was  8.2  dB  at  30,000  rpm  and  7.3  dB  at  40,000  rpm.  Sound  power  level  at  the 
BPF  tone  was  reduced  by  4.4  dB  at.  30.000  rpm  and  2.9  dB  at  40.000  rpm.  respectively. 
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Another  significant  contribution  of  this  dissertation  lied  in  the  microphone-based  con¬ 
troller  which  provided  a  potential  practical  TEB  approach  for  the  implementation  in 
realistic  engines.  Microphones  were  flush  mounted  near  the  rotor  behind  the  IG\  s  to 
sense  the  tonal  noise  at  BPF.  In  the  first  control  scheme,  optimum  control  algorithm 
was  used  to  find  the  optimum  blowing  rate  by  searching  for  the  minimum  of  the  BPF 
magnitude.  However,  this  approach  had  several  drawbacks  with  the  main  one  to  be  the 
slow  time  response  of  the  system,  i.e.,  about  40  seconds  to  converge  to  the  optimum. 
The  novelty  of  the  second  control  scheme  was  using  both  the  BPF  tone  magnitude  and 
phase  to  form  the  error  signal  for  a  PID  controller.  The  phase  information  of  the  BPF 
tone  was  used  to  distinguish  between  the  underblowing  or  overblowing  cases.  The  con¬ 
vergence  rate  was  accelerated  by  about  another  order  of  magnitude  as  compared  to  the 
first  controller.  The  acoustic  results  obtained  from  both  schemes,  including  far-field 
directivity  and  sound  power  reduction,  agreed  well  with  those  obtained  in  the  Pitot- 
probe  based  experiments.  The  maximum  noise  reduction  of  7.9  dB  was  achieved  in  the 
far  field  and  the  sound  power  level  was  reduced  by  4.9  dB.  These  results  demonstrated 
the  effectiveness  and  potential  of  this  non-intrusive  acoustic  sensing  strategy  for  active 
TEB.  It  is  also  worth  to  mention  that  a  novel  analog  adaptive  BPF  power  extractor 
was  proposed  and  employed  successfully  in  the  control  systems. 


Preliminary  analytical  studies  were  also  conducted  on  the  modeling  of  TEB.  An  ana- 
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lvtical  model  for  TEB  was  first  proposed  based  on  Lakshminarayana-Davino  pure  wake 
model.  The  momentumless  wake  profiles  computed  from  this  mass  flow  rate  conserved 
model  were  self-similar  downstream  of  the  IGY .  The  TEB  model  was  then  introduced 
into  a  2D  aeroacoustic  code  to  investigate  effect  of  TEB  in  the  tonal  sound  generation. 
The  effects  of  double-hump  wake  profiles  on  the  inlet  sound  field  were  then  examined. 
A  key  result  showed  that  more  blowing  air  than  the  amount  leading  to  maximum 
noise  reduction  was  actually  needed  to  achieve  a  momcntumless  wake  in  a  converged 
microphone-based  control  system,  if  the  filled  wake  profile  had  a  double-hump  shape. 

7.2  Future  Work 


Since  active  TEB  research  is  still  in  its  first  stage,  it  is  far  from  being  well-developed 
for  practical  implementations.  It  needs  further  investigations  and  analyses.  In  this 
section,  a  potential  improvement  on  the  TEB  controller  with  acoustic  sensing  approach 
will  be  first  discussed.  Then,  the  future  work  including  other  sensing  approaches  will 
be  addressed. 

As  the  first  recommendation,  an  inlet  flow  control  device  (ICD)  should  be  used  in  the 
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experiments.  The  purpose  of  the  ICD  is  to  minimize  the  spurious  effects  of  ground 
testing  on  acoustic  measurements  by  breaking  up  incoming  vortices  and  eliminating 
inlet  flow  distortions.  As  stated  before,  the  lack  of  the  ICD  has  the  adverse  effects 
on  the  acoustic  sensing  approaches  developed  in  Chapter  5.  These  effects  include  the 
severe  fluctuation  in  both  BPF  tone  magnitude  and  phase.  Use  of  the  ICD  would 
bring  an  immediate  benefit  to  the  TEB  control  systems  developed  in  Chapter  5.  i.e.. 
the  system  response  of  the  first  controller  can  speed  up  by  reducing  the  average  time 
while  the  BPF  tone  phase  signal  in  the  second  controller  would  exhibit  less  variability. 

The  second  recommendation  involves  using  a  mode  extraction  technique  for  advanced 
acoustic  sensing  approaches  to  separate  IGV-rotor  interaction  from  other  noise  sources. 
This  mode  separation  technique  would  help  improve  the  second  controller  in  Chapter 
5  so  that  it  could  work  for  all  simulator  speed.  As  described  at  the  end  of  Chapter  5,  a 
new  control  algorithm  was  attempted  to  improve  the  second  controller.  However,  the 
preliminary  testing  results  were  not  successful.  This  is  because  the  sound  field  in  the 
inlet  is  generated  not  only  by  IGV-rotor  interaction  but  by  other  noise  sources  such  as 
rotor-rotor,  rotor-EGV  interactions,  and  inlet  flow  distortion. 

The  basic  idea  of  mode  separation  technique  is  to  perform  wavenumber  transform  using 
an  array  of  casing  microphones  to  extract  out  the  mode  amplitude  and  phase  of  the 
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significant  acoustic  modes  induced  by  the  IGY-rotor  interaction.  These  significant 
acoustic  inodes  induced  by  IGV-rotor  interaction  can  be  obtained  from  the  analytical 
prediction  of  the  aeroacoustic  code  CUP2D.  An  array  of  casing  microphones  provides 
the  spatial  information  of  the  inlet  sound  field  while  the  wavenumber  transform  can 
be  viewed  as  a  technique  to  design  “a  spatial-domain  digital  filter"  as  compared  to  the 
regular  Fourier  transform  to  design  “a  time-domain  digital  filter” .  In  the  following,  use 
of  a  circular  microphone  array  is  briefly  described.  The  purpose  here  is  to  introduce  the 
mode  separation  concept  instead  of  the  implementation  details,  i.e.,  a  stepping  stone 
toward  further  investigation. 


In  order  to  extract  out  the  circumferential  mode  rn.  a  circular  array  of  Nc  microphones  is 
placed  circumferentially  on  the  inlet  case  as  shown  in  Figure  7.1.  If  the  nth  microphone 
signal  is  denoted  as  r(n)  (n  =  0, 1, 2,  •  •  •  ,  Nc),  the  discrete  wavenumber  transform  in 
the  azimuth  direction  can  be  expressed  as 

A  4-1 

e(Ae>)  —  ^  r(n)e~JnA<t>  (7.1) 

n—0 

where  j  =  y/— T;  A<£,  which  is  equal  to  is  the  phase  delay  applied  to  the  micro¬ 
phones  to  extract  the  circumferential  mode  m;  The  complex  wavenumber  coefficient 
e(A<£),  which  is  referred  as  the  error  signal  for  control  systems,  can  be  written  in  the 
form  of  e(A(p)=|e(A<£)|Ze(A<£),  where  |e(A^)|  is  the  mode  amplitude  and  Ze(A^)  is 
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Figure  7.1:  Mode  magnitude  and  phase  extraction  using  a  circular  microphone  array. 

the  mode  phase.  Note  that  multiple  circumferential  modes  can  be  extracted  out  by  us¬ 
ing  different.  A0  in  Eqn.  (7.1).  Once  the  mode  magnitude  and  phase  of  the  significant 
acoustic  modes  are  obtained,  they  can  be  carried  over  to  Figure  5.19  to  substitute  the 
BPF  tone  magnitude  and  phase,  respectively  to  implement  the  PID  control  algorithm. 

In  addition  to  the  acoustic  sensing  approaches,  other  practical  sensing  techniques  can 
be  developed  as  well.  These  techniques  may  consist  of  mounting  eddy  current  sensors 
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or  light  probes  in  the  casing  radially  outside  a  blade  row  [G7]  to  directly  measure  the 
blade  vibrations.  Eddy  current  sensors  work  in  the  principle  of  detecting  magnetic 
flux  changes  caused  by  the  spinning,  blades  while  light  probes  involve  using  a  laser 
beam  to  detect  blade  vibration  [68.  69].  Since  the  relatively  small  blade  vibration  is 
embedded  in  the  large  blade  rigid-body  motion  as  the  rotor  passes  by  the  sensors, 
the  most  challenging  part  of  the  research  may  be  the  extraction  of  the  rotor  blade 
vibration  information  in  the  sensor  output  signal.  In  conjunction  with  these  sensors, 
the  appropriate  TEB  control  algorithms  should  be  also  developed. 

Other  future  work  may  include  using  a  3D  aeroacoustic  code  to  investigate  TEB  span- 
wise  optimization  to  minimize  the  amount  of  blowing  air.  The  use  of  the  minimum 
amount  of  TEB  air  is  very  important  since  in  real  engines  it  needs  to  be  bled  from 
the  compressor.  Excessive  use  of  bleeding  air  will  impair  the  compressor  efficiency. 
Minimization  of  TEB  air  may  be  obtained  by  blowing  less  air  to  the  blade  hub  while 
more  air  to  the  blade  tip  since  the  blade  tip  exhibits  relatively  large  vibrations.  In 
addition,  TEB  control  to  re-energizc  the  turning-IGY  wakes  should  be  also  pursued 
since  turning  IGVs  or  stators  are  commonly  used  in  real  aircraft.  The  wake  profile  of 
the  turning  IGVs  is  not  symmetrical  and  the  axis  of  the  wake  profile  may  change  with 


the  inlet  flow  conditions. 
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Appendix  A 


Adaptive  BPF  Tone  Extractor 


This  appendix  describes  the  two  major  components  used  in  the  adaptive  BPF  >,one 
extractor  in  Chapter  5.  In  addition,  some  implementation  issues  of  the  extractor  are 
discussed. 

The  proposed  adaptive  BPF  tone  extractor  consists  of  two  main  electronic  components 
which  are  switch-capacitor  filter  (SCF)  and  phase-locked  loop  (PLL)  as  shown  in  Figure 
5.6.  The  SCF  functions  as  a  bandpass  filter  while  the  PLL  provides  adaptability  in 
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the  system.  Besides  functioning  as  a  bandpass  filter,  an  SCF  can  also  be  configured 
as  a  notch,  low-pass  or  high-pass  filter.  They  have  been  commercially  available  in  the 
market  in  the  recent  years.  The  cost  for  each  chip  is  around  a  couple  of  dollars  which 
is  much  cheaper  than  a  DSP  which  is  required  in  the  implementation  of  conventional 
adaptive  filters.  SCFs  need  no  external  capacitors  or  inductors,  which  allows  consistent, 
repeatable  filter  designs.  Furthermore,  SCFs  can  have  little  sensitivity  to  temperature 
changes.  This  is  because  SCFs  are  clocked,  sampled-data  systems;  the  input  is  sampled 
at  a  high  rate  and  is  processed  on  a  discrete-time,  rather  than  continuous,  basis.  This 
is  a  fundamental  difference  between  SCFs  and  conventional  active  and  passive  filters, 
which  are  also  referred  to  as  “continuous  time”  filters.  The  chips  from  Liner  Technology 
Corporation  implement  the  bandpass  filters  with  the  transfer  functions  in  the  s-domain 
as  (take  2nd  order  as  examples)  [70]: 


P bandpass  —  HqBP  ^ 


/  Q 


S2  4*  s{iOn/ Q)  +  UJ 2 


Wn 


a.’c/A/ 


(A.l) 
(A. 2) 


where  Q  is  the  quality  factor.  H0bp  is  the  gain  at  the  center  frequency  u/„;  uc  is  the 
clock  frequency  and  M  normally  is  equal  to  100  or  20  which  depends  on  the  specific 
chips.  It  is  noted  from  Equation  A.l  that  the  phase  shift  is  zero  at  the  center  frequency. 


One  the  other  hand,  the  PLL  is  used  to  increase  the  input  signal  frequency  by  M  times 
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in  this  application  to  provide  the  clock  signal  for  SCF.  Phase  locked  loops  are  actually 
widely  used  in  wireless  communication.  Basically,  a  PLL  is  a  feedback  system  which 
synchronizes  an  oscillator  in  phase  and  frequency  to  an  incoming  signal.  As  can  be 
seen  in  Figure  A.l,  a  PLL  has  three  main  parts  of  a  phase  comparator,  a  low-pass 
filter  and  a  voltage  controlled  oscillator  (VCO).  The  phase  comparator  measures  the 
phase  difference  between  the  input  and  output  signals  and  produces  an  error  signal 
proportional  to  the  measured  phase  difference.  If  the  low-pass  filter  is  disregarded 
for  the  moment,  the  error  signal  will  drive  the  VCO,  changing  its  frequency  so  as  to 
minimize  the  phase  difference  between  the  input  and  output  signals.  Ultimately,  the 
frequency  of  the  output  signal  will  be  identical  to  that,  of  the  input  signal  and  will 
follow  every  change  of  the  latter.  Since  the  output  signal  is  obtained  by  a  frequency 
divider  with  a  dividing  factor  of  M,  the  frequency  of  the  clock  signal  is  M  times  of 
the  input  signal.  To  enable  a  PLL  acquire  lock  in  the  first  place,  it  should  of  course 
be  designed  so  that  the  difference  between  the  input  and  output  frequencies  is  not  so 
large  that,  after  filtering,  there  is  no  control  signal  for  the  VCO. 


The  static  performance  of  the  extractor  mainly  depends  on  the  SCF  while  the  dynamic 
performance  is  primarily  determined  by  the  performance  of  PLL.  The  performance  of 
SCFs  varies  from  chips  to  chips.  The  most  important  parameter  of  our  interests  is  the 
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Figure  A.l:  The  block  diagram  for  a  basic  phase-locked  loop. 

maximum  center  frequency  error,  Ce,  which  is  defined  as: 

Ce  =  (A.3) 

/o 

where  f'0  is  the  measured  frequency  in  the  output  while  /0  is  the  ideal  frequency  of 
output  signal.  The  SCF  used  here  in  the  application  is  LTC1264  manufactured  by 
Linear  Technology  Incorporation.  The  center  frequency  error  of  this  chip  is  0.2%. 
This  can  be  considered  as  a  tiny  error.  The  center  frequency  can  be  up  to  230kHz. 
The  performance  of  PLL  is  mainly  governed  by  the  loop  components,  i.e.,  the  phase 
comparator,  the  loop  filter,  the  voltage  controlled  oscillator  used  in  the  loop.  Especially 
the  loop  filter  plays  the  most  important  role  among  them.  An  active  filter  can  be  used  as 
the  loop  filter  to  obtain  a  better  loop  performance.  The  PLL  used  here  is  74HCT4046N 
from  National  Semiconductor  Incorporation.  It  has  three  phase  comparators  which  are 
exclusive-OR,  edge-trigger  JK  flip-flop  and  edge  trigger  RS  flip-flop.  In  addition,  it  has 


Jimvei  Feng 


Appendix  A.  Adaptive  BPF  Tone  Extractor 


200 


a  built-in  YCO  with  an  excellent  frequency  linearity. 


The  adaptive  BPF  power  extractor  is  the  key  for  the  control  system  to  success.  As 
a  matter  of  fact,  there  are  several  important  implementation  issues  which  need  to  bo 
addressed  here.  The  original  rpm  signal  from  the  built-in  magnetic  pickup  should  be 
pre-conditioned  before  it  can  be  used  by  the  adaptive  BPF  extractor.  There  are  several 
reasons  for  it.  First,  the  original  rpm  signal  is  noise-contaminated.  Especially,  it  is 
strongly  disturbed  by  the  line  frequency,  i.e.,  50  Hz  disturbance.  Therefore,  a  high 
pass  filter  with  the  cut-on  frequency  of  100  Hz  is  used.  Secondly,  the  amplitude  of 
rpm  signal  is  at  a  very  low  voltage  level  and  varies  with  simulator  speed  as  well.  In 
the  experiments,  the  rpm  signal  should  be  converted  to  a  TTL  (transistor  to  transistor 
logic)  signal.  Lastly,  the  rpm  signal  should  be  divided  by  six.  Use  of  the  undivided 
signal  in  frequency  would  lead  the  extractor  to  work  onlv  in  several  discrete  narrow 
frequency  ranges.  This  is  because  the,  six  magnetic  sensors  are  not  mounted  equally  in 
space  in  the  duct  circumference.  In  consequence,  the  pulses  in  the  rpm  signal  are  not 
equally  spaced  in  time  leading  to  the  broadband  frequency  content  in  the  input  signal. 
Experiments  show  that  the  pulses  in  the  rpm  signal  is  strictly  equally  spaced  in  time 
after  the  rpm  signal  is  frequency-divided  by  six  and  then  frequency-multiplied  by  six 
using  PLL.  Another  issue  is  that  enough  gain  should  be  provided  to  the  microphone 
output  signals  before  they  are  fed  to  SCFs.  This  helps  to  introduce  less  broadband 
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electronic  noise  to  the  circuits.  As  we  docs  to  a  traditional  ADC.  the  input  signal  range 
should  be  amplified  or  attenuated  so  that  it  matches  the  ADC  range  in  order  to  make 
full  use  of  the  ADC  resolution.  The  same  rule  applies  to  an  SCF  because  an  SCF.  like 
an  ADC,  also  works  on  the  discrete  signals. 

Lastly,  it  is  also  worth  to  mention  that  the  adaptive  BPF  tone  extractor  can  find 
applications  in  the  areas  of  active  noise  and  vibration  control,  condition  monitoring 
and  fault  detection  system  for  rotatory  machinery,  and  communication  system,  etc. 


Appendix  B 


Phase- Voltage  Converters 


This  appendix  describes  the  working  principles  of  the  phase-voltage  converter  (PYC) 
used  in  Chapter  5.  The  block  diagram  of  the  PVC  implementation  is  presented  and 
the  use  of  the  kernel  component  of  XOR  is  also  introduced. 


Phase-voltage  converter  (PVC)  is  an  important  component  in  the  microphone-based 
TEB  control  system.  The  phase  signal  of  BPF  tone  needs  to  be  converted  to  a  voltage 
signal  so  that  it  can  be  used  by  the  digital  signal  processor.  The  basic  requirement  for 
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this  conversion  is  that  the  controller  should  be  able  to  determine  whether  it  is  under 
or  over  blowing  by  detecting  the  voltage  signal. 


Without  loss  of  generality,  negative  voltages  are  specified  to  represent  underblowing 
while  positive  voltages  to  represent  overblowing.  Therefore,  in  Table  5.1.  90°  through 
180°  in  microphone  #  1  and  #  3  should  be  converted  to  negative  voltages  since  90" 
through  180°  are  known  for  underblowing.  Similarly,  0°  through  90°  in  microphone  # 
1  and  #  3  should  be  converted  to  positive  voltages  because  they  occur  in  overblowing 
case.  On  the  other  hand,  0°  through  90°  of  microphone  #  2  and  #  4  should  be  converted 
to  negative  voltages  while  90°  through  180°  of  microphone  # 2  and  #  4  should  be 
converted  to  positive  voltages  for  the  same  reason;  90"  in  all  microphones  should  be 
converted  to  zerp  voltage.  To  avoid  the  ambiguity,  in  the  real  implementations  the 
signals  fror>"  microphone  #  1  and  #  3  are  inverted  (i.e..  a  phase  change  of  7t)  so  that 
the  phase  range  of  180"  -  0°  is  changed  to  0°  — 180".  In  consequence.  0"  through  ±90" 
in  all  the  microphones  should  be  converted  to  negative  voltages  while  90"  through  180" 
or  -90°  through  -180°  arc  to  positive  voltages. 

The  phase-voltage  converter  with  the  characteristics  shown  in  Figure  5.20  can  meet 
the  above  requirements.  The  implementation  block  diagram  is  shown  in  Figure  B.l.  It 
is  important  to  remark  that  this  application  requires  the  PVC  to  be  insensitive  to  the 
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amplitude  of  the  input  signal  since  the  BPF  tone  fluctuates  severely  in  the  testing.  As 
it  can  be  seen,  the  BPF  tone  extracted  from  a  casing  microphone  and  the  rpm  signal 
are  converted  to  the  TTL  signal  iy.e..  high-level  of  5  volts  and  a  low-level  of  0  volts) 
before  they  are  input  to  the  exclusive-OR  gate.  The  frequency  multiplier  is  needed 
because  the  frequency  of  rpm  signal  from  the  magnetic  pickup  is  only  one  third  the 
blade  passing  frequency.  The  phase  shifter  is  needed  to  set  the  relative  phase  to  0°  at 
noblowing. 

The  exclusive  OR  gate  is  the  kernel  component  in  the  PYC.  Figure  B.2  depicts  the 
truth  table  that  defines  the  gate  operation.  The  output  goes  high  when  only  one  input 
goes  high,  and  the  gate  output  remains  at  a  logic  low  for  all  other  inputs.  Also,  note 
that  the  figure  depicts  v,  and  v0  as  having  the  same  frequency.  Thus,  the  mean  value  tv, 
obtained  using  a  low  pass  filter  (LPF)  following  XOR,  displays  a  triangular  function  of 
the  phase  difference  A d>.  Finally,  the  exact  plot  shown  in  Figure  5.20  can  be  obtained 
if  the  LPF  output  is  offset  by  a  DC  voltage. 
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Figure  B.l:  The  block  diagram  for  the  PYC  implementation. 


Figure  B.2:  The  operation  of  the  exclusive  OR  gate. 


Appendix  C 


Introduction  of  TEB  into  CUP2D 
Code 


This  appendix  first  describes  the  basic  idea  of  introducing  TEB  into  the  CUP2D  code. 
The  FFT  coefficients  of  the  TEB  model  is  then  derived  which  are  used  in  Chapter  6. 


The  basic  idea  behind  the  formulation  of  the  aeroacoustic  CUP2D  code  is  to  obtain 
the  stator  and  rotor  loadings  first,  then  to  compute  the  sound  field  according  to  them. 
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A  summary  of  the  theory  is  as  follows.  The  stator  loading.  LS.  and  rotor  loading.  LR. 
are  obtained  by  solving  Eqn.  (C.l).  It  should  be  pointed  out  that  each  variable  in  Eqn. 
(C.l)  denotes  a  matrix.  ITi?  and  ITS  mean  the  wash  velocity  on  the  rotor  and  stator. 
KRR  ( KSR ,  etc.)  matrix  is  called  influence  coefficient  matrix  which  relate  upwash 
velocities  at  the  stator  and  rotor  control  points  as  functions  of  the  unknown  loading 
on  the  rotor  and  stator  panels.  KRR  denotes  the  effect  of  rotor  on  rotor,  while  KSR 
means  the  effect  of  stator  on  rotor,  and  so  on. 


KRR 

KSR 

LR 

WR 

KRS 

KSS 

LS 

ITS 

In  addition,  two  boundary  conditions  should  be  satisfied  when  solving  Eqn. (C.l).  The 
flow  tangency  requirement  on  the  first  row,  say  if  it  is  stator,  is  shown  in  Eqn.  (C.2). 

lTS(n, i)  =  WSS(nJ)  +  WRS{n,i)  =0  for  all  n,  i  (c-2) 

where  loading  harmonics  are  counted  by  n  and  control  points  are  counted  by  i.  Eqn. 
(C.2)  states  that  the  sum  of  the  upwash  at  the  stator  due  to  the  stator  and  rotor, 
denoted  by  IT SS(n,  i)  and  I VRS(n,i),  respectively,  is  zero  for  all  harmonics  and  at  all 
control  points.  For  the  second  row  in  the  2D  model,  say  if  it  is  rotor,  the  condition  is 
different  since  it  is  excited  by  the  viscous  wake  of  the  upstream  stators,  as  indicated  in 
Figure  6.6.  If  we  call  the  wake  upwash  harmonics  HTTS(n,  i).  then  the  rotor  boundary 
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condition  is: 

WR{n.i)  =  WRR(n,  i)  +  WSR(n,  i)  =  -W\YS(n,j)  for  all  n.  i  (c-3) 

It  indicates  that  the  net  effect  of  the  waves  from  the  rotor,  \VRR{n,  i)  and  the  waves 
from  the  stator.  \VSR(n,  i)  must  be  equal  and  opposite  to  the  wake  upwash  for  flow 
tangency. 

The  active  TEB  model  should  be  transformed  into  frequency  domain  using  Fourier 
integral  as  required  by  CUP2D.  The  analytical  derivation  to  compute  the  Fourier  co¬ 
efficients  conducted  here  is  similar  to  the  work  by  Hanson  [48].  The  pure  wake  profile 
of  upstream  row  is  rewritten  here  according  to  Eqn.  (6.1): 

5^  =  H’dc(si)  •  f{y,  (C.4) 

where  tfit  streamwise  and  pitchwise  coordinates  si  and  y  are  shown  in  Figure  6.6; 
/(iMi)  is: 

f(y.si)  —  e_0'693['^]2  (C.5) 

In  Eqn.  (C.4)  and  (C.5),  Hdc(si)  and  <J(si)  are  governed  by  the  Lakshminarayana- 
Davino  model  which  is  represented  by  Eqn.  (6.6)  and  Eqn.  (6.7).  Reference  [48]  shows 
that  y  can  be  expressed  as: 


V  —  V'iCosOi  —  00362  •  (fli  —  )Rt  —  xsindi 


(C.6) 
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where  t  denotes  time.  With  the  use  of  Eqn.  (C.G)  and  relations  of  x  =  xs  +  s2cos02 
and  y2  =  s2sin02,  Eqn.  (C.5)  can  be  expressed  in  terms  of  t  as: 

'?in(02  -  0x)s2  —  xssin6i  —  cosf?i(^i  ~  Q2)RC2^ 


f(t,  si)  =  exp  -0.693 


<5(si) 


(c.?; 


Since  it  is  recognized  that  this  one  pulse  repeats  itself  at  the  blade  passing  period  of 
T  =  27r/£i(fii  -  fl2),  we  can  express  f(t.  Si)  in  a  Fourier  series  of  the  following  form: 

OO 

fits 0=  X]  Fn(Sl)  ■  einB^-^)l  (C.8) 


.  n—~oc 

The  usual  Fourier  series  formula  leads  to: 


2^/n  IVdc(si) 

F"(si)  =  -q  •  -i?rexp 


V. 


exp 


/ inBi 


—s2sin(9i  —  62)  +  xssin9i] 
cos0l 


(C.9) 


where. 


s/FG\cos9x 


2.13<5(.s1) 

Similarly,  the  Fourier  coefficients  associated  with  the  TEB  jet  can  be  obtained: 

'  inB\  [— s2sin{9]  —  92)  +  xssin9 


(C.10) 


n[])  QT  Wi  P 


where. 


_(™)2 

.  [QT  J 


exp 


cos6\ 


(C.ll) 

^/ttG  [COs6[ 

2.13a05(si) 

(C.12) 

M*(Sl)  =  bo-WM 

(C.13) 
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Now,  only  the  components  of  the  wake  and  TEB  jet  normal  to  the  stator  chord  are  re¬ 
quired.  This  is  given  by  multiplying  sin(91  —  #2)-  Hence,  the  desired  upwash  harmonics 
to  be  used  in  the  code  can  now  be  written  as: 

H^rma/j  =  H  dc(Sl)  COS02  ^  +  .  e<nBi(fl,-n2)0  (C>14) 

W2  J„  Ml  cosQi 

where  for  convenience,  the  relation  of  \\\/W2  =  cos02/ cos9\  is  used;  H  „orma/  is  the 
component  of  the  mixed  flow  normal  to  the  stator  chord.  The  role  of  TEB  is  reflected 
by  -Fj(si);  Fn{s i)  is  associated  with  the  pure  wake  as  shown  in  Eqn.  (C.9)  while  Fj 
can  be  obtained  from  Eqn.  (C.ll).  This  equation  should  be  input  into  the  subroutine 
GETWAKE  to  compute  WTTS(n,  i)  as  needed  for  inversion  of  the  matrix  equation  of 


(C.l). 


Appendix  D 


Flow  Condition  Estimation  for 
Simulator  Running  at  30,000  rpm 


This  appendix  derives  the  flow  condition  parameters  for  the  simulator  running  at  30,000 
rpm.  These  parameters  including  the  axial  flow  velocity,  sound  speed  and  fluid  density 
upstream  IGY,  in-between  IGY  and  rotor  and  downstream  rotor  are  required  to  input 
the  CUP2D  code  to  conduct  analytical  studies  in  Chapter  6. 
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Row  of  IG  V  Row  of  rotor 


Figure  D.l:  Re-depict  of  the  2D  cascade  model  of  the  simulator. 

The  2D  cascade  model  is  re-depicted  here  as  shown  in  Figure  D.l  for  convenience.  Three 
regions  in  the  model  are  defined  which  arc  upstream  IGY  (denoted  by  ‘‘a”),  in-between 
IGY  (“b")  and  downstream  rotor  (“c”).  The  axial  flow  mach  number,  sound  speed  and 
fluid  density  in  region  “a’:  are  denoted  as  Mxa,  Aa  and  pa,  respectively.  Similarly,  Mxb. 
At,,  pi,  and  Mxc,  Ac.  pc  represent  those  in  region  ‘b':  and  'lc“ ,  respectively;  Mc  is  the 
Mach  number  of  absolute  flow  velocity  in  region  V.  Since  the  non-turning  IGYs  are 
used  in  the  first  row.  the  flow  conditions  in  the  region  “a”  are  identical  to  those  in 
region  “b”  except  that  IGY  wakes  are  present  in  region  "b" .  The  one-stage  simulator 
is  estimated  to  provide  a  pressure  ratio  of  1.5  with  an  efficiency  of  80%;  Mxa  and  Mxc 
are  assumed  to  be  0.3  and  0.6,  respectively.  The  stagger  of  rotor  blades  is  assumed  to 
be  45°.  The  ambient  temperature  is  300  Kelvin. 
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Since  Mxa  is  known,  the  ratio  of  the  static  pressure  Pa  to  the  total  pressure  P0a  and 
the  ratio  of  the  static  temperature  Ta  to  the  total  temperature  T$  can  be  found  in  the 
well-established  fluid  property  look-up  table  [71]: 


pa 

—  =  0.94 

TDa 

M) 

'Ta 

—  =  0.98 
Ta 

J  o 


(D.l) 

(D.2) 


where  P0n  is  equal  to  the  ambient  pressure  which  is  101,000  Pa  and  P0o  is  300  K.  It 
follows  that  P“  =  94,940  Pa  according  to  Eqn.  (D.l)  and  Ta  =  294  K  according  to 
Eqn.  (D.2).  Using  these  results,  the  fluid  density  pa  can  be  obtained  according  to  the 
ideal  gas  law: 

pa 


Pa  = 


RTa 

=  1.125  kg/m 3 


(D.3) 

(D.4) 


wheie  R  is  the  gas  constant  which  has  the  value  of  287  N  ■  m./kg  ■  K.  Similarly,  the 
sound  speed  Aa  can  be  found  [71]: 


sJlRT0 

(D.5) 

343.7  m/s 

(D.6) 

where  7  is  the  specific  heat  capacity  ratio,  i.e.,  1.4.  Up  to  here,  the  flow  conditions 
in  region  ;‘a”  and  “b”  have  been  obtained  since  both  regions  share  the  same  flow 


parameters.  They  are 
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2^1  xa  —  ^fxb 

A  a  —  ^6 

/':  .  Pb 


0.3 

343.7  m/s 
1.125  kg/m 3 


Since  the  pressure  ratio,  Pq/Pq,  produced  by  the  rotor  is  1.5,  the  ideal  total  tempera¬ 
ture  in  the  region  “c” ,  T/'1  can  be  obtained,  if  the  process  is  isentropic  and  no  efficiency 
is  lost,  by  solving  [71] 


7—1 

/  pc\  7^,2 

[  r0  |  _  _  0 

\pbJ  n 


(D.7) 


Tq'  is  computed  to  be  336.8  K.  Since  the  efficiency  is  estimated  to  be  80%,  the  actual 
total  temperature  Tq  should  be  calculated  according  to  the  following  relation 


(Tq*  —  To)  =  (To  -  Tq)  ■  80% 


(D.8) 


It  follows  that  Tq  —  346  K.  To  get  the  static  temperature  Tc,  the  absolute  flow  velocity 
Mc  should  be  obtained  first.  Since  the  rotor  blade  stagger  02  is  assumed  to  be  45°  and 
Mxc  is  assumed  to  be  0.6.  Mc  is  calculated  to  be  0.84  according  to  the  velocity  triangle 
in  Figure  D.l.  With  the  value  of  A/c,  the  following  relations  can  be  found  again  in  the 


■*0 

'J'C 


0.62 

(D.9) 

0.87 

(D.10) 

look-up  table 
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where  the  total  pressure  P0C  has  a  value  of  152,000  Pa  since  the  pressure  ratio  is  1.5: 
the  total  temperature  T0C  is  346  K  as  found  from  Eqn.  (D.7).  It  follows  that  Pc  is  then 
computed  to  be  93,600  Pa  according  to  Eqn.  (D.9)  and  Tc  is  301  K  according  to  Eqn. 
(D.10).  The  fluid  density  pc  is  calculated  by  using  the  ideal  gas  law  once  again: 


Pc  — 


pc 

RTC 


(DTI) 


=  1.08  kg/m3 


(D.12) 


Similarly,  the  sound  speed  Ac  can  be  found: 


Ac  =  \fWc 

=  347.8  m/s 


Therefore,  the  flow  conditions  in  region  “c”  have  been  obtained: 


Mxc  =  0.6 
Ac  =  347.8  m/s 
pc  =  1.08  kg/m 3 


(D.13) 

(D.14) 


As  a  final  note,  the  drag  coefficient  of  the  IGV,  Cp  is  selected  as  0.01  in  CLP2D  code. 
This  value  is  obtained  by  using  Figure  9.24,  Page  600  of  the  reference  [55]  since  the 
Reynolds  number  calculated  for  the  case  is  around  7.0  •  10°. 


